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OPTOELECTRONIC  MATERIALS  AND  COMPONENTS :  MINIATURIZED 
and  thin  ft tji  laSERS  axd  modulator S 


SUMMARY 

The  research  being  carried  out  under  the  subject  contract  is  directed 
toward  the  development  of  miniaturized  lasers  and  modulators  compatible 
with  dielectric  waveguide  and  fiber  optic  transmission  elements.  Specific 
objectives  include: 

1.  development  of  rare  earth  doped  fluoride  ana  oxide  materials 

suitable  for  miniaturized  and  thin  film  lasers; 

2.  development  of  techniques  for  the  epitaxial  growth  of  KTN 

compositions  in  structural  configurations  appropriate  for 
electro-optical  light  modulation; 

3.  development  of  light -emitting  dWes  as  pump  sources  for 

miniaturized  lasers; 

4.  design,  fabrication  and  evaluation  of  miniaturized  and  thin 

film  laser -source  "packages"  which  include  the  associated 
•  pump  and  modulator. 

The  following  listing  summarizes  our  achievements  during  the  first 
six  months  of  the  program: 

1.  We  have  demonstrated  laser  action  at  2  ym  in  a  single  mode 

Ho:YLF  minilaser  0.2  x  0.6  x  0.9  mm  in  dimension.  Lasing 
occurred  at  a  threshold  power  of  20  mw;  an  output  power  of 
*  j  mw  was  obtained  at  150  mw  pump  power.  The  laser  was 

operated  in  a  dielectric-cavity  mode  which  is  particularly 
well  adapted  for  coupling  to  optical  waveguides  (Section  II). 

2.  Films  of  KTN  have  been  successfully  grown  on  KTa03  substrates  by 
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carried  out  (Section  V) . 

We  expect  by  the  end  of  the  first  year  of  the  program  to  have  an 
LED  pumped  2  pm  Ho:YLF  laser  operating  cw  at,  or  near,  room  temperature. 

We  believe  that  in  this  minilaser  we  will  have  a  prototype  of  a  class  of 
lasers  that  will  be  extremely  useful  for  integrated  optical  systems.  The 
only  serious  competitive  source  on  the  scene  is  the  diode  laser,  and  its 
future,  as  a  room-temperature  source,  appears  clouded  by  problems  involving 
reliability  and  lifetime. 


I.  INTRODUCTION 


When  we  began  this  program  we  believed  that  we  had  sufficient  know¬ 
ledge  of  energy  transfer  and  fluorescence  in  lithium  yttrium  fluoride  (YLF) 
to  allow  us  to  predict  with  some  confidence  that  the  development  of  a 
miniaturized  2  pti  Ho:YLF  laser  was  a  realistic  prospect.  Our  wo^  to  date, 
described  in  Section  II,  bears  out  this  initial  confidence.  We  have 
demonstrated  lasing  action  in  a  crystal  of  volume  0.2  x  0.5  x  0.6  mm  The 
threshold  for  lasing  was  low,  about  20  mw,  and  we  were  able  to  get  a  cw 
output  of  about  5  mw  at  a  pump  power  of  150  mw.  This  laser  appears  to  be 
the  smallest  dielectric  laser  ever  to  operate  at  such  low  thresholds.  Our 
goal  is  to  achieve  lasing  action  in  samples  of  even  smaller  dimensions. 

The  laser  is  side  pumped  and  operates  in  an  internal  reflection,  two- 
dimensional,  dielectric  cavity  mode.  A  mode  analysis  is  given  in  the 
Appendix*.  Output  occur*,  via  yvaneSiuai  £ icxud  csioulouw 

outside  the  dielectric  boundaries  by  the  totally  reflected  internal  waves. 
This  electromagnetic  configuration  appears  particularly  attractive  for 
integrated  optics  because  it  allows  coupling  to  a  waveguide  structure 
merely  by  placing  the  cavity  in  close  proximity  to  the  guide.  Coupling 
structures  such  as  prisms,  wedges  or  gratings  are  not  required.  We  are 
continuing  this  work  with  emphasis  on  obtaining  a  Ho:YLF  composition 

operable  at  or  near  room  temperature  in  the  near  term,  and  are  investi- 

3+  34- 

gating  the  feasibility  of  utilizing  other  activators  such  as  Tm  and  Er 
in  YLF  as  well  as  improvirg  the  pumping  efficiency  and  other  factors  by 
going  to  a  high  refractive  index  oxide  host.  Further  miniaturization  is 
also  a  goal  in  order  to  match  the  size  of  the  diode  pumps. 

The  work  to  date  on  the  diode  pumps  is  described  in  Section  III. 

The  thrust  of  the  program  has  been  the  design  and  fabrication  of  specialized 
diode  geometries;  the  growth,  testing  and  evaluation  of  Si-doped  SaAs 
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LED's  with  emission  tailored  to  the  Yb:YLF  minilaser  absorption  band  and 
with  efficient  high  intensity  output;  and  the  growth  and  characterization 
of  Zn  and  Si  doped  GaAsSb  LED's  with  the  objectives  of  understanding  the 
system  for  use  vs  more  sophisticated  heterostructure  diode  geometries. 

The  technology  development  to  permit  the  growth  of  desired  geometries 
in  the  GaAstSi  system  is  virtually  complete  and  several  runs  have  been 
made  to  explore  the  emission  wavelength  and  bandwidth  as  a  function  of 
various  parameters.  It  is  expected  that  useful  pump  diodes  from  this 
system  will  be  available  befcre  the  end  of  the  next  reporting  period. 

The  growth  of  thin  KTN  films  for  electro-optical  waveguide  modulators 
is  described  in  Section  IV.  Films  have  been  grown  on  KTaO^  substrates  by 
a  liquid  phase  epitaxy  technique  and  have  shown  waveguiding  in  preliminary 
experiments.  Masks  are  being  fabricated  for  appropriate  electrodes  to 
permit  evaluation  of  a  wave  deflection-modulation  scheme  as  well  as  a 
simple  step  phase  modulator.  Chemical  vapor  deposition  growth  experiments 
have  been  made  and  a  method  for  generating  and  transporting  KC1  vapor  has 
been  developed.  An  improved  two-zone  furnace  is  being  constructed  and  we 
expect  by  the  end  of  the  next  reporting  period  to  be  able  to  evaluate  this 
technique  for  tl  a  growth  of  KTN. 

Instrumentation  for  accurate  measurement  of  refractive  indices  of 
crystals  is  described  in  Section  V,  together  with  a  description  of  the 
newly  installed  clean  air  chemistry  laboratory. 

A  mode  calculation  for  the  dielectric  resonator  minilaser  is  given 
in  the  Appendix. 
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II.  MINILASEd  DEVELOPMENT 
2.1  Prototype  Mlnilaser 

To  obtain  an  efficient  minilaser  it  is  necessary  for  the  laser 

crystal  to  have  a  high  absorption  coefficient  at  the  pump  wavelength. 

This  objective  can  be  achieved  by  employing  a  high  concentration  of  the 

active  ion  as  in  the  NdLa  pentaphosphate  laser*-  or  by  incorporating  into  the 

host  a  high  concentration  of  sensitizer  ions.  The  latter  approach,  which 

we  are  following  in  our  work,  has  proved  successful  in  the  2.06  lim  HosYLT 

2 

laser  sensitized  with  Er  and  Tm. 

Additional  problems  encountered  in  the  minilaser  involve  the  design 

of  the  resonator  and  output  coupling  scheme.  One  resonator  configuration 

3 

that  has  been  employed  is  simply  a  scaled  down  version  of  the  usual  rod 
laser  with  reflecting  end  plates.  A  laser  of  this  type  would  have  to  be 

more  attractive  route  is  to  ure  a  configuration  capable  of  sustaining 

resonant  modes  confined  to  the  cavity  volume  by  total  internal  reflection 

(Fig.  1).  The  evanescent  field  associated  with  total  internal  reflection 

can  be  used  to  couple  to  an  output  structure.  Using  such  a  scheme  we 

have  succeeded  in  obtaining  single  mode  cw  laser  action  at  2.06  ym  in  a 

single  crystal  of  LiY^  462ErQ  s^o  035H°0  003*4'  The  laser  transition  is 

3+5  5 

from  the  lowest  level  of  Ho  1^  to  the  upper  Ig  level. 

The  experimental  details  are  shown  ir.  Fig,  2.  The  laser  crystal  with 
Its  four  long  sides  polished  and  with  dimensions  0.94  mm  x  0.60  mm  and  5  mm 
long  was  mounted  over  an  undoped  LiYF^  prism,  A  tapered  gap  was  formed  by 
resting  one  end  of  the  crystal  on  top  of  a  20  micron  thick  evaporated  indium 
pad  while  the  other  end  rested  on  the  prism.  A  0.2  mm  wide  region  of  the 

O 

crystal  was  pumped  with  the  output  of  a  cw  krypton  ion  laser  (6471  A).  By 
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Fig.  1  Internal  refl  ction  mode  in  cavity  coupled  to  dielectric 

waveguide  via  evanescent  field  in  coupling  gap. 

I  ' 
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Internal  reflection  modes  in  k-cpace  for  the.  F.o:\’LF 
oinilaser.  Total  internal  reflection  occurs  for 
those,  modes  within  the  sector  60.  The  mode  marked 
0  is  the  one  appearing  at  threshold;  the  nearby 
modes  marked  x  appear  above  threshold  (see  Fig. 4). 


varying  the  position  of  the  pumped  region  along  the  length  of  the  crystal, 
the  output  could  be  observed  as  a  function  of  coupling  strength.  The 
laser  crystal  and  prism  were  mounted  in  vacuum  in  a  liquid  dewar. 

The  resonator  ir  closed  on  four  sides  by  total  internal  reflecting 
dielectric  boundaries  and  open  on  the  other  two  sides.  Since  only  a 
small  segment  of  the  active  material  is  pumped,  only  those  modes  with 
propagation  vector  transverse  to  the  long  axis  of  the  crystal  are  excited, 
j't  is  shown  in  the  Appendix  that,  in  a  rectangular  crystal,  the  condition 
for  total  internal  reflection  modes  to  exist  is  n^  >  /  2  where  n^  is  the 
refractive  index  of  the  crystal  and  n2  is  the  index  of  the  surrounding 
medium.  Figure  3  shows  the  region  of  k^-space  available  for  internal 


reflection  modes  in  a  rectangular  crystal  wi  .h  n^  =  1.46  n2»  a  situation 
which  applied  tc  a  crystal  of  LiYF^  (n, )  in  vacuum  (.n^)  •  The  angle  60 
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50  *=  2  tan'1  l  A)  -  1]*-  y  ,  (1) 

2  * 

and  for  our  experimental  situation  is  numerically  equal  to  3.5°. 

Light  is  coupled  out  of  the  resonator  into  the  prism  through  i/'e 
evanescent  field  in  the  gap. 

Output  spectra  of  the  laser  are  shown  in  Fig. 4.  Threshold  was 
reached  with  approximately  20  mW  pump  power.  Near  thr behold  laser  action 
took  place  in  only  one  mode,  as  is  evident  from  the  sirgle  line  in  the 
spectrum.  The  linewidth  is  determined  by  the  resolution  limit  of  our 
spectrometer.  With  pump  power  about  twice  threshold,  two  other  modes 
appear.  T1  ;  unequal  spacing  of  the  modes  can  be  explained  by  considering 
the  modes  of  a  two-dimensiona-L  resonator  as  is  done  in  the  Appendix.  ’’’he 
theoretical  anaLysis  shows  that  «.  modes  which  appear  at  increasing  pump 


2.064  2.065  2.066 


2.064  2.065  2.066  2.067  2.068  2.069 

Wavelength  (/zm) 

Fig.  4  Output  spectra  of  Ho:YLF  minilaser:  (a)  near  threshold; 
(b)  twice  threshold. 
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power  are  separated  from  the  Initial  mode  by 

‘  7^  ^  +  &  m 

where  Ap,  Aq  *  0,  ±1,  ±2,  ±3,  etc.  The  modes  appearing  in  Fig. 4  corres- 
pond  to  AXq  (Ap  =  0,  Aq  =  1)  =  -1.1  x  10  pm  and  AXq  (Ap  =  -1,  Aq  =  1)  = 

_3 

0.62  x  10  pm.  The  experimental  values  are  in  precise  agreement  with 
theoretical  predictions.  The  relative  positions  of  these  modes  in  k-space 
are  shown  in  the  inset  in  Fig.  3. 

Output  powers  of  5  mW-  were  measured  with  150  mW  pump  power  and  a  gap 
size  of  5  microns .  The  threshold  did  not  increase  monotonically  with 
decreasing  gap  size,  indicating  that  losses  associated  with  edge  chipping 
and  other  crystal  defects  partially  determined  the  threshold.  Lower 
thresholds  are  expected  with  more  careful  sample  preparation. 

The  transmission  through  the  coupling  gap  is  largest  for  modes  whose 
propagation  vector  is  incident  on  the  base  of  the  crystal  at  the  critical 
angle,  and  decreases  rapidly  with  increasing  incident  angle.  Figure  5 
shows  the  calculated  transmission  loss  versus  mode  angle  for  a  TE  cavity 
mode  coupled  across  a  5  pm  gap.  Only  the  lowest  loss  modes  will  reach 
laser  threshold.  Note  that  internal  reflection  modes  exist  only  in  the 
interval  between  43.3°  and  46.7°.  When  the  mode  angle  increases  above 
43.3°  (the  critical  angle)  we  enter  the  region  of  internal  reflection,  and 
transmission  of  energy  out  of  the  cavity  modes  occurs  via  the  evanescent 
field  at  the  output  coupling  face.  In  our  experiments  this  was  the 
narrow  face  of  the  crystal  (dimension  a).  With  increasing  angle  there  is 
progressively  less  transmission  across  the  gap.  At  40.7°  evanescent  fields 
still  exist  at  the  coupling  face,  where  the  condition  for  internal  reflection 
is  still  maintained.  However,  we  no  longer  satisfy  this  condition  on  the 
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Fig.  5  C, 
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orthogonal  b  face  and,  consequently,  the  transmission  loss  rapidly 
increases  as  energy  is  radiated  from  the  two  b  faces. 

We  see  in  Fig.  3  that  the  cavity  modes  which  are  excited  above  thres¬ 
hold  have  an  angle  of  incidence  9  (measured  from  the  normal  to  the  a  face) 
which  is  less  than  the  angle  for  the  mode  with  lowest  threshold.  This 
result  is  consistent  with  loss  curve  shown  in  Fig.  5, 

As  a  source  for  integrated  optical  networks,  this  type  of  minilaser 
has  the  advantage  of  much  higher  efficiency  than  proposed  thin  film  rare 

4 

earth  lasers.  The  largest  absorption  coefficient  of  a  rare  earth  ion 
transition  that  can  be  pumped  by  a  light-emitting  diode  is  about  a  = 

200  cm  l.  Thus,  for  a  thin  film  1  micron  thick,  only  27.  of  the  incident 
light  can  be  absorbed  when  transversely  pumped.  With  the  dielectric 
resonator,  the  dimension  in  the  pumping  direction  can  be  made  about  1/a, 
so  that  most  of  the  pump  light  wiJ 1  be  absorbed,  ’fading  Lv  au  increase  le 
pumping  efficiency  by  a  factor  of  about  50. 

The  number  of  possible  resonant  modes  can  always  be  made  arbitrarily 
small  in  this  scheme  by  proper  control  of  the  indices  n.  and  n2> 

Decreasing  the  size  of  the  resonator  increases  the  mode  separation  in 
k-space.  Thus,  with  smaller  resonators  single  mode  operation  should  be 
obtainable  even  at  pump  powers  many  times  threshold. 

The  Ho: alphabet  YLF  laser  materials  studied  so  far  have  run  cw  only 
at  liquid  N2  temperature.  Future  work  will  be  directed  toward  a  room 
temperature  miniature  cw  holmium  laser.  To  obtain  room  temperature  cw 

'  3+ 

operation  the  Ho  concentration  must  be  low  (<  0.5%)  so  that  the  terminal 

laser  level  will  have  a  small  population,  but  there  are  two  factors 

3+ 

limiting  how  low  the  Ho  concentration  can  be  made:  (1)  res  ator  loss, 
and  (2)  energy  transfer  rates. 
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If  laser  resonators  with  losses  of  0.01%  can  be  made,  materials  with 
3+ 

less  than  0.01%  Ho  would  have  sufficient  gain  for  cw  laser  operation. 
Conventional  Fabry-Perot  resonators  have  a  mininun  loss  of  about  1%  from 
mirror  scattering  and  diffraction  losses.  The  miniature  total  internal 
reflection  resonator  should  have  much  less  loss.  Future  work  will  be  done 
on  measuring  and  reducing  the  losses  of  this  resonator.  Edge  chipping  is 
a  problem  in  fabricating  these  small  crystals;  scattering  losses  from 
these  defects  will  be  less  from  resonators  of  different  geometries  than 
those  previously  used. 

In  the  Ho-alphabet  YLF  laser  material  energy  transfer  between  the 
3  5 

Tm:  and  the  Ho:  1^  levels  is  very  fast  compared  to  the  fluorescent 

lifetimes  of  the  levels  so  that  the  populations  of  the  levels  can  be 

3+ 

described  by  a  thermal  equilibrium.  Therefore,  a  low  Ho  concentration 
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Ho  fluorescence  will  remain  high.  It  is  not  yet  known  how  low  the  Tm 

3+  3+ 

and  Ho  concentration  can  be  made  before  energy  transfer  from  Er  is  - 

inefficient . 

Single  mode  laser  operation  is  needed  for  integrated  optics 
applications.  The  total  internal  reflection  resonator  has  been  operated 
single-mode  for  pump  powers  up  to  twice  threshold.  Further  work  will  be 
done  to  obtain  single  mode  operation  over  a  larger  range  of  pump  powers  by 
making  the  resonator  smaller,  which  will  also  make  it  more  compatible  with 
the  LED  pump. 


2.2  Fluoride  Crystal  Growth 

Alphabet  YLF  samples  of  the  following  compositions  have  been  utilized 


for  this  program: 
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#159f 

L1Y0 . 48  7Er 0 . 5Tm0 . 01Ho0 . 003F4 

#171f 

LiY0 . 4  62ErO .  S^O .  01H°0 . 003F4 

#193f 

LiY0 . 628Er0 . 2OYb0  .lTta0 . 067Hoi 

#213f 

LiYb0 . 495Er0 . 4945Ho0 . OOD5TmO 

#2l7f 

LiErF.  :  0.005  Tm 

4 

#219f 

LiErF.  :  0.01  Tm 

4 

Sample  #171f  was  used  for  the  minilaser  described  in  Section  2.1,  while 

sample  #213f  has  yet  to  be  tested.  Measurements  on  samples  #2l7f  and 

219f  indicated  that  the  Tm  fluorescence  was  almost  completely  quenched  by 
3  4 

back  transfer  (Tm  Er  *^3/2^  ant*  subsequent  impurity  quenching  of  the 

A 

Er  *13/2  via  energy  migration  to  an  impurity  center.  The  impurity  center 
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3+ 

of  Dy  .  It  is  evident  that  higher  purity  ErF^  than  is  currently 
commercially  available  will  have  to  be  used  when  low  activator  concen¬ 
tration?  are  necessary.  Zone  refining  experiments  are  in  progress. 


2.3  New  Materials 

A  considerable  amount  of  exploratory  research  has  been  devoted  to 
finding  suitable  oxide  rare  earth  hosts  with  high  refractive  indices  to 
provide  a  better  opt  .1  match  to  the  pump  LED  and  to  give  more  flexibility 
for  coupling  into  modulator  and  .vave  guide  structures.  These  have 
included  CAMGAR  (CaY2Mg2Ge3012) ,  =  1.83,  and  Y^i^,  *  2.35.  The 

latter  material  is  of  the  most  interest,  not  only  for  the  higher  index  but 
also  because  this  pyrochlore  structure  places  the  rare  earth  ion  in  a  site 
with  a  center  of  inversion  symmetry  (3m  D^)  .  This  results  in  only 
magnrtic  dipole  transitions  being  allowed  and  may  provide  better  branching 
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ratios  for  some  transitions. 


q  i 

Y^2Ti2°7  cr>’stals  with  Er  compositions  of  0.1,  1  and  2  atomic  % 
have  been  grown.  The  refractive  index  as  a  function  of  wavelength  is 
shown  in  Section  V,  Fig.  19.  As  with  Tm3+  in  LiErF^  the  Er3+  fluorescence 
in  the  0.1%  sample  was  completely  quenched  by  impurities  and  much  purer 
feed  materials  will  have  to  be  used.  We  propose  to  use  solvent  extraction 
as  well  as  conversion  to  chlorides  for  zone  refining  as  methods  of  puri¬ 
fication.  The  new  clean  air  chemistry  laboratory  described  in  Section  V 
makes  these  procedures  practical. 
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III.  DIODE.  PUMPS 


3.1  Diode  Pump  Requirement 

The  solid  state  minilaser  excitation  source,  or  pump,  will  be  a  semi¬ 
conductor  light-emitting  diode  which  has  been  designed  specifically  to  pump 
the  laser  geometries  already  discussed.  The  diodes  to  be  used  will  be 
designed  to  have  planar  emitting  surfaces  to  which  the  lasers  can  be 
intimately  attached.  Two  possible  structures  are  illustrated  in  Fig.  6. 

The  diode  geometry  will  further  be  such  that  the  junction  can  be  well  heat 

2 

sunk  to  permit  operation  at  high  current  levels,  3-400  A/cm  ,  and  the 
junction  area  will  be  slightly  greater  than  that  of  the  laser  to  maximize 
the  pump  energy  entering  the  minilaser.  General  diode  structures  have  been 
designed  and  fabrication  techniques  are  being  developed. 

The  diode  emission,  as  well  as  geometry,  is  being  tailored  for 

*'  —  *•  1  ^  '  -«.i— wl  n«p  J  mV  nwQ  n  £  -i  n  n  fNwl  An  ^  A  f  p#>»iv*pn 
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but  the  emission  spectrum  too  must  be  adjusted  to  match  the  pump  band  of 
the  minilaser.  This  requires  as  narrow  an  emission  band  as  possible  peaking 
at  the  laser  absorption  maximum.  To  achieve  this  two  approaches  are  being 
followed:  Cl)  the  use  of  silicon-doped  gallium  arsenide  diodes,  and 

(2)  the  use  of  silicon  and/or  zinc-doped  gallium  arsenide  antimonide  diodes. 
Both  systems,  as  will  be  discussed  in  detail  in  the  following  section,  emit 
further  into  the  infrared  than  do  conventional  zinc-doped  GaAs  LED's,  which 
is  what  is  required,  and  each  has  its  own  particular  advantages.  The  GaAs;Si 


system  has  been  the  subject  of  considerable  study  so  that  our  work  here 
is  as  much  development  as  research.  We  have  been  concerned  in  this  program 
with  developing  our  own  capability  in  growing  the  required  diode  layers, 
with  tailoring  the  pump  diode  emission  spectrum,  and  with  maximizing  the 

emission  intensity  and  overall  pump  efficiency.  The  GaAsSb  system,  on  the 
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other  hand,  is  much  more  poorly  understood**  and  the  work  on  this  system  is 
considerably  more  basic,  starting  with  the  determination  of  parameters  in 
the  growth  system  and  the  behavior  of  dopants  in  the  alloy  system. 

.  ■»  summarize  then,  the  diode  minilaser  pumps  program  has  three  areas 
of  emphasis,  or  thrusts:  (1)  the  design  and  fabrication  (including 
technology  development)  of  specialized  diode  geometries;  (2)  the  growth, 
testing,  and  evaluation  of  Si-doped  GaAs  LED's  with  emission  tailored  to 
pump  the  Yb:YLF  minilaser  absorption  band  and  with  efficient,  high  intensity 
output;  and  (3)  the  growth  and  characterization  of  Zn  and  Si  doped  GaAsSb 
LED's  with  the  objectives  of  fully  understanding  this  system  and  evaluating 
it  for  ultimate  use  in  the  pump  diodes.  Initial  pump  diodes  and  minilaser- 


diode  structures  will  evolve  from  a  combination  of  the  results  of  programs 
1  and  2  above.  Future  structures  will  involve  more  sophisticated,  hetero- 


\t  materials  developed  in  the  third  effort. 


3.2  Emission  Tailoring 

The  minilaser  we  intend  to  build  has  a  strong  pump  band  at  0.97  pm 
and  it  is  this  band  that  the  pump  diode  emission  will  be  matched  to. 

Because  conventional  GaAs  diodes  emit  at  0.90  pm*^  this  requires  the  use  of 
either  silicon-doped  GaAs  diodes,  which  with  sufficient  silicon-doping  will 
emit  at  0.97  pm,  or  the  use  of  an  alloy  compound  system  such  as  GaAsSb  with 
silicon  doping  in  homostructures  or  with  zinc  doping  in  heterostructures, 
i.e.,  between  GaAs  layers. 

Silicon  doping  has  attractive  features  for  the  pump  diode  application. 
First,  silicon  can  enter  the  GaAs  lattice  either  at  gallium  sites  where  it 
acts  as  a  donor  or  at  arsenic  sites  where  it  behaves  as  an  acceptor.^ 

Which  position  it  enters  depends  on  the  temperature  the  crystal  is  grown  at, 
and  by  growing  an  epitaxial  layer  through  a  transition  temperature  (,v  900°C 
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There  is  some  evidence  in  the  literature  that  this  is  the  case,  but 
other  factors  such  as  the  width  of  the  n  to  p  transition  region,  that  is, 
the  junction  grading  constant,  must  also  be  considered.  One  of  our  main 
objectives,  ir.  fact,  is  to  determine  the  effect  of  the  grading  constant  on 
the  emission  spectrum  and  diode  efficiency.  It  is  definitely  the  case 
that  the  Si-doped  GaAs  emission  peak  is  broader  than  the  emission  peak  of 
Zn-doped  GaAs,  roughly  40  nm  vs  25  nm.  To  get  the  desired  0.97  pm 
emission  using  Zn-doping  or  lower  levels  of  Si-doping,  however,  a  lower 
bandgap  semiconductor  must  be  used,  and  for  this  purpose  the  alloy  GaAsSb 
is  being  studied.  The  room  temperature  bandgap  of  GaAs  Sb..  decreases 
from  1.425  eV  in  GaAs  to  0.725  eV  in  GaSb^,  and  for  x  0.97,  for  example, 
should  be  sufficiently  lowered  to  give  the  desired  emission  for  a  low  0.5% 
Si  doping.  In  homostructure  diodes  there  will  clearly  be  a  compromise 

»omn'ro/1  K !»f-T.roeir.  t-Vio  i.ii d tV  rtf  *-Vio  om-fcc-ipo  Hnp  and  fHo  leirpl  rtf  ciliVnn 

doping  needed  to  get  emission  sufficiently  far  below  the  absorption  edge  so 
that  reabsorption  is  not  a  problem.  If  heterostructures  are  used,  the 
reabsorption  problem  can  be  minimized  by  surrounding  the  emitting  GaAsSb 
layer  by  wider  bandgap  GaAs  layers,  and  that  problem  can  be  effectively 
separated  from  the  problem  of  the  emission  linewidth. 

3 . 3  Liquid  Phase  Epitaxy  of  Diode  Layers 

The  pump  diodes  are  fabricated  in  epitaxial  layers  grown  by  liquid 
phase  epitaxy,  LPE.  A  typical  device  might  require  a  p-type  layer  on  an 
n-type  layer  all  on  an  n-type  or  semi-insulating  GaAs  substrate.  As  was 
described  above,  when  silicon  is  used  as  a  dopant  only  one  growth  is 
necessary  to  produce  both  of  the  layers.  If  zinc,  a  p-type  dopant,  is 
used,  however,  an  n-laver  may  have  to  be  grown  first  from,  for  example,  a 
tin-doped  melt ,  and  then  the  p-type  layer  grown  from  a  second  zinc-doped 
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melt.  Multiple  growths  would  also  be  required  if  more  complicated  hetero¬ 
structure  devices  are  being  fabricated. 

We  currently  have  two  LPE  growth  systems  dedicated  to  this  program. 
Both  are  hor_zontal,  sliding  furnaces  similar  to  the  one  illustrated  in 
Fig.  7.  In  such  a  system,  the  substrate  and  melt  are  contained  in  a 
graphite  boat  like  that  in  Fig.  3.  As  can  be  seen  in  the  cutaway  sketch 
of  another  boat  in  Fig.  9,  the  substrate  rests  in  a  depression  in  the 
bottom  of  the  boat.  The  melt  consisting  of  gallium  saturated  (at  the 
growth  initiation  temperature)  with  arsenic  (and  antimony  when  growing 
GaAsSb)  is  contained  in  the  well  in  the  graphita  slider.  The  slider  illus¬ 
trated  in  Fig.  8  has  several  wells  and  thus  can  be  used  to  grow  multiple 
layers;  the  boat  in  Fig.  9  has  only  one  well. 

To  see  how  these  boats  are  used  in  an  LPE  growth,  it  is  illustrative 
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substrate.  Careful  preparation  of  the  substrate  is  one  of  the  singularly 
most  important  parts  of  a  successful  epitaxial  .growth.  High  quality 
material  must  be  used,  although  grown  layers  tend  to  have  higher  crystalline 
perfection  than  the  substrate,  and  the  substrate  surface  must  be  highly 
polished  and  chemically  etched.  Our  substrates  are  given  a  five  minute 
etch  in  5:1:1  etch  (five  parts  by  volume  H2S0^,  one  part  H2O2,  and  one 
part  H2O)  rinsed,  dried,  and  placed  in  the  boat.  The  boat  at  this  stage 
contains  the  gallium  for  the  melt.  The  gallium  is  loaded  into  the  boat 
first  where  it  is  baked  at  1000°C  in  hydrogen  for  twelve  hours  to  deoxidize 
it  prior  to  loading  the  other  items.  The  arsenic,  antimony,  and  dopants 
are  added  after  the  substrate  and  then  the  boat  is  sealed  in  the  furnace. 
Growth  is  performed  in  a  hydrogen  atmosphere.  With  the  boat  loaded  and  in 

position  the  system  is  flushed  with  hydrog. n  overnight.  It  is  then  hcatud 
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to  approximately  20° C  above  the  growth  initiation  temperature  and  allowed 
*o  equilibrate.  The  temperature  is  then  lowered  to  the  growth  initiation 
temperature  and  once  it  stabilizes  the  slider  is  moved,  using  a  quartz 
push  rod,  to  bring  the  melt  inco  contact  with  the  substrate.  It  is  often 
desirable  to  etch  cff  a  thin  layer  of  the  substrate  prior  to  initiating 
the  growth  and  this  can  be  accomplished  by  either  making  the  melt  purposely 
undersaturated  at  the  initiation  temperature  or  by  increasing  the  tempera¬ 
ture  a  few  degrees  after  the  sliding.  The  later  technique  affords  very 
good  control  over  the  etch  back  particularly  if  a  second  "dummy"  suv  trate 
is  used  to  equilibrate  the  melt  at  the  growth  initiation  temperature  prior 
to  moving  it  onto  the  "good"  substrate.  In  either  case,  the  temperature 
is  then  lowered  at  a  constant  rate,  25° C/hour  in  our  typical  growths,  until 
the  desired  layer  thickness  is  grown.  This  is,  of  course,  readily  deter- 
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literature. 


Once  the  proper  thickness  layer  has  been  grown,  the  melt  is  removed 
from  the  substrate  by  moving  the  slider  ove'  again  with  the  quartz  push 
rod.  The  furnace  is  then  cooled  to  room  temperature  and  the  substrate  is 
removed  and  used.  The  growth  of  multiple  layers  is  done  in  a  similar 
manner  although  etch  back  of  the  grown  layers  is  not  necessary,  nor  in  fact 
desirable,  and  the  use  of  a  dummy  substrate  is  essential. 

The  materials  used  in  the  growth  are  all  obtained  commercially. 

GaAs  substrates  are  purchased  from  Electronic  Materials  Corporation  in  the 
form  of  (100)  oriented  wafers,  0.015"  thick,  with  a  chemically  polished 
face.  These  are  cut  into  substrates  approximately  7  mm  square  on  a  wire 
saw  in  our  laboratory.  The  gallium  used  is  nominal  6-9' s  purity  and  the 
antimony  and  arsenic  arc  added  to  the  L : s  in  t:.-.-  form  of  undoped  poly¬ 
crystal  line  C  As  and  GaSb .  The  silicon  is  obtained  in  the  form  of  standard 
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arsenic-doped  waters  from  Monsanto  that  are  scribed  and  broken  into  small 
dice  in  the  MIT  Microelectronics  Laboratory. 

As  described  the  LPE  technique  appears  to  be  a  straightforward 
process,  as  it  should  be.  It  is  our  experience  that  with  certain 
precautions,  essentially  by  being  careful  at  the  proper  times,  the  commonly 
reported  difficulties  with  LPE  such  as  constitutional  supercooling  are  not 
in  fact  encountered.  The  key  areas  of  corcern  appear  to  be  (1)  that  the 
gallium  be  prebaked  and  the  system  be  carefully  flushed  prio?  co  growth  to 
avoid  oxides,  (2)  that  the  substrate  be  carefully  prepared  and  etched  to 
remove  surface  damage  prior  to  growth,  and  (3)  that  the  melt  be  removed 
completely  and  cleanly  after  growth.  The  latter  is  accomplished  by  cutting 
the  substrate  to  completely  fill  the  boat  substrate  well  and  by  keeping  the 
substrate  and  layer  thickness  such  that  there  is  a  small  clearance  between 
the  slider  bottom  and  the  grown  surface. 

A  modification  of  our  standard  technique  described  above  which  we  have 
found  significant  for  the  GaAsSb  growth  is  the  use  of  a  graphite  plug  in 
the  melt  well.  This  plug  keeps  the  celt,  which  in  the  GaAsSb  growths  con¬ 
tains  a  xarge  amount  of  Sb  and  has  a  large  surface  tension,  from  tailing  up 
and  improves  the  surface  coverage  and  layer  uniformity.  Most  importantly, 
it  reduces  the  crystallization  on  the  upper  surface  of  the  melt  which,  if 
not  controlled,  can  interfere  with  the  slider  motion. 

3.4  Results  to  Date 

In  the  initial  six  months  of  this  program  we  have  been  very  successful 
in  getting  two  LPE  grower.  systems  functioning  which  are  currently  producing 
both  epitaxial  GaAs:Si  and  GaAsSb’.Si  layers.  The  grown  samples  are 
divided  into  several  pieces,  one  of  which  is  angle-lapped  and  stained  to 
check  junction  depth  uniformity  and  junction  planarity.  Small  gold  contacts 
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are  evaporated  through  a  metal  mask  on  another  piece  and  it  is  cut  into 
test  diode  dice  which  are  mounted  as  illustrated  in  Fig.  10.  Spectral 
emission  and  efficiency  as  a  function  of  operating  current  are  then  obtained. 
A  typical  emission  spectrum  for  a  GaAs  device  with  1.7%  Si  doping  is  shown 
in  Fig.  11. 

In  the  GaAsSb:Si  program  we  have,  with  the  system  modification  des¬ 
cribed  earlier,  been  able  to  grow  high  quality  appearing  layers  with  good 
surface  coverage  and  flat  surfaces.  Me  are  currently  refining  this  tech¬ 
nique  to  improve  our  reproducibility  and  simultaneously  are  collecting  data 
on  the  n  to  p  transition  temperature  as  a  function  of  Si  and  Sb  concen¬ 
trations  . 

The  development  cf  pump  diode  geometries  and  fabrication  techniques 
has  also  proceeded  well.  A  David  Mann  photorepeater  has  recently  been 

-in  nnr  Mi  rmol  erfronirs  T.ahoratorv  and  will  be  used  to  Droduce  the 
necessary  mask  sets  for  the  pump  diode  patterns .  Actual  diode  fabrication 
must  await  growth  of  the  LPE  layers  but  work  on  development  of  fabrication 
techniques  has  proce  Jed  using  bulk  GaAs  samples,  i.e.,  bare  substrates. 
Techniques  for  etching  mesas  using  patterned  pyrolytic  SiO^  as  a  mask  and 
for  patterning  gold  contact  patterns  using  standard  photoresists  and  photo¬ 
lithographic  methods  have  been  developed.  We  are  currently  developing 
techniques  to  do  the  deep  n+  diffusions  in  GaAs  needed  for  the  device 
structure  shown  in  Fig.  6. 


Ball-bonded 


Figure  10-  Artist's  sketch  of  mounted  diode  chip  for  use  in 
testing  the  diodes  for  emission  spectrum  and  effi¬ 
ciency. 


Typical  diode  emission  spectrum 
doped  GaAs  diode  containing  1.7) 
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IV.  MODULATOR  PROGRAM 

We  describe  in  this  section  that  part  of  our  program  concerned  with 

the  growth  of  thin  films  for  electro-optical  waveguide  modulators.  We 

have  concentrated  on  materials  in  the  solid  solution  series  KTa,  Nb  0_. 

1-x  x  3 

Compositions  in  this  system  are  attractive  because  they  possess  a 
particularly  high  figure  of  merit  for  electro-optical  modulation.  Two 
distinct  growth  techniques  are  being  explored:  (I)  chemical  vapor  deposition 
(CVD) ;  and  (2)  liquid  phase  epitaxy  (LPE) . 

4.1  CVD  Growth 

A  CVD  system  for  the  growth  of  KTN  in  bulk  and  thin  film  form  has  been 
constructed.  The  syscem,  shown  in  Fig.  12,  has  provision  for  six 
independent  gas  flows.  Two  chlori.ie  flows  are  used  for  the  production  of 
TaCl,.  and  NbCl,.  vapor  by  the  controlled  chlorination  of  Ta  and  Nb  metal. 
Monitored  flows  are  provided  for  the  oxygen  and  hydrogen  necessary  for  the 
reaction 


2KCl+2(l-x)TaCl5+(2x)NbCl5+6H2+302  2KTa1_xNbx03+12HCl 

A  third  chlorine  flow  is  used  to  create,  between  the  chloride  streams  and 
the  H2  and  02  flows,  a  sheath  which  prevents  the  reaction  from  occurring  at 
the  nozzle,  where  it.  could  cause  clogging,  and  allows  the  reaction  instead 
to  occur  downstream  in  the  reaction  zone.  Finally,  helium  can  be  introduced 
to  suppress  back-diffusion  of  the  metal  vapors,  and  to  allow  for  purging 
the  system  during  start-up  and  termination  of  the  run. 

The  production  of  the  TaCl  and  NbCl  vapors  is  accomplished  without 
difficulty.  The  Ta  and  Nb  are  chlorinated  in  external  reactors  heated  to 
240°C  and  the  vapors  are  transported  to  the  furnace  in  heated  lines.  The 
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growth  of  KTN  by  chemical 
ration. 
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production  of  KC1  vapor  is  much  more  difficult.  After  a  number  of  false 
starts  we  ended  up  with  an  arrangement  in  which  liquid  KC1  is  vaporized 

l 

within  a  Knudsen  cell.  The  cell  is  a  cylindrical  quartz  vessel  3"  in 
length,  1''  in  diameter,  with  a  0.25  mm  diameter  hole  at  the  exit  end.  To 

obtain  a  controlled  flow  the  cell  is  positioned  in  the  furnace  at  a  point 
where  the  temperature  is  slightly  above  the  melting  point  of  KC1 
(ca.  770°C) .  Unfortunately,  in  our  present  single  zone  furnace  there  is 

a  substantial  gradient  at  the  point  where  we  must  position  the  Knudsen  cell. 

We  are  in  the  process  of  constructing  a  new  furnace  which  will  contain  two 
zones,  one  a  reaction  zone,  the  other  a  region  for  the  controlled  ■volatiliza¬ 
tion  of  KC1. 

Work  with  the  present  furnace  has  been  concerned  with  finding  the 
appropriate  conditions  for  KTN  growth.  Nineteen  runs,  lasting  from  1  to  8 
hours,  have  been  made  thus  far.  Runs  1  through  8  involved  general  nesting 
of  the  system  as  well  as  checks  on  TaCl^  and  NbCl^  production  and  flow. 

Runs  9  through  1^  were  concerned  with  finding  a  suitable  scheme  for  the 

a, 

production  of  a  controlled  flow  of  KC1  vapor.  Runs  15  and  16  were  tests 
to  see  whether  the  pentoxides  of  Ta  and  Nb  could  be  grown.  Runs  17  through 

19  were  unsuccessful  attempts  to  grow  KTN.  It  appears  that  if  we  are  to 

t 

be  successful  in  growing  KTN  by  chemical  vapor  deposition  we  need  greater 
flexibility  in  control  of  KC1  flow  than  we  have  in  our  present  furnace. 

Our  new  two-zone  furnace  should  provide  this  needed  flexibility. 

In  our  first  runs  with  the  new  furnace  we  plan  to  concentrate  on  the 
growth  of  KTaO^.  Only  after  this  composition  is  successfully  grown  shall 
we  proceed  with  the  growth  of  solid  solutions  in  the  KTN  system. 

4 . 2  Growth  by  Liquid  Phase  Epitaxy 

Using  the  method  of  liquid  phase  epitaxy  we  have  successfully  grown 
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thin  films  of  KTN  on  a  substrate  of  KTaO^.  We  are  employing  a  dip  method 
in  which  a  suitably  prepared  substrate  of  KTaO^  is  dipped  into  a  molten 
mixture  of  ^CO^,  Ta2°5  and  Ni32°5’  held  in  contact  witlt  the  melt  for  a 
brief  period,  and  upon  withdrawal,  is  rapidly  spun  to  remove  any  attached 
droplets . 

The  KTaO^  substrates  are  cut  from  crystals  grown  by  the  top-seeded 

solution  technique.  They  are  prepared  as  circular  discs  0,5  -  1.0  mm 

thick  and  8  -  18  mm  in  diameter  with  the  face  of  the  substrate  parallel  to 

a  (100)  crystallographic  plane.  Thet surface  on  which  we  desire  to  grow 

the  epitaxial  film  is  polished  on  a  GMI,  Model  1000  polishing  machine 

using  successively  finer  grades  of  diamond  paste  ending  with  0.25  um  grit 

size.  Flatness  is  maintained  to  within  one  quarter  wavelength  of  sodium 

light.  The  sample  is  degreased,  cleaned  in  a  warm  chromic-sulfur ic  acid 

SoIuL’.ou.^ai'd  .  c inset*  in  distilled  w<ftu. 

The  substrate,  attached  to  platinum  seed  rod  with  three  platinum  wire 

prongs,  is  slowly  lowered  to  within  6  cm  of  the  surface  of  the  melt  (melt 

temperature  1250°C)  where  it  is  allowed  to  stabilize  for  about  15  minutes 

and  then  lowered  into  the  melt  while  rotating  at  50  60  rpm.  After  a 

growth  period  varying  from  10  -  AO  minutes  the  sample  is  withdrawn  and  spun 

at  1000  rpm  foe  10  seconds.  The  seed  rod  is  then  lifted  at  a  rate  of  about 

2  cm /min  until  the  substrate  is  fully  removed  from  the  furnace. 

The  growth  rate  of  the  epitaxial  films  produced  by  this  procedure  is 

quite  rapid,  ca  5  -  10  ym/min,  suggesting  that  growth  occurs  under  super- 
* 

saturated  cond'  ions.  The  melt  used  for  the  runs  we  have  made  to  date 
Consists  of 

57.6  mole  percent  KjCOj 

25. A  mole  percent  NbjO. 

17.0  mcie  percent  TajO^ 
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The  "as  grown"  epitaxial  films  had  rippled  surfaces  and  when 
examined  in  polarized  light  showed  birefringence  despite  the  fact  that 
films  were  in  a  compositional  range  where  the  material  is  structurally 
cubic  at  room  temperature.  Presumably  the  birefringence  is  strain- 
induced  . 

The  compositional  variation  of  the  films  was  examined  with  an 
electron  microprobe.  Typical  results  are  shown  in  Figs.  13  and  14  for  a 
film  having  the  composition  KTa^  ggNbQ  The  scatter  in  the  data 

obtained  in  the  traverse  parallel  to  the  surface  lies  within  the  limits  of 
instrumental  accuracy.  The  degree  to  which  we  can  determine  the  abrupt¬ 
ness  of  the  transition  from  film  to  substrate  is  limited  by  the  positional 
accuracy  of  the  microprobe  system,  about  *  4  Jim. 

In  our  future  work  we  plan  several  modifications  in  film  preparation 
procedure.  In  order  *-c  slow  dw.-.  tne  grout*  rate:  *-»c  p'O’..  to  diivte  tne 
melt  somewhat  by  the  addition  of  a  suitable  flux.  We  will  carry  out  a 
systematic  investigation  of  the  way  in  which  film  growth  and  quality  is 
affected  by  the  various  steps  in  the  soak,  dip,  withdrawal  and  spin 
sequence.  To  reduce  film  strain  we  shall  explore  the  possibilities  of 
slightly  modifying  the  substrate  composition  in  an  attempt  to  improve  the 
lattice  constant  match  between  substrate  and  film.  We  also  hope  to  reduce 
strain  effects  by  ion  milling  the  substrate  subsequent  to  our  final 
diamond  polish. 

4.3  Film  Evaluation 

The  films  as  grown  were  not  suitable  for  light  guiding  because  of 
their  rippled  surfaces.  These  irregularities  were  eliminated  by  careful 
polishing.  The  resultant  films  were  optically  flat,  scratch-free,  and 
varied  in  thickness  from  50  -  70  pm.  Light  guiding  was. observed  using  a 
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Fig.  13  Microprobe  scan  of  KTN  film  on  KTaO^  substrate; 
Sv.a.i  is  parallel  to  surface. 
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Fig.  14  Microprobe  scan  of  KTN  film  on  KTaO^  substrat 
scan  is  perpendicular  to  film  edge. 


rutile  prism  to  couple  into  the  film.  By  using  a  double-sided  prism  we 

were  able  to  resolve  a  series  of  mode  pattern  "m-lines"  on  an  output 
12 

screen.  The  patterns  appeared  and  disappeared  at  approximately  the 
angular  patterns  predicted  by  our  theoretical  calculations.  We  are  now 
in  the  process  of  more  carefully  examining  the  guided  modes  in  these  films. 

To  evaluate  the  electro-optical  properties  we  will  initially  check 
the  performance  of  the  films  in  a  waveguide  deflection-modulation  scheme 
employing  an  interdigital  electrode  structure  deposited  on  the  film 
surface  (Fig. 15).  We  are  in  the  process  of  fabricating  a  suitable  mask 
for  the  electrode  configuration.  Electrode  fingers  will  be  3  mm  long, 

25  pm  wide,  spaced  50  pm  center  to  center.  A  mask  is  also  being  made  for 
a  phase  modulator  which  employs  a  simple  strip  configuration  (Fig.  16) 
consisting  of  two  parallel  electrodes  3  mm  long  separated  by  3  mils.  The 
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film  thickness;  accordingly,  the  3-mil  separation  was  chosen  to  match  the 
75  pm  thickness  of  our  present  films. 


modulation. 
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V.  SPECIAL  INSTRUMENTS  AND  FACILITIES 
5.1  Refractive  Index  Instrument 

In  order  to  make  precise  measurements  of  refractive  index  we  have  built 

13 

a  spectrometer  patterned  after  a  design  by  Rond.  The  instrument  employs 
the  method  of  minimum  deviation.  The  crystal,  ci.*t  xu  the  form  of  a  prism, 
is  mounted  on  a  precision  goniometer  table  and  placed  in  a  collimated  beam 
of  monochromatic  light  so  that  the  orism  deflects,  the  beam.  This  prism  is 
then  rotated  about  an  axis  parallel  to  its  edge  until  the  angular  deflection 
D  is  a  minimum.  The  refractive  indev  is  calculated  from  the  relation 

n  *  sin  I  (A  +  D)/2J/s.in  (A/2) 


where  A  is  the  prism  angle. 

The  instrument  is  now  operating  on  a  routine  basis.  Typical  results 
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An  =  ±  0.002. 


5.2  Clean  Air  Chemistry  Laboratory 

As  noted  in  sections  2.2  and  2.3,  room  temperature  cw  operation  of 
our  sensitized  rare  earth  rainilasers  will  require  low  activator  concentra¬ 
tions  and  will  place  extreme  constraints  on  the  purity  of  the  host  and 
sensitizer  combination.  Our  previous  efforts  to  purify  feeds  have  resulted 
in  the  development  of  a  variety  of  methods  for  purification,  but  most  of 
them  have  been  of  limited  effectiveness  due  to  recontamination  from  the 
laboratory  atmosphere. 

During  this  reporting  period,  we  have  completed  the  installation  of 
a  class  100  clean  air  chemistry  laboratory  designed  and  constructed  by 
Environmental  Air  Control,  Inc.  of  Hagerstown,  Md .  This  laboratory  is  a 
room  within  a  room,  with  an  area  of  approximately  250  sq.  ft.  As  shown  ir 
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Refractive  index  i»f  Yb-Ti-O,  as  a  function  of  wavelength 


Fig.  18  Refractive  index  >f  KTaO„  as  a  function  of  wavelength 


Fig,  20,  there  are  two  front  air  curtain  type  clean  air  hoods,  one  standard 
laboratory  hood,  and  tv?o  eight-foot  laminar. flow  benches,  all  equipped  with 
the  usual  chemistry  bench  services.  This  laboratory  will  be  used 
exclusively  for  the  preparation  of  the  very  high  purity  fluoride  and  oxide 
feed  materials  needed  for  this  program. 


chemistry  room 
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APPENDIX 

Modes  in  a  Two  Dinenslonal  Dielectric  Resonator 

Our  objective  is  to  find  the  normal  modes  that  exist  in  the  two  dimen¬ 
sional  dielectric  cavity  shown  in  Fig.  A-l.  The  cavity  has  a  rectangular 
cross-section  with  a  length  a.  in  the  x-dimension,  b.  in  the  z-dimension.  The 
cavity  index,  n^,  is  less  than  n the  index  of  the  surrounding  medium.  An 
"exact"  solution  to  this  electromagnetic  problem  can  only  be  carried  out 
numerically;  however,  a  reasonably  accurate  closed  form  solution  is  attain¬ 
able  if  we  introduce  some  simplifications  based  on  the  physics  of  the  pro¬ 
blem.  It  id  evident  that  for  a  cavity  with  dimensions  large  compared  to  a 
wavelength  most  of  the  energy  will  be  confined  to  region  1  when  the  condi¬ 
tions  for  total  internal  reflection  are  satisfied.  The  energy  stored  outside 
tb°  nsnity  -Iq  email  hscsiKP  In  this  reeion  the  fields  are  evanescent.  Within 
the  cavity  we  expect  the  fields  to  be  essentially  harmonic.  If  we,  in  fact, 
take  the  solutions  to  be  harmonic,  we  can  match  the  fields  on  the  four  sur¬ 
faces  of  the  cavity  to  exponentially  decaying  fields  outside.  While  this 
solution  will  match  the  boundary  conditions  on  the  faces,  it  will  fail  at 
the  corners.  Despite  this  failure  the  solutions  should  exhibit  little  error 
as  long  as  most  of  the  energy  is  confined  to  the  cavity. 

Proceeding  along  these  lines  we  shall  construct  solutions  which  behave 
in  the  following  way:  in  region  1  the  field  varies  sinusoidally  in  x  and  z; 
in  region  2  and  4  the  variation  is  sinusoidal  in  z  and  exponential  in  x;  in 
3  and  5  we  have  sinusoidal  variation  along  x  and  exponential  dependence  on  z. 
These  assumed  fields  fail  as  an  exact  solution  because  they  are  incompatible 
in  the  "corner"  regions  indicated  by  the  shaded  boundaries  in  Fig;  A-l. 

Our  representation  shall  be  in  ter.r.s  of  modes  TE  and  T?l  to  z.  For  modes 


Fig.  A-l  Two-dimensional  dielectric  cavity 
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The  ratio  of  the  first  equation  to  the  second  gives 
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Equations  A-13  and  A-14  determine  the  k  '  s  and  k  's  and  resonant  frequencies 

X  z 

of  the  even  TM  modes.  These  equations  can  be  transformed  to  a  more  useful 

form  by  using  equations  A-4  through  A-6  to  express  u  and  v  in  terms  of  k  ,  k  , 

X  z 


cr.ci  We  '■'Stair. 


-  kK  cot2  <¥>  +  "  kx+  kl 


-  (~)  >2  cot2  <rf-)  +  (|i)  k2  -  k2  +  k2 
Z  2  E2  X  X  Z 


(A-15) 

(A-16) 


Similarly,  the  characteristic  equations  for  TM  modes  of  symmetry  types 
(2),  (3)  and  (4)  can  be  found.  There  is  complete  duality  between  the  TM  and 
TE  modes  of  the  resonator;  so  the  characteristic  equations  must  be  dual.  The 
results  are  presented  in  Table  A-I.  Although  the  characteristic  equations 
cannot  be  solved  exactly,  simple  graphical  solutions  do  exist.  Consider  the 
even  TM  modes.  From  Eqs.  A-15,  A-16 
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because  in  Eqs.  A-15  and  A-16  the  first  term  on  the  left  in  both  equations  is 
a  negative  quantity. 

Thus  ve  obtain 
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Therefore  there  are  no  solutions  (i.e.,  no  resonant  modes)  for  — <  /2. 

_  n2 
nl  /— 

Only  for  — >  v2  do  resonant  modes  exist  and  the  range  of  k  and  k  is  limited 
^2  x  z 

by  Eqs.  A-19and  A-20  to  the  region  of  k-space  shown  in  Fig.  A-2  where 

.  -1  r,n^2  V2  TT 

‘•“i‘  lv:,2'  -  i  ~  2' 

Since  cot2  )  goes  from  +°°  to  0  to  +00  for  <_  k  <  where 

Z  a  X  a 

p  is  any  positive  integer,  there  are  at  most  two  solutions  of  Eq.  A-15  in  the 

interval  -^-^<k  <  -2— for  k  fixed.  Furthermore,  only  one  of  these 
a  —  x  —  a  z 

solutions  satisfies  Eq.  A-13,  since  u  is  positive.  The  same  considerations 
are  true  for  Eq.  A-6.  Therefore,  there  is  only  one  resonant  mode  in  the 
interval 


am  <  k  <  aifetu  and  am  <  k  <  aiSs+n 

a  x  a  b  2  b 


(A-21) 


where  p  and  q  are  any  positive  integers  and  k  and  k  also  satisfy  Eqs.  A-19 

x  z 

and  A-20.  So  the  mode  density  is  -^5-  in  the  allowed  region.  Finally,  the 

4tt 

number  of  resonant  modes  in  a  frequency  interval  Aw  is 

,  ,  r  n,  2  -1/2  ,  n2w 

All  =  4{-l  +  -  tan  1  r  (— )  -  ll  }  ™  ~  Aw 

L  n2  '  A;'  c2  (A-  -22) 
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Mg.  A-2  Totally  reflected  nodes  occur  within  the 

region  of  k-space  subtended  by  the  angle  60 
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where  the  factor  4  accounts  for  the  four  types  of  TM  mode  symmetry.  The  result 
Is  the  same  for  TE  modes. 

The  preceding  analytical  treatment  of  modes  ir.  a  two  dimensional  rec¬ 
tangular  dielectric  cavity  is  succinctly  summarized  by  the  geometric  construc¬ 
tion  shown  in  Fig.  A-3.  In  k-space  draw  the  two  circles  k,  *=  n,k  ,  k,.  =  n0k 

1  1  o  l  l  o 

(k2<ki').  Construct  the  horizontal  and  vertical  tangents  to  the  inner  circle  k^. 
In  Figure  A-3  these  are  1-1'  and  2-2'.  Find  the  points  A,  B  where  the  tangents 
intersect  the  outer  circle  k^.  The  radial  lines  drawn  from  the  origin  to  A 
and  B  subtend  the  sector  in  k-sprc^  in  which  the  totally  re  lecting  modes  are 
contained  (the  shaded  region  in  the  figure).  It  should  be  evident  from  the 
construction  that  the  angular  width  of  the  sector  decreases  as  k£  increases 
relative  to  k^  and  the  sector  collapses  to  zero  width  when  the  ratio  k^/k^  is 
such  that  the  points  A  «m<2  L>  coincide.  The  critical  ratio,  evident  trom  the. 
geometry,  is  ^/k^  =  1//2.  For  larger  ratios  the  cavity  will  support  no  in¬ 
ternally  reflecting  modes. 

In  Fig.  A-2  the  modes  the  shown  equally  spaced  in  intervals  tr/a  in  k^ 

and  ir/b  in  k  .  This  representation  while  not  exact  is  an  excellent  approxi- 
z 

mation  to  the  mode  spacing  when  the  cavity  dimensions  are  very  large  compared 
to  the  free  space  wavelength  Aq.  In  this  domain,  which  applies  to  our  experi¬ 
mental  situation,  we  can  write 

.2  2,  2  .  pTT.  2  ,  qTT\2 

k  -  n  k  “  (  “0  +  (  *-) 

1  1  o  a  a 


t-  *  ^ 


(A-23) 
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where  k  ■  2ir/X  is  the  free  space  wave  vector  corresponding  to  the  wavelength 
o  o 

X  and  p  and  q  are  positive  integers.  An  expression  for  the  incremental  mode 

o  5  - 

spacing  at  large  wave  numbers  can  be  obtained  if  we  take  the  derivative  of 


Eq.  A-23: 


"lVko  "  t  <?><>»  +  <f>  <?>*>» 


(A-2A) 


For  U£  and  chosen  so  that  the  sector  of  allowed  inodes  in  k-space  is  small 
we  have 

(A-25, 

Upon  substituting  this  result  into  Eq.  A-2A  we  obtain,  after  some  simplifica¬ 
tion,  . 


AX 


-  ^  ) 


<x 


bv  2s) 


where  Ap  and  Aq  are  positive  ani  negative  integers  ar.d  may  include  zero. 

We  note  in  conclusion  that  each  point  in  Fig.  A-3  represents  eight  nodes: 
one  set  of  A  symmetry  types  for  TE  modes  plus  a  similar  sec  for  TM.  These 
8  modes  are  not,  in  general  all  degenerate,  but  in  the  limit  a,b»X  they  tend 
to  coalesce  to  a  single  point  in  k-space. 
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.Summary 

This  project  has  two  objectives:  to  unravel  the  complex  problem 
of  superconductivity  of  high-Tc  A15  materials  by  tunneling  experiments; 

i 

and  to  develop  materials  for  practical  superconducting  microwave 
cavities.  In  this  report  we  describe  our  progress  towards  the  first 
objective.  The  importance  of  understanding  the  A15  group  stems  from 
the  fact  that  they  have  the  highest  known  transition  temperatures  and 
other  technologically  useful  properties;  and,  such  knowledge  should 
show  the  way  to  still  higher  transition  temperatures  and  greater 
technical  utility.  Our  approach  to  the  problem  is  to  perform  high- 
quality  tunneling  experiments  A15's  with  the  highest  Tc's  and  to 
develop  computer  software  which  is  capable  of  deconvoluting  the  data 
into  phonon  spectra. 

Progress  in  the  past  six  months  has  centered  on  two  areas:  (a) 
tunneling  into  Nb^CAl  gGe  2)  and  (b)  deconvolution  of  phonon  spectra 
from  S- I -S'  tunnel  junctions. 

(a)  We  have  achieved,  in  our  initial  exploratory  attempts, 
superconductive  tunneling  into  Nbj(Al,Ge)  ribbon  substrates.  Although 
ribbons  of  the  U5  material  offer  good  flat  surfaces  for  barrier 
formation,  and  thus  a  high  probability  fox  successful  tunneling,  the 
metallurgical  condition  of  these  particular  ribbons  is  not  sufficiently 
bulk-like  that  our  results  can  be  considered  representative.  On  the 
other  hand,  our  work  with  the  ribbons  indicates  that  the  oxides  of 
Nb^CAljGe)  are  very  fragile,  deteriorating  under  sustained  bias  or 
temperature.  The  most  durable  barriers  were  formed  after  glow  discharge 
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cleaning  of  the  surface  and  thermal  oxidation  of  the  resulting  substrate. 
T*  e  best  quality  tunneling  data  resulted  when  the  ribbon,  before 
barrier  formation,  was  not  subject  to  the  usual  low  temperature 
ordering  anneal  (725  C) . 

(b)  Our  modification  of  the  McMillan- Rowe 11  Inversion  Routine 
has  resulted  in  a  more  (numerically)  accurate  set  of  computer  programs 
which  can  be  used  to  deconvolute  the  phonon  spectra  from  S-I-S' 
tunnel  junctions  where  both  S  and  S'  can  simultaneously  be  non-BCS 
superconductors.  The  only  constraint  on  the  deconvolution  program  is 
that  the  phonon  spectrum  of  one  electrode  must  be  known;  i.e.,  both 
electrodes  can  not  be  total  unknowns. 

One  implication  of  the  above  is  that  we  can  look,  with  greater 
precision  than  we  thought  possible,  at  more  kinds  of  junctions  than 
we  thought  possible.  The  init.ial  successes  at  tunnel  junction 
fabrication  on  ribbon  substrates  indicate  that  the  goal  of  this 
research  is  definitely  feasible,  and  we  are  now  tackling  the  substrate 
surface  problem,  with  the  same  optimism.  We  should  remark  at  this 
point  that  recent  developments  in  Nb^Ge  in  other  laboratories 
indicate  that  the  potential  for  still  higher  T£  A15  compounds  is 
very  far  from  exhausted;  thus  the  original  premise  of  this  research, 
the  importance  of  fundamental  understanding  of  A15  superconductivity, 
takes  on  still  greater  importance. 
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I.  Introduction 

In  the  last  decade,  electron  tunneling  has  come  to  be  recognized 
as  the  most  important  technique  for  studying  the  properties  of  metals 
in  both  the  superconducting  and  normal  metal  states.  In  a  tunnel 
junction,  if  one  (or  both)  of  the  metals  is  in  the  superconducting 
state  at  least  three  important  properties  of  the  metal(s)  can  be 
determined  from  I-V  and  dl/dV  -  V  characteristics:  the  electronic 
density  of  states  of  the  metal  in  its  superconducting  state,  the  phonon 
density  of  states,  and  the  electron-phonon  coupling  constant. 

The  purpose  of  the  present  research  is  to  determine  these  three 
parameters  for  representative  members  of  the  A15  class  of  compounds . 
These  materials  are  especially  interesting  for  two  reasons.  First,  the 
three  highest  known  superconducting  transition  temperature  (Tc) 
compounds  ar£A15's:  Nb3Ge,  NbjGa,  Nb3(Al  75Ge  2$) .  (In  addition, 
most  materials  with  T  ■  >  15K  are  A15  compounds.)  Secondly,  many  of  the 
high  Tc  A15's  are  structurally  unstable.  It  is  important,  for  future 
research,  as  well  as  for  eventual  applications  of  superconductivity, 
to  know  if  the  high  Tc's  attained  in  these  materials  are  in  any  way 
connected  with  their  instabilities.  There  is  even  a  question 
as  to  whether  the  "usual"  BCS  pairing  mechanism  adequately  describes 
the  superconducting  properties  of  these  alloys. 

In  this  report  we  discuss  two  different  aspects  of  our  research 
to  date.  The  first  section  contains  a  description  and  analysis  of  the 
tunne'.ing  data  we  have  obtained  from  Nb3(Al,Ge)  ribbon-substrate  tunnel 
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junctions  with  both  oxide  and  carbon  barriers.*  We  have  found  that 
whereas  it  is  possible  to  fabricate  such  junctions  using  a  lead- 

( 

bismuth  counterelectrode  and  to  obtain  superconductive  tunneling,  they 
are  relatively  fragile  and  deteriorate  rapidly  under  temperature 
recycling  and/or  sustained  app)!ed  bias.  The  methods  of  barrier 
fabrication  are  extensively  reviewed  here.  The  same  material  used 
to  make  the  ribbon  substrates  was  also  rf  sputtered  into  a  supercon¬ 
ducting  film  by  the  IBM,  Thomas  J.  Watson  Laboratory.  A  film  junction 
with  a  carbon  barrier  and  lead-bismuth  counterelectrode  was  then 
fabricated  ar.d  tested.  The  results  from  this  junction  are  surprisingly 
similar  to  carbon  junctions  made  on  the  ribbon  substrate.  All  the 
carbon  barrier  samples  have  a  characteristic  structure  in  their  deriva¬ 
tives  which  we  have  not  yet  identified. 

The  second  section  is  a  reproduction  of  the  thesis  of  Alan  Jay 
Dubin,  S.B.,  M.I.T. ,  August  1973.  In  this  thesis  the  concepts  behind 
superconducting  tunneling  aie  explored  and  a  description  of  the  method 
of  unravejxng  actual  tunneling  data  to  reveal  the  superconducting 
electron  spectrum,  the  phonon  spectrum,  and  the  electron-phonon  inter¬ 
action  constant  is  given.  The  techniques  described  are  a  modification 
of  the  original  McMillan- Rowe  11  phonon  inversion  routines.  Major 
changes  have,  been  made  in  two  of  the  three  programs.  The  first  modifi¬ 
cation  is  a  program  introduction  comprised  of  all  of  the  software 
necessary  to  convert  the  output  of  conventional  tunneling  electronics 


*  A  brief  communication  on  this  work  will  appear  in  Physics  Letters 
(in  press) . 
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into  absolute  data  characteristics.  The  second  stage  of  this  routine, 
which  attempts  to  isolate  the  characteristics  of  the  unknown  (or 
unexplored)  counter-electrode,  has  been  extensively  revised.  The 
numerical  accuracy  o.t  this  stage  has  been  improved  by  more  than  an 
order  of  magnitude  and  the  new  program  will  enable  double  superconductor 
junctions  in  which  neither  electrode  is  a  BCS  (or  weak-coupling)  material 
to  be  analyzed.  The  latter  is  absolutely  necessary  for  tunneling 
investigations  on  transition  metal  elements  or  alloys  since  these  are 
not  weak-coupling  materials,  and  because  it  is  found  experimentally 
that  they  form  their  best  tunnel  junctions  with  other  non-BCS  materials. 
The  net  result  of  these  modifications  is,  we  feel,  an  inversion  routine 
which  has  greater  accuracy  and  flexibility  than  the  original  inversion 
routine,  is  snore  specific  m  the  raw  data  handling  techniques,  and, 
finally,  has  a  far  wider  range  of  application  than  the  original 
program. 

II.  Exploratory  Tunneling  Results  on  Nb-Al-Ge 

We  have  succeeded  in  producing  Nb_(Al  0Ge  0)  -  barrier  -  Pb  _Bi  , 
tunneling  junctions  on  ribbon  substrates  of  the  A15  compound.  In 
this  discussion  we  shall  discuss  steps  taken  to  date  to  obtain  tunneling 
in  Nbj(Al,Ge)  -  Pb  ?Bi  3  junctions  and  attempt  to  relate  these  steps 
to  the  observed  junction  behavior. 

(a)  Ribbon  Fabrication 

The  substrates  were  ribbons  of  Nb^(Al,Ge)  fabricated  by  a  process 
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into  absolute  data  characteristics.  The  second  stage  of  this  routine, 
which  attempts  to  isolate  the  characteristics  of  the  unknown  (or 
unexplored)  counter-electrode,  has  been  extensively  revised.  The 
numerical  accuracy  of  this  stage  has  been  improved  by  more  than  an 
order  of  magnitude  and  the  new  program  will  enable  double  superconductor 
junctions  in  which  neither  electrode  is  a  BCS  (or  weak-coupling)  material 
to  be  analyzed.  The  latter  is  absolutely  necessary  for  tunneling 
investigations  on  transition  metal  elements  or  alloys  since  these  are 
not  weak-coupling  materials,  and  because  it  is  found  experimentally 
that  they  form  their  best  tunnel  junctions  with  other  non-BCS  materials. 
The  net  result  of  these  modifications  is,  we  feel,  an  inversion  routine 
which  has  greater  accuracy  and  flexibility  than  the  original  inversion 
routine,  is  more  specific  in  the  raw  data  handling  techniques,  and, 
finally,  has  a  far  wider  range  of  application  than  the  original 
program. 

II.  Exploratory  Tunneling  Results  on  Nb-Ai-Ge 

We  have  succeeded  in  producing  Nb,(Al  QGe  ..)  -  barrier  -  Pb  _Bi 
tunneling  junctions  on  ribbon  substrates  of  the  A15  compound.  In 
this  discussion  we  shall  discuss  steps  taken  to  date  to  obtain  tunneling 
in  Nbj(Al,Ge)  -  Pb  ^Bi  ^  junctions  and  attempt  to  relate  these  steps 
to  the  observed  junction  behavior. 

(a)  Ribbon  Fabrication 

The  substrates  were  ribbons  of  Nbj(Al,Ge)  fabricate!  by  a  process 
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developed  in  our  laboratory  .  Hie  ribbon  geometry  is  especially 
suited  for  tunneling  work  (as  films  are)  because  of  the  existence  and 
accessibility  cf  flat  surfaces.  The  ribbons  are  made  by  first  enclosing 
powders  of  Nb,  Al,  and  Ge  in  the  appropriate  composition  in  a  cupro¬ 
nickel  tube;  then  reducing  the  tube  in  size  by  a  series  of  swaging  and 
rolling  steps;  and  then  finally  removing  the  cupronickel  sheath  with 
nitric  acid.  (The  resulting  ribbons  in  our  laboratory  were  0.1  mm  x 
2.5  nun  in  cross  section.)  The  ribbon  is  stored  in  this  condition 
until  a  series  of  junctions  is  to  be  made.  It  is  then  heated  resistively 
to  1700°C  in  an  Argon  atmosphere  to  produce  the  A15  compound.  X-ray 
diffraction  and  metallcgraphic  studies  on  several  of  these  samples 
showed  that  90  -  95%  of  a  given  ribbon  possessed  the  A15  structure.  The 
typical  Tc  of  such  ribbons  was  17.5  K  2, 

Prior  to  junction  fabrication  at  least  half  of  the  samples  in 
the  present  investigation  were  then  annealed  at  725  C  for  90  hours 
in  less  than  10  Torr.  This  is  done  in  order  tc  increase  atomic 
order.  After  this  anneal,  there  was  a  general  increase  of  about  1.5  K 
in  a  sample's  3. 

The  last  step  before  the  barrier  fabrication  was  to  etch  the 
surface  of  the  ribbons  to  remove  the  outside  layers  of  the  ribbon  and, 
consequently,  any  depleted  layers  at  the  surface(and  also  any  surface 
contaminants  such  as  nickel  atoms  remaining  after  the  bath  in  nitric  acid) . 

Both  chemically  etched  samples^  and  those  cleaned  by  glow  discharge 
appeared  smooth  down  to  1  micron.  Figures  1  and  2  are  the  SEM 

*  References  and  pictures  for  this  section  of  the  report  are  found  on  pp  16  ff. 
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photographs  of  the  ribbon  surfaces  after  each  of  the  two  etching 
procedures.  Both  show  a  considerable  porosity  which  is  a  remnant 
of  the  sintering  process. 

In  the  glow  discharge  cleaning  prior  to  barrier  formation,  the 
ribbon  was  maintained  at  a  negative  potential  difference  greater  than 
1000  V  in  3  x  10  2  Torr  of  Argon  so  that  oxidized  and  cold-worked 
material  was  sputtered  off.  (See  Fig.  3).  This  cleaning  (or  that 
with  the  potential  reversed)  proved  to  be  essential  to  obtaining 
junctions  which  exhibited  tunneling  structure. 

(b)  Barrier  Formation 

As  is  well  known,  the  quality  of  the  barrier  is  crucial  for 
obtaining  high  resolution  tunneling  data.  Excess  currents  and 
barrier-related  structure  in  tunneling  characteristics  can  be  so 
severe  as  to  totally  mask  superconducting  characteristics.  An 
additional  problem  arises  for  superconductors  of  extremely  short 
coherence  lengths,  such  as  A15  alloys:  the  interface  region  between 
barrier  and  bulk  material  constitutes  a  significant  fraction  of  the 
region  sampled  by  the  tunneling  electrons,  even  for  interfaces 
which  would  ordinarily  be  considered  quite  sharp. 

In  the  present  work  two  different  barrier  types  were  tested: 
oxide  and  carbon.  The  oxide  harrier  was  formed  either  by  anodization 
in  the  glow  discharge  system  (Fig.  3)  without  exposure  to  atmosphere, 
or  by  placing  the  samples  in  a  free  flow  of  oxygen  at  150  C  (thermal 
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oxide).  In  the  former  procedure  the  sample  was  held  a  few  hundred 

_2 

volts  above  the  cathode  in  5  x  10  Torr  of  oxygen.  Although  it 

f 

was  not  difficult  to  oxidize  the  samples,  and  they  were  of  an  insulating 
nature,  we  later  found  that  they  are  apparently  both  physically 
and  electrically  fragile  when  made  in  the  above  ways.  Junctions  were 
made  with  carbon  film  barrier  on  4  ribbons  which  had  previously  had 
a  725  C  anneal.  The  carbon  barriers  were  fabricated  in  two  different 

O 

ways.  In  one  instance  paracrystallini  C  films  of  about  150  A 
thickness  were  arc-evaporated  onto  the  cleaned  ribbon  as  barriers^. 

7 

In  the  other,  an  electron  beam  deposition  technique  was  used  to 

O 

fabricate  barriers  100  +_  30  A.  Although  the  resistances  and  physical 
characteristics  of  these  two  types  of  carbon  barriers  are  significantly 
different,  their  tunneling  characteristics  are  almost  indistinguishable. 
The  lack  of  durability  found  in  the  oxide  barriers  extends  to  the  carbon 
barriers  as  well. 

After  barrier  fabrication  by  whatever  technique,  junction 

O 

boundaries  were  then  out- lined  by  a  formvar  mask,  and  a  300.  A  thick 
Pb  ^Bi  j  counterelectrode  deposited  to  complete  the  tunneling  sandwich. 
Conventional  four  point  probe  connections  were  then  made,  and  the 
entire  junction  covered  with  a  protective  layer  of  formvar.  Two  or 
three  junctions  -/ere  usually  made  cn  any  one  ribbon  sample, 

(c)  Nbj(Al,Ge)  Film  Junctions 

A  number  of  unreacted  ribbon  samples  wi.re  given  to  Mr.  Walter 
Kateley  at  the  IBM,  Thomas  J.  Watson  Laboratory,  from  which  to  rf 
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jputter  NbjCAl.Ge)  films  .  We  sought  to  compare  the  quality  of 

tunneling  of  the  ribbon,  to  the  film  substrates  The  one  film  on  which 

we  have  completed  tests  had  ^  average  thickness  of  3000  A,  was 

metallic  in  color  (not  quite  opaque) ,  and  had  a  room  temperature 

-2 

resisti'  ity  ^3  x  10  ft-cm.  The  transition  temperature  of  the 
unannealed  film  was  -  1.6  K  with  a  transition  interval  M).5  K. 

Four  carbon  barrier  junctions  (electron  beam  deposition)  were 
made  on  this  film,  only  two  of  which  were  operable. 

(d)  Results 

Our  best  quality  results,  to  date  have  been  obtained  with  an 
up annealed  ribbon  which  had  been  mounted  as  the  anode  in  the  Ar  glow 
discharge  cleaning,  and  which  had  a  thermal  oxide  barrier.  The 

Q 

tunneling  characteristics  are  shown  in  Figs,  d  and  5  .  At  4.2  K 

this  junction  had  a  resistance  of  270  ft.  It  broke  down  shortly  after 
the  second  sum  peak  bias  was  reached. 

The  value  of  2A  of  Nbj(Al,Ge)  obtained  from  this  junction  was 
0.70  +_  .05  meV  at  4.2  K  and  1.21  meV  +_  .05  at  2.6  K.  If  the  energy 
gap  followed  the  usual  empirical  temperature  dependence1*1,  wj  would 
obtain  the  values:  T£  =  4.6  K,  2A(0)  =  1.3  meV  and  (2A/kTJ  =  3.3. 

A  second  junction  on  the  same  ribbon  yielded  a  value:  2A  =  0.65  +0.1  meV 
at  4.2  K.  It  deteriorated  even  more  quickly,  however,  disintegrating 
before  it  was  possible  to  pump  from  4.2  to  2.5  K.  Strangely  enough 
it  had  a  junction  resistance  of  1200ft  at  4.2  K.  Several  other  ribbons 


Co< 


cleaned  via  glow  discharge  and  having  a  thermal  o^ide  barrier, 

were  obtained  with  structure  above  4  meV,  but  the  attendent  high  excess 

currents  made  an  accurate  analysis  questionable. 

As  discussed  previously,  for  tunneling  data  to  reflect  the  bulk 

properties  of  the  sample  the  depth  of  the  surface  region  must  be 

less  than  the  coherence  length  for  superconductivity.  Since  in  Al5's 

0  11 

the  coherence  length  is  20-80  A  ,  we  believe  this  study,  as  well 

12 

as  all  previous  ones  ,  gave  surface  limited  results.  Our  results 
are  probably  influenced  by  inadequate  cleaning  of  the  surface, 
the  presence  of  segregated  A1  or  other  inhomogeneities  in  the  ribbon 
surface,  and/or  contamination  from  the  cupronickel  tubing  used  in 
fabricating  the  ribbons. 

(e)  Discussion 

In  all  we  have  investigated  approximately  30  ribbons  and  70 

junctions.  All  but  five  of  these  junctions  had  oxide  barriers.  For 

the  oxide  junctions  only  4  gave  rise  to  analyzable  superconductive 
14 

tunneling  .  "ihree  oxide  junctions  had  recognizable  characteristics 
and  all  three  had  thermal  oxide  barriers;  none  were  subjected  to  the 
low  temperature  (725  C)  anneal.  The  other  junctions  made  on  unannealed 
ribbon  were  generally  of  poor  quality  but  had  sharp  structure  in  the 
low-bias  portion  of  the  derivative  curve  (see  Fig.  6).  This  might  be 
explained  either  by  A1  segregation  in  the  substrate  or  by  drastic 
reduction  of  T£  and  A  for  material  near  the  ribbon  surface.  Subsequent 
evaluation  of  X-ray  data  on  these  ribbons  indicated  a  rather  high 
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percentage  of  unreacted  Germanium.  This  would  lead  to  the  formation 
of  Al-rich  non-equilibrium  phases  in  the  material.  In  any  case,  the 
best  tunneling  data  was  observed  on  the  ribbon  subjected  to  the 
reverse  glow  discharge  cleaning. 

Junctions  fabricated  on  annealed  ribbons  tended  to  be  of  even 

poorer  quality  than  on  the  unannealcd  substrates.  Instead  of 

having  the  low  bias  structure  shown  in  Fig.  6,  they  tended  to  exhibit 

structure  at  biases  M  meV  above  the  gap  edge,  consistent  with 

13 

phonon  structure  observable  in  lead  and  lead  alloys.  Figure  7  is  one 
such  derivative  trace  taken  at  1.2  K.  This  particular  trace  and  others 
with  such  clear-cut  structure  all  corresponded  to  thermal  oxide 
barriers, although  most  of  the  others  had  even  higher  excess  currents 
and  thus  a  high  parabolic  background  in  the  derivative  scan. 

Figure  8  from  an  anodic  barrier  sample  is  c  typical  junction  with 
the  high  excess  current.  The  high  excess  currents  of  most  junctions 
could  be  caused  by  poor  surfaces  (e.g.,  Al-depletion  during  annealing) 
or  an  inadequate  tunneling  barrier.  Thus  far,  we  have  not  been  able 
to  pin  down  tl  -  cause  of  this  tehavior.  Micro-shorts  and  hysteresis 
have  been  present  in  many  junctions  but  no  evidence  of  Josephson 
behavior  has  been  observed.  All  junctions  have  exhibited  breakdown 
after  two  or  three  derivative  scans  at  liquid  helium  temperatures; 
many  have  deteriorated  in  the  time  it  takes  to  pump  the  system  down 
from  4.2  to  1.2  K.  The  cause  of  this  behavior  is  unknown. 
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All  of  the  carbon  barrier  junctions  had  unusual  structure  at 
low  temperatures.  First  derivative  curves  for  both  a  ribbon  substrate 
junction  and  the  sputtered  Nb3(Al,Ge)  film  junction  are  given  in  Fig. 

9.  The  characteristic  "three  bump  structure"  appears  in  all  of  the 
dl/dV  -  V  curves  for  carbon  barrier  junctions  below  4.2  K.  (In  the 
film  data  it  may  appear  at  4.2  K  but  the  derivative  is  not  clean 
enough  to  say  unequivocally  that  it  is  there.)  Just  as  in  these  two 
traces,  however,  the  bias  at  which  these  "bumps"  appear  is  not 
related  in  any  obvious  way  to  the  usual  derivative  analysis  procedures. 

It  does  not  seem  likely  that  this  structure  is  due  to  the  carbon  alone 
'despite  the  fact  that  there  is  no  evidence  of  its  appearonce) 
in  the  oxide  tunnel  junctions)  because  we  have  made  a  number  of 
other  S-I-S'  junctions  with  carbon  barriers  sandwiched  between  In 
and  Al,  Sn,  and  In  counterele.ctrodes  and  see  absolutely  no  evidence 
for  its  existence  in  the  derivative  scans.  We  are  presently  looking 
into  the  possibility  of  interaction  between  the  carbon  and  either  Nb3(Al,Ge) 
or  the  lead -bismuth  counterelectrode.  Nonetheless,  it  would  have  to 
be  a  somewhat  special  type  of  interaction,  if  its  characteristic  energy 
is  that  fickle.  Clearly  if  there  is  any  kind  of  chemical  bonding,  or 
spreading  of  the  barrier  -  material  interface,  then  different  methods 
of  obtaining  artificial  barriers  will  have  to  be  explored.  Unfortunately, 
these  barriers  appear  to  be  as  fragile  as  their  oxide  counterparts. 
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(f)  Summary 

We  have  achieved  superconducting  tunneling  on  ribbon  substrates 
of  Nb3(Al,Ge) ;  of  all  the  samples  exhibiting  non-normal  metal  tunneling 
only  two  had  ideal  junction  characteristics.  Junction  barriers  were 
initially  oxides  prepared  by  one  of  a  variety  of  methods;  e.g., 
anodization  by  glow  discharge  and/or  thermally  grown  in  oxygen.  Prior 
to  oxidation  the  ribhon  surfaces  were  cleaned  by  glow  discharge  in  an 
Argon  atmosphere  and/or  chemically  etched.  Irrespective  of  the 
combination  of  tnese  procedures  and  whether  the  ribbon  was  annealed 
or  not,  junctions  were  invariably  of  very  poor  quality.  Typically, 
when  structure  did  appear  in  the  derivative  characteristic,  it  implied 
unusual  gap  values  for  the  Nbj(Al,Ge)  or  the  Pb-Bi  alloy,  or  both. 

All  characteristics  with  distinctive  tunneling  structure  were  thermally 
oxidized  unannealed  ribbon  substrates.  Although  many  of  these  junctions 

were  not  of  good  tunneling  quality,  the  data  from  them  indicated  a 
prob.lem  of  segregated  A1  in  the  ribbons,  particularly  after  a  low 
temperature  anneal.  Subsequent  X-ray  examination  confirmed  the  presence 
of  unreacted  Ge  in  our  samples. 

In  an  effort  to  obtain  tougher  barriers,  we  investigated  carbon 
as  a  barrier  material.  Although  the  tunneling  data  obtained  from 
such  junctions  was  clearly  not  that  of  two  normal  metals  separated  by  a 
semiconductor,  we  'nave  not  yet  been  able  to  unravel  what  structures 
we  have  observed.  Film  junctions  made  by  rf  sputtering  our  unreacted 
ribbon  material  and  then  depositing  both  a  carbon  barrier  and  a 
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Pb  _Bi  _  counter-electrode  have  characteristics  almost  identical  to 
the  analogous  ribbon- substrate  carbon  junctions.  In  addition,  the 
carbon  junctions  have  proved  to  be  no  more  durable  than  similar  ones 
with  oxide  barriers.  The  initial  success  wiu\  both  the  oxide  and  the 
carbon  barriers,  however,  encourages  us  that  f>’"ther  research  into 
available  fabrication  techniques  will  probably  result  in  more  adequate 
barriers . 

(gj  Direction  of  Future  Work 

Clearly  the  success  we  have  had  with  our  preliminary  investigations 
on  the  ribbon  substrates  dictates  that  in  our  continuing  work,  special 
care  must  be  taken  to  optimize  the  homogeneity  and  smoothness  of  the 
A15  surfaces.  The  major  emphasis  in  this  research  is  on  fabricating 
suitable  barriers  for  the  junctions.  Since  the  film  geometry  also 
provides  smooth  surfaces  and  dur  results  indicate  a  behavior  analogous 
to  that  of  the  ribbons,  we  shall  continue  our  present  collaboration 
with  IBM  in  the  A15  film  investigation. 

The  best-characterized  form  of  Nbj(Al,Ge)  is  the  arc-melted  bulk 
material;  with  this  in  mind,  we  have  begun  to  devise  techniques  for 
exposing  surfaces  of  such  material  which  promise  to  be  characteristic 
of  the  bulk.  Typical  alternatives  under  investigation  are  controlled 
brittle  fracture,  controlled  etching  and  various  metal lographic 
techniques.  Figure  10  is  a  typi  al  SEM  picture  of  the  surface  of  one 
of  our  recently  prepared  samples.  The  presence  of  such  cleavage 
steps  in  quantity  makes  this  particular  surface  unsuitable  for  our 
purposes.  We  hope  to  produce  better  surfaces  sorn. 
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Figure  Captions 

Figure  1  SEM  photograph  of  chemically  etJied  Nb^fAl.Ge)  ribbon  surface. 

Figure  2  SEM  photograph  of  the  glow  discharge  etched  ribbon  surface. 

Figure  3  Glow  discharge  configuration. 

Figure  4  I-V  characteristic  of  a  Nb_(Al,Ge)  -  Ox  -  Pb  7Bi  ,  junction 
of  270.0.  at  4.2  K.  ./  .3 

Figure  5  dl/dV  -  V  characteristic  of  the  same  junction  as  in  Fig.  4. 

Figure  6  Derivative  scan  of  an  unannealed  ribbon  junction  with  a  thermal 
oxide  barrier  (4.2  K  resistance  =  12A)  at  1.2  K. 

Figure  7  Derivative  scan  of  an  annealed  ribbon  junction  with  a 

thermal  oxide  barrier  (4.2  K  resistance  =  10TT)  at  1.2  K. 

Figure  8  Derivative  scan  of  an  annealed  ribbon  junction  with  an  anodic 
oxide  barrier  (4.2  K  resistance  =  4/L)  at  1.2  K. 

Figure  9  A  first  derivative  scan  of  two  different  carbon  barrier 

junctions.  The  upper  curve  (A)  is  a  ribbon  substrate  sample 
at  1.2  K.  The  lower  curve  (B)  is  an  rf  sputtered  film  junction 
at  1.2  K. 

Figure  10  SEM  photograph  of  a  fragment  of  an  arc-cast  Nb3(Al,Ge)  ingot. 
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A  MODIFICATION  OF  THE  McMILLAN- ROV7ELL 


TUNNELING  DATA  INVERSION  PROGRAM 


ABSTRACT 

The  McMil1 in-Rowell  Inversion  Program  for 
extracting  phonon  spectra  and  other  normal  state 
properties  from  superconducting  tunnel  junction 
data  has  been  modified  for  more  efficiency  and 
greater  ease  of  data  handling.  In  addition ,  this 
modified  program  will  enable  double  superconductor 
junctions  in  which  both  electrodes  are  strong¬ 
coupling  superconductors  to  be  analyzed.  The 
sensitivity  cf  the  program  to  deviations  from  weak- 
coupling  superconductivity  has  been  increased  so 
that  previously  unanalyzed  materials  such  as  Niobium 
can  now  be  investigated.  Having  obtained  the  phonon 
spectrum  of  a  given  material,  ether  transport  prop¬ 
erties  which  depend  upon  the  electron-phonon  inter¬ 
action  can  then  be  calculated;  various  examples  of 
such  properties  are  given. 
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CHAPTER  I 


A.  Introduction 

In  the  last  decade,  electron  tunneling  has  be¬ 
come  recognized  as  one  of  the  most  important  tech¬ 
niques  for  studying  the  properties  of  conductors  in 
both  the  superconducting  and  normal  metal  states. 

A  tunr  -1  junction  is  a  metal-insulator-metal  sand¬ 
wich  through  which  current  flows  via  a  quantum  me- 
canical  tunneling  process.  Such  devices  are  partic¬ 
ularly  important  in  the  study  of  superconducting 
materials*-  If  one  (or  both)  of  the  metals  is  in 
the  superconducting  state  at  least  three  important 
properties  of  the  metal  (metals)  can  be  determined 
from  I-V  and  dl/dV  -  V  characteristics  of  the  junc¬ 
tion.:  the  electronic  density  of  states  of  the  metal 
in  the  superconducting  state,  the  phonon  density  of 
states,  and  the  electron-phonon  coupling  constant. 

The  BCS  theory  of  superconductivity  (named  after 
the  formulators :  Bardeen,  Cooper,  and  Schrieffer)  (3) 
predicts  the  existence  of  an  energy  gap  at  i_he  Fermi 
surface  of  a  superconducting  metal.  In  a  normal 
metal,  at  T=0,  the  Fermi  surface  in  energy  space  is 
sharply  defined:  all  states  below  the  Fermi  energy. 


Ep,  are  occupied  while  above  they  are  empty.  In  a 
superconductor  there  exists  an  energy  gap  of  2A  in  the 
electronic  excitation  spectrum  centered  about  the 
Fermi  energy.  This  gap  is  caused  by  the  binding  of 
electrons  to  form  Cooper  pairs ;  a  phonon  of  energy  2A 
is  required  to  break  these  pairs,  that  is,  to  excite 
single  electron-like  particles. 

The  purpose  of  this  thesis  is  to  describe  how 
tunneling  data  can  be  processed  to  obtain  the  elec¬ 
tronic  spectrum(s),  the  phonon  spectrum (s ) ,  and  the 
electron-phonon  interaction  constant (s)  appropriate 
to  the  metals  comprising  a  superconductor  junction. 

The  initial  "inversion  process"  was  originated  by 
McMillan  and  implemented  by  McMillan  and  Rowell  on 
lead  tunnel  junctions  (15) .  It  has  since  been  used 
by  many  others  on  various  materials  (1,5,12,19).  Its 
basic  limitation  is  that  it  assumes  a  BCS-model  of 
electron  pairing  due  to  electron-phonon  interactions 
in  the  metal.  Although  we  have  used  the  basic  outline 
for  inversion  initiated  by  McMillan,  substantial 
accuracy  has  been  added  to  the  program  by  more  effi¬ 
cient  and  flexible  data  handling. 


In  this  chapter  the  tunneling  process  in  both 
normal  and  superconductor  junctions  is  described. 

The  emphasis  is  on  how  the  superconducting  electronic 
density  of  states  is  reflected  in  the  conductance 
data  of  tunnel  junctions.  In  Chapter  II  the  strong¬ 
coupling  modification  of  the  BCS  theory  is  briefly 
discussed  and  various  parameters  needed  for  the  sub¬ 
sequent  calculation  of  the  phonon  spectrum  and  elec¬ 
tron-phonon  coupling  constant  are  introduced.  Chap¬ 
ter  III  contains  a  discussion  of  the  various  stages 
in  the  computer  program  which  enable  the  parameter 
calculations.  It  is  hoped  that  the  reader  will  ob¬ 
tain  an  overview  of  the  total  inversion  process  in 
this  chapter.  The  appendix  contains  a  computer  print¬ 
out  of  the  inversion  program  along  with  introductory 
notes  and  comments  intended  to  help  the  expei imenter 
implement  this  program  in  his  own  laboratory.  Finally/ 
a  brief  section  devoted  to  ideas  for  improving  the 
present  program  has  been  included. 

B.  Tunneling 

A  simple  theory  of  single  particle  tunneling  ex¬ 
plains  how  junction  curves  reflect  the  electron  density 
of  states  of  the  metals  comprising  the  junction  (31/10). 
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In  thermal  equilibrium,  the  Fermi  level  of  the 
two  metals  across  the  junction  are  the  same.*  A  volt¬ 
age,  V,  is  applied  so  that  the  left  side  is  at  a 
higher  potential  than  the  right  side.  The  current 
through  tha  junction  can  be  calculated  as  if  the  sin¬ 
gle  particles  were  electrons.  •  he  numfcer  of  electrons 
available  to  flow  across  the  barrier  from  left  to 
right  in  the  energy  interval  IE  will  be  equal  to 
N^(E  -  eV)f(E  -  eV)dE  where  is  the  electron  density 
of  states  on  the  left  side  and  f  is  the  Fermi  function 
which  gives  the  probability  that  a  given  state  of 
energy  (E  -  eV)  is  occupied.  On  the  other  hand,  there 
must  be  unoccupied  states  on  the  right  hand  side  to 
receive  these  electrons.  The  number  of  these  states 
equals  NjiEHl  -  f(E)]  where  ^  is  the  electron  den¬ 
sity  of  states  of  the  metal  on  the  right  hand  side. 
Therefore  the  current  flowing  from  right  to  left  is 

X- )  - 1 v) f (E- e V )] [ - ( (O]]  diE  CO 

where  equals  the  probability  that  an  electron  of 
energy  E  will  tunnel  through  the  barrier  from  left  to 
right . 


See  Fig.  1,  p.  20. 
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By  particle-hole  symmetry  {N (E)  =  N(-E)}  and 
since  f(E)  =  {1  -  f(-E)}  the  current  flowing  from 
right  to  left  is 

where  (E)  is  the  probability  that  an  electron  will 
tunnel  through  the  barrier  from  right  to  left. 

Assuming  that  an  electron  has  just  as  much  chance 
to  go  either  way  across  the  barrier  tor  P^2  (E)  =p21  1 

the  total  current  through  the  barrier,  I,  equals: 

X  *  1,+  Tx=  ^  -  c.V  }  M2(Ob  (c--eV)-i  ( l) 

— 1 

Assuming  P12(E)  is  a  constant»  A,  and  that  for  a  nor¬ 
mal  metal  the  density  of  states  remains  constant  for 
a  small  applied  bias  voltage,  then  (E  -  eV)  =  (0) 

and  N2(E)  =  N2  (0)  where  energies  are  now  measured  from 
the  Fermi  energy.  For  small  applied  V, 

til  Dirac  deWa.  function 

-eV  cktV  <*«.V 

Hence,  in  the  limit  of  small  applied  voltages, 

X  =  AM/O'i  Nx(o)cV  =  0«,V  (.3} 

such  that  D  =  AN1(0)N2(0),  a  constant.  In  other  words, 
the  tunnel  junction  obeys  Ohm's  law.  In  practice, 
this  rule  is  only  approximately  obeyed  (20,17,31)  so 
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that  (dI/dV)N,  the  current  derivative  for  the  "normal 
state  junction",  is  not  a  constant,  as  predicted  by 
Eqn. (3) . 


If  one  of  the  metals  of  the  junction  is  super¬ 
conducting  then  one  normal  metal  density  of  states, 
N(E),  must  be  replaced  by  a  superconducting  density 
of  states.  Using  the  BCS  theory,  the  current  be¬ 


comes  : 


L*  (t-tv)-  hlY3  dlt 


where  the  BCS  density  of  states  is  given  by 


(<0 


MCo-)-  ■  (s') 

and  N(0)  is  the  density  of  the  electrons  of  the  metal 
in  the  normal  state  at  the  Fermi  level  (3) .  Upon 
integration  (10) , 


T*2.D-£-Z(-^ 


rr*  A, 

ACT 


where  is  a  modified  Bessel  function  of  the  second 
kind;  k  is  Boltzman's  Constant;  and  T  is  temperature. 
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It  is  particularly  convenient  to  evaluate  the 
current  Eqn.(*i)  at  T=o«  Here 

f  CO  =  <L  \  $(0-0  t>o 

therefore, 

fCt-e-V)-  $(0  - 

•-0,  O  >  E  t  E  >e.V 

Taking  into  account  that  for  a  superconductor,  due 
to  the  gap,  there  are  no  electron  states  for  0<E<Ag, 
then 


Clearly,  then,  the  conductance  curve  for  the  junction 
does  reflect  the  superconducting  density  of  states.  In 
practice,  however,  because  D  is  not  really  constant  (see 
discussion  page  11)  it  is  usual  to  take  conductance 
curves  for  both  states  of  the  metal,  superconducting 
and  normal.  The  ratio  of  these  two  quantities  at  a 
given  bias  V  is  then: 


-  14  •• 

2f  both  metals  of  the  junction  are  in  the  super¬ 
conducting  state,*  Eqn.  (2)  becomes: 

T=  \  ...  ^Mto)  —  rte-evVHOldit 


4: 


(S.-cV)  E 


(«0 


where  Aqi  is  the  energy  gap  value  of  one  of  the  metals 
and  Aq2  that  of  the  other.  When  T=0  and  A0i=A(p=Aq,  (31) 


X  -  0  Uao  -r  eVK-O  ~  H  *  K  CO]  (  \o") 

A0-t  eV 


and  0  9or  eV  <  2A 


where 


^  _  teV  ~  I  A0\  and  V^C-O  E(-0>  are. 
~  C  eV-vlA6') 


complete  elliptic  integrals  of  the  first  kind.  Then 
dl/dV  at  T=0  is 
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Therefore,  if  both  metals  of  the  junction  are  in  the 
superconducting  state  the  conductance  curve  of  the 
junction,  dl/dV  -  V,  represents  a  convolution*  of  the 
densities  of  states  of  both  metals.  Ir  practice,  be¬ 
cause  D  is  not  a  constant,  ic  is  usual  to  consider  the 
ratio 


m 

mt 


*4  -e.vy- 


-  at 


'i. 

ox 


(\0 


In  principle  if  both  A^  and  A^p 


ire  known  the  normal¬ 


ized  experimental  data 


should  equal 


the  calculated  value  of  this  ratio  as  given  by  Eqn. (11) . 

On  the  other  hand,  the  superconducting  behavior  of 
most  metals  does  deviate  somewhat  from  that  predicted 
by  the  BCS  theory  (28).  Thus,  depending  on  the  state 
of  the  metals  on  both  sides  of  the  junction,  the  exper¬ 
imental  "normal  derivative"  data  will  not  agree  with 
that  calculated  using  either  Eqn. (8)  or  Eqn.  (11).  The 
amount  by  which  the  experimental  and  calculated  values 
differ  is  then  used  as  a  measure  of  how  "weak- coup ling" 
or  BCS-like  the  given  material (6)  are. 


*  In  the  mathematically  rigorous  sense  the  current  in¬ 
tegral  is  not  a  convolution  integral  because  the  upper 
limit  of  the  integral  is  a  function  of  V.  We  will  refer 
to  Eqn.  (11)  as  a  convolution,  hov*e'rer. 
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In  fact,  there  exists  a  whole  class  of  super¬ 
conductors  which  are  not  well-described  by  the  BCS 
theory  of  superconductivity,  the  so-called  "strong- 
coupling"  superconductors.  For  these  materials  the 
superconducting  electronic  density  of  states  is  no 
longer  given  by  Eqn. (5) ,  but  rather  by  (25): 

where  the  energy  gap  parameter  A (E)  is  now  a  complex 
quantity: 

(re) 

and  the  solution  of  the  transcendental  equation 
A (E=Aq)  =  Ag  corresponds  to  the  simple  BCS  energy  gap. 

In  order  to  obtain  phonon  spectra  for  the  mate¬ 
rials  composing  the  tunnel  junction,  it  is  necessary 
first  to  obtain  their  superconducting  densities  of 
states.  This  can  be  done  quite  simply  for  a  normal¬ 
superconducting  junction  using  Eqn.  (8)  and  for  a 
superconducting-superconducting  junction  where  one 
electrode  is  a  BCS-like  material  using  Eqn. (11).* 

*  In  either  of  these  cases,  if  the  "unknown"  is  a 
strong  coupling  material  the  appropriate  superconduct¬ 
ing  .density  of  states  is  replaced  by  Eqn. (12). 
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In  the  latter  case  the  normalized  experimental  conduct¬ 
ance  is  used  along  with  the  BCS  density  of  states  with 
A0  known,  to  "deconvolute"  the  integral  expression  im¬ 
plicit  m  Eqn.(ll)  and  obtain  the  normalized  "unknown” 
density  of  states.  The  deconvolution  process  is  not 
an  easy  one,  but  the  main  procedure  has  previously 
been  outlined  by  McMillan  and  others  (13,15,20). 

Although,  in  essence,  we  follow  this  outline,  a 
different  procedure  for  the  deconvolution  has  been  de¬ 
veloped  in  this  thesis.  Ami  in  particular,  it  can 
also  be  used  when  neither  of  the  junction  electrodes 
is  a  BCS-like  superconductor.  This  modification  has 
been  necessary  for  a  number  of  reasons.  First,  the 
"inversion  program"  as  received  (13)  required  a  rather 
large  density  of  states-deviation"  from  the  BCS  value 
(certainly  well  abovu  1%)  in  order  to  be  able  to  ac¬ 
curately  determine  a  realistic  phonon  spectrum  for  the 
"unknown”  material.  And  second,  many  strong- coupling 
Superconductors  with  small  deviations  from  BCS-behavior 
make  the  best  tunnel  junctions  with  other  non-BCS  super¬ 
conductors.  An  example  of  a  material  where  both  these 
problems  arise  is  Niobium.  It  has  been  shown  that  the 


cleanest  possible  single  crystal  Nb  junctions  are 
made  with  strong-coupling  materials  Indium  (2,9)  and 
the  lead-bismuth  alloy  Pb  ^Bi  ^  (8).  The  expected  de¬ 
viations  in  the  Nb  density  of  states  are  less  than 

1%.* 

The  remainder  of  this  thesis  will  concentrate  on 
th''  mechanics  of  the  Inversion  Program  and  discussions 
of  problem-areas  in  using  the  programs  with  experimen¬ 
tal  data. 


*  M.  A.  Macvica*:,  private  communication. 


Fig.l* 

Energy  diagram  for  a  metal-insulator-metal  junction 

(a)  atthermal  equilibrium 

» 

(b)  at  potential  difference  V. 


^reproduced  from  Ref. (31) 
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'The  energy  diagram  for  an  NS  junction. 

(b)  V  >  A/e,  (c)  I  -  V  characteristic  at 


^Reprinted  from  Ref.  (31) . 


(a)  V  —  0 , 
T  =  0. 
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Fig.  3. 


The  energy  diagram  for  a  SS  junction  at  (a)  eV  =  A^  -  A&, 
and  (b)  eV  =  Afa  +  Aa. 

(c)  The  energy  diagram  for  an  I!S  junction  at  eV  =  A. 

The  I  -  V  character- sties  of  (a)  an  SS  junction 
(b)  NS  junction. 
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CHAPTER  II 

A.  Strong  Coupling  Theory 

The  BCS  theory  is  based  on  an  interaction  be¬ 
tween  pairs  of  conduction  electrons  which  is  trans¬ 
mitted  by  the  ionic  lattice.  The  effective  potential 
associated  with  this  process  is  called  the  electron- 
phonon  interaction,  V.  In  the  BCS  treatment,  it  is 
assumed  that  the  interactions  between  electrons  are 
instantaneous  and  affect  only  those  electrons  within 
an  energy  ED  of  the  Fermi  surface  (ED  is  the  Debye 
energy) .  Thus ,  V  is  taken  to  be  a  constant  for  elec¬ 
tron  energies  les~  than  EQ  and  zero  for  larger  elec¬ 
tron  energies.  The  weak  coupling  condition  is  also 
assumed  in  this  treatment;  i.e.,  that  N(0)V<<1  where 
N(0)  is  the  density  of  states  at  the  Fermi  surface. 

These  assumptions  resulted  in  a  set  of  universal  laws 
which  all  superconductors  obey  (3) .  The  most  impor¬ 
tant  of.  these  is  an  equation  for  the  critical  temper¬ 
ature  where  superconductivity  occurs: 

JkVEoe®  (\3ol) 

Host  real  superconductors,  however,  are  observed 
to  have  properties  which  deviate,  at  least  slightly, 
from  thesvi  laws  (27,28).  In  particular,  it  was  observed 
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in  tunneling  experiments  that  the  BCS  density  of  states 
was  not  adequate  to  explain  the  experimental  conductance. 
Hence ,  a  more  accurate  treatment  of  the  electron-phonon 
interaction  was  sought  which  would  take  into  account 
both  the  non-instantanecus  nature  of  this  interaction 
(phonons  have  finite  velocities)  and  the  fact  that  the 
electronic  energy  states  at  energies  of  the  order  of 
Ep#  are  no  longer  well  defined  (electrons  have  short 
lifetimes  whi_h  in  turn  give  rise  to  damping  effects  in 
their  transport  properties) .  Using  the  thermodynamic 
Greenf function  approach  Eliashberg  (7)  included  both 
of  these  effects  into  a  more  general  theory  of  super¬ 
conductivity.  While  this  approach  is  still  based  on 
the  pairing  concept  introduced  by  BCS,  it  is  more  gen¬ 
eral  in  that  it  relates  the  superconducting  state  direct¬ 
ly  to  details  of  the  normal  state  such  as  the  phonon 
spectrum  of  the  material  and  the  electron-phonon  coup¬ 
ling  constant  (15,24). 

Scalapino,  Schrieffer,  and  Wilkins  (SSW) *  added 
the  residual  Coulomb  interaction  effects  that  screen 
the  electron-phonon  interaction  which  were  also  neglected 
in  the  BCS  theory  (23)  to  the  Eliashberg  treatment  and 
found  the  following  equations  properly  describe  the 
superconducting  state  for  a  one-band  isotropic  model 
at  zero  temperature: 

*  Ref.  (23)  by SSW  contains  most  of  the  following  dis¬ 
cussion;  also  helpful  is  Ref.  (17) - 
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<pii)  [V  -  ]  y,[\CE;t)-/} 


CM) 


[\-  z(fcS]t  =  CdE'Ret' {.E'1- ('4\rf)’']  /a  [^. (£',£)] 


where 


Kt  Ct,E')  =  £_  C  dEc^  (E,,)  Fp  (Eo)  x 

[(e*  +E*"E^ *  iS")  i  Ce'-eyE^vS) 


(\G> 


and.  Ci.(0  =  'He)  -  &_  (E)  + 
1(0  * 


IAT1  e) 


(n) 


Here  E  is  the  electron  energy  measured  from  the  Fermi  sur¬ 
face,  Eg  is  the  phonon  energy,  F  is  the  phonon  density  of 
states  function,  (?  labels  the  phonon  polarization  and  a^(Eq) 
is  an  effective  phonon  coupling  function.  <p  is  the  super¬ 
conducting  electron  pairing  function?  Z  is  the  electron- 
phonon  renormalization  function?  and  A  is  the  energy  gap 
function  analogous  to  the  simple  BCS  energy  gap  value 
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Now  (HE)  is  a  measure  of  the  pairing  self-energy 
of  the  Cooper  pairs ,  it  gives  a  measure  of  how  much  the 
energy  of  the  quasi-particles  (or  excited  single  particle 
states)  is  changed  because  pairing  exists  in  the  system. 

Z(E)  is  a  measure  of  how  much  the  mass  of  the  quasi¬ 
particle  differs  from  the  electronic  mass  because  of  the 
electron-phonon  interactions  in  the  system  (i.e.,  a 
quasi-particle  is  a  "dressed  electron") .  It  is  taken 
to  be  the  same  in  the  conducting  state  as  it  is  in  the 
normal  state.  The  gap  function  A(e)  has  two  parts.  The 
reed  part  giyes  the  amount  of  energy  needed  to  break  a 
Cooper  pair;  i.e.,  the  pair  binding  energy  a?  a  runction 
of  energy.  The  imaginary  part  corresponds  to  the  shorten¬ 
ing  of  the  lifetimes  of  the  excited  quasi-particles;  i.e., 
it  is  proportional  to  the  probability  that  a  virtual  excited 
guasi^-particle  will  decay  from  an  energy  E  to  the  gap  edge 
Aq  with  the  emission  of  a  phonon. 

Roughly,  the  physics  of  the  pairing  process  is  that 
a  spin-up  electron  at  the  Fermi  surface  exchanges  a  phonon 
with  a  spin-down  electron  also  at  the  Fermi-surface  with  a 
momentum  transfer  "q".  a(E)  is  a  measure  of  the  "strength 
of  the  interaction"  or  probability  that  this  exchange  will 
be  made.  F(E)  indicates  the  kind  of  "particle"  (a  phonon) 
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that  will  be  exchanged. * 

Instead  of  using  an  additional  self-consistent 
equation  which  would  take  care  of  the  Coulomb  repulsion 
between  electrons,  SSW  used  a  Coulomb  pseudo-potential 
method  that  accounts  for  higher  energy  excitations 
(E>EC)  by  cutting  off  their  integrals  at  E=EC  where 
EcslOEjj  and  inserting  a  dimensionless  parameter  y*  = 
N(0)Uc  into  their  equations.  This  factor  is  a  number 

i 

and  not  a  function  of  energy  because  the  effective 
Coulomb  potential  varies  on  a  scale  of  the  order  of  1 
eV  while  the  pairing  potential  varies  only  in  an  energy 
interval  of  meV's  about  the  Fermi  energy.  The  effect 
of  the  Coulomb  term  on  the  superconducting  system  is, 
basically,  to  reduce  the  ability  of  the  phonons  to  re¬ 
spond  to  the  electrons;  and  thus  to  reduce  the  contribu¬ 
tion  of  the  phonons  to  the  electron's  self  energy. 

SSW  were  also  able  to  show  that  the  density  of 
states  of  a  superconductor  is  (23,25) : 

Ws(t>N<cO*e  ?}  Cifl 

*  If  the  phonon  is  emitted  when  the  first  electron  goes 
from  state  then  it  must  be  absorbed  by  the  otj^er 

partner  of_>th§  pair  which  will  scatter  from  -k^  to  -k+ ; 
and  q  is  (k'-k) .  We  shall  refer  to  either  of  the  two 
electrons  mentioned  as  the  virtual  excited  quasi-particles 
or  electrons. 
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which  agree* with  the  BCS  expression  for  the  density  of 
states  (Eqn.(5))  in  the  limit  of  an  energy  independent 
gap  function.  Although  this  expression  for  the  normal¬ 
ized  density  of  states  (NDS)  mekes  clear  that  structure 
in  A(E)  will  be  reflected  in  the  NDS,  this  fact  is  more 
clearly  brought  o".t  by  expanding  Ng(E)  in  powers  of 
(A/E)  : 

E  »  mb')  2 

There  exists  a  simple  model  (23,24)  for  describ¬ 
ing  how  the  solution  of  Eqns. (14) - (17)  relate  to  the 
"picture"  we  have  used.  Consider  the  case  where  the 
material  has  a  phonon  density  of  states  with  a  single 
particle  peak  at  Eq  (See  Pig.  4.).  At  excitation 
energies  less  than  EQ,  since  most  of  the  phonons  which 
can  be  exchanged  by  the  electrons  have  a  higher  natural 
frequency  than  Eq  they  tend  to  over-respond  to  the 
•driving  frequency"  causing  an  attractive  electron- 
electron  interaction  while  at  energies  greater  than  E 

o 

the  phonons  are  driven  above  the  natural  frequencies  and 
tend  to  respond  out  of  phase  with  the  driving  field.  Thi 
corresponds  co  a  repulsive  electron-electron  interaction. 
For  excitation  energies  near  Eg,  the  real  part  of  A  peaks , 
or, the  pair  binding  energy  increases  strongly  because  of 
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Fig.  4 . * 

(a)  A  phonon  peak  at  (to  =  a)Q)  causes  A(u)  to  look  as 

% 

(b)  and  produces  density  of  states  (c) . 


♦reprinted  from  Ref.  (11) . 
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the  resonant  exchange  of  phonons  that  have  the  same  re¬ 
sponse  frequency  as  that  of  the  real  lattice.  The  Re^ffc) 
decreases  to  negative  values  for  E>E^  since  the  electron- 
electron  interaction  is  repulsive.  The  ImA(E)  which  is 
a  measure  of  the  damping  effects  will  be  zero  at  very 
low  energies  since  the  quasi-particles  are  well  defined 
there;  rise  to  a  peak  near  the  resonance  because  of  the 
increase  in  phonon  emission;  and  finally  it  will  vanish 
because  the  excited  electrons  can  no  longer  fall  to  the 
gap  edge  by  omitting  a  phonon. 

By  considering  Eqn. (19)  for  the  density  of  states 
it  is  clear  how  the  phonon  structure  is  indirectly  mir¬ 
rored  in  Ns (E)  through  the  gap  function  and  why  there  are 
shoulders  and  drops  in  Ns (E)  with  respect  to  the  BCS-NDS. 
Adding  a  Coulomb  term  to  this  picture  would  only  change, 
qualitatively,  the  ReA(E)  after  resonance.  In  this  case, 
it  would  remain  negative  instead  of  going  to  zerp,  in¬ 
dicating  the  net  Coulomb  repulsion  of  the  virtual  unbound 
electrons  (23,24} . 

Another  important  quantity  of  interest  is  the 
effective  electron-phonon  coupling  constant, X  ,  which 
equals  (16,23) 

\  -  z  J  (26) 

o 
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and  represents  some  "averaged"  electron-phonon  parameter 
for  the  whole  system.  In  fact,  the  renormalization 
function,  Z(E',  near  the  gap  edge  can  be  shown  to  be 
equal  to : 

z  -  \  +  x  -  UO 

tA 


where  m*  is  the  effective  electron  mass.  This  expression 
brakes  sense  since  X  is  a  measure  of  the  phonon  screening 
,'of  electrons . 

McMillan  (16)  has  given  a  very  clear  discussion  on 
how,  once  X,  Z(E),  and  a  F  are  known,  many  other  elec¬ 
tronic  properties  of  the  normal  state  of  metals  can  be  de 
termined  su«-h  as  the  Fermi  velocity,  the  electronic  co¬ 
efficient  of  the  specific  heat,  the  electronic  band  struc 
ture,  and  the  normal  metal  density  of  states.  It  is 
through  a  knowledge  of  such  properties  that  new  uses  for 
these  materials  can  be  found.  These  normal  state  prop¬ 
erties  can  also  be  used  to  predict  what  the  critical 
temperature  for  the  superconducting  state  will  be. 
McMillan  (16)  has  derived  the  formula: 


.\.qa(u>D 

X  V  £  -MV 


(22) 


When  compared  with  the  BCS  Eqn. (13a)  it  is  clear  that 
the  "effective  interaction"  for  superconductivity  be¬ 
comes:  x  C'a.'B) 

just  as  our  crude  picture  would  predict. 
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Using  the  fact  that  A(E)  is  reflected  in  the 

2 

density  of  states  of  superconductor  and  that  a  (E)F(E) 
is  related  through  the  gap  equations  to  Ng(E)  Me'  illan 
and  Rowell  (17,13)  developed  a  program  no  determine 
the  ct^F  of  a  metal  from  the  experimental  first  deriva¬ 
tive  curve  of  a  tunnel  junction.  The  method  by  which 
this  "inversion" is  accomplished  is  described  in  Chapter  3. 

B.  Practical  Details 

It  would  be  instructive  to  point  out  several 

2  2 

properties  about  dl/dV  and  d  I/dV  curves  and  the  junc¬ 
tion  samples  which  influences  the  correctness  of  the  in¬ 
put  data  and,  thus,  the  correctness  of  the  output  of 
this  program. 

First,  the  junctions  used  must  be  good  ones. 

Rowell  (17,20)  has  laid  down,  in  several  places,  the 
basic  criteria  for  determining  the  perfection  of  junc¬ 
tions.  Rowell  and  MacVicar,  et  al. ,  have  laid  down 
rules  for  determining  if  the  dl/dV  curves  of  a  given 
junction  do  reflect  bulk  material  characteristics  (17,9). 

From  the  BCS-NDS  given. in  Eqn. (5)  we  see  that  as 
E  approaches  infinity  the  BCS-NDS  approaches  one*.  Even 
for  non-BCS  superconductors  for  larger  bias  the  NDS 
equals  the  BCS-NDS.  Therefore,  for  higher  voltages  the 
vaJue  of  the  normal  first  derivative  should  equal  the 
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value  of  the  superconducting  first  derivative.  It  has 
been  found,  however,  that  this  does  not  always  happen 
(12,11).  A  possible  explanation  of  this  phenomenon  is 
that  in  the  normal  state  the  resistance  of  the  metal 
films  of  the  junction  is  significant  in  comparison  with 
the  junction  resistance  (12) .  To  compensate  for  errors 
of  this  type,  in  the  past,  a  constant  has  been  added  to 
all  values  of  the  normal  first  derivative  data  such  that 
this  data  is  made  equal  to  the  superconducting  data  at 
high  biases.  The  phonon  curves  generated  from  this  data 
have  been  unlike  the  phonon  curves  for  other  metals  (17) . 
The  effect  of  such  a  procedure  on  the  generated  a^F 
spectrum  is  unknown. 

McMillan's  program  requires  only  that  the  NDS  of 
a  metal  be  known  in  order  to  determine  its  phonon  spec¬ 
trum.  Still,  it  is  preferable  to  modify  the  first  deriva¬ 
tive  data  by  the  second  derivative  data.  The  reason  for 
this  is  that  critical  points  in  the  phonon  spectrum  of  a 
metal  (the  places  where  (dF(2)/dE)  equals  infinity)  are 
reflected  as  square  root  singularities  and  logarithmic 
singularities  in  the  d^I/dV^  curve  (21).  However,  due  to 
thermal  smearing  and  phonon  broadening,  singularities 

appear  as  peaks  and  valleys  rather  than  mathemat: cally 

U.'k'tc 

sha^'p  singularities.  The  point  is  thatAthe  conductance 
curve  contains  very  accurate  data,  the  second  deriva- 

tive  curve, resolves  the  "small"  structure  in  the  first 
derivative  data. 
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Finally,  it  should  be  pointed  out  that  the 

McMillan  inversion  gives  a*F,  the  phonon  spectrum 

multiplied  by  the  electron-phonon  coupling  function, 
o 

a  F  should,  nevertheless,  bo.  a  good  reflection  of. 

2 

F(E)  since  a  (E)  is  thought  to  be  a  slowly  changing 

at 

continuous  function  with  no  critical  points.  Thus  one 

2 

expects  the  "general  shape"  of  a  F(E)  to  be  that  of 

2 

F(E),  and  the  critical  points  of  a  F(E)  to  be  those  of 
F(E) . 

Neutron  scattering  is  an  alternative  way  of 
determining  a  phonon  spectrum  of  a  metal.  Comparisons 
of  the  two  methods  for  various  metals  have  shown  that 
the  critical  points  of  the  phonon  spectrums  are  in  gen¬ 
eral  agreement  (17) .  The  discrepencies  in  the  magnitude 
of  the  peaks  can  be  accounted  for  by  the  fact  that 
tunneling  inversion  generates  a^(E)F(E)  while  neutron 
scattering  is  thought  to  generate  F(E). 

C.  Temperature  Dependence (23) 

Finite  temperature  dependence  has  also  been  intro¬ 
duced  into  the  strong  coupling  theory  of  superconductiv¬ 
ity  iEqns. (14) - (17)  become  very  much  more  complicated.]. 
Here  we  merely  want  to  point  out  that  qualitative  changes 
that  temperature  brings  into  our  crude  picture  of  the 
superconducting  state.  Increasing  T  is  equivalent  to 
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putting  thermal  energy  into  cur  system.  This  thermal 
energy  will  tend  to  weaken  or  break  pairs  which  are  not 
tightly  bound.  Thus,  just  as  in  theBCS  theory,  the 
gap  edge  function  Aq  becomes  a  function  of  temperature, 
decreasing  as  T  increases.  In  an  analogous  way  we 
expect  the  entire  real  part  of  A(E)  to  be  reduced  from 
the  T=0  value  of  A(E).  In  addition  the  virtual  excited 
quasi-particles  in  this  state  can  now  recombine  with  the 
'  thermally  generated  quasi-particles ,  so  we  expect  more 
structure  (or  resonances)  in  both  real  and  imaginary 
parts  of  the  gap  functions .  Clearly ,  the  damping  effects , 
reflected  in  the  ImA(E)  will  also  be  decreased,  energy  for 
energy,  over  the  zero  temperature  case.  Finally,  it  has 
been  shown  that  (24) 

&fE,T)aL^r]  &(e.,C>)  (.2-0 

that  is,  the  temperature  dependent  gap  function  is, 
roughly,  a  scaled  down  version  of  the  A(t’)  curve  where 
the  scale  factor  is  the  BCS  temperature  dependent  gap 
edge. 
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Chapter  III 


"The  Programs" 

Three  programs  are  required  to  produce  a  phonon 
spectrum  from  the  raw  tunneling  data.  Stage  I  smoothes 
and  interpolates  the  experimental  input.  Stage  II  re¬ 
moves  the  influence  of  the  counter-eiectrode  and  produces 
data  characteristic  only  of  the  unknown  metal.  Stage 
III  manipulates  theoretical  functions  which  should  re¬ 
produce  the  output  of  Stage  II.  A  detailed  discussion  of 
each  stage  is  given  below. 

Stage  I 

The  goal  of  the  first  program.  Stage  I,  is  to  take 
the  raw  data  and  compute  the  deviation  of  the  experimental 
conductance  values  of  the  ju  ction  from  tKat  predicted 
by  the  BCS  tunneling  model  iGiven  in  Egn.(ll)].  The 
quantity  a (V)  determined  at  all  applied  voltages,  V, 
where 


O  ='  (normalized  experimental  data) _  _  ^ 

"(BCS  predicted  normalized  function) 

> 

is  used  as  the  measure  of  that  deviation.  (One  is  sub¬ 
tracted  from  the  quotient  to  emphasize  the  deviations 
which  are  of  the  order  of  6%,  at  most  (17,  24))  # 


I 


Stage  I  proceeds  in  the  following  steps.  dl/dV 
with  both  metals  superconducting,  dl/dV  with  both  metals 
normal,  and  d2I/dV2  with  both  metals  superconducting 
are  computed  from  the  input  data  which  is  taken- directly 
from  vile  xy-recorder.  Then  the  normalized  density  of 
states,  (dI/«3V)supER/(dx/dv)N0RWAL  is  computed.  If 
second  derivative  data  is  available,  the  i-tegral  of 
the  second  derivative  data  is  least  square  fitted  to  the 
normalized  first  derivative  data,  and  this  fitted  norma¬ 
lized  first  derivative  is  smoothed.  Then,  the  deviation 
of  this  experimental  normalized  first  derivative  from 
the  BCS  predicted  normalized  first  derivative  is  computed. 

Several  steps  are  required  to  compute  (dI/dV)crTO 

SUPER 

and  (dI/dV)NQRMAL  from  the  laboratory  input: 

The  input  is  given  in  pairs,  an  X  value  and  a  corres¬ 
ponding  Y  value.  It  would  be  preferable  if  X  equalled 
a  voltage  value,  V,  and  Y  equalled  dl/dV  at  V;  instead, 
X^aV+b  and  Y=c(dI/dV) +d.  Therefore  the  actual  values  of 
V  and  dl/dv  must  be  computed.  To  do  this  several  para¬ 
meters  must  be  supplied  to  the  program  to  compute  a,b, 
c,  and  d. 

Secondly,  there  may  exist  in  the  input  stream  several 
anomalous  prints  which  repeat  previously  recorded  input 
pairs;  that  is,  the  X  values  of  these  points  is  le^s  than 
a  previous  X  value  of  tne  input  stream.  The  program 
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assumes  a  monotonically  increasing  X  value;  hence,  these 
anomalous  points  are  removed. 

Thirdly,  the  values  of  dl/dV  will  not  usually  be 
recorded  experimentally  at  a  specific  voltage  interval, 
i.e,  at  V=0,H,2H, . . . ,nH  with  H  being  the  interval.  How¬ 
ever,  the  programs  Stage  II  and  Stage  III  require  that 
<r(V) ,  and  therefore  dl/dV  values,  be  tuken  at  constant 
intervals.  Therefore,  the  values  of  dl/dV  for  specified 
increments  of  the  bias  voltage,  V  must  be  computed.  This 
is  done  either  by  polynomial  interpolation  using  the 
points  whose  voltage  values  are  in  the  vicinity  of  the 
voltage  at  which  the  values  of  dl/dV  are  required,  or  by 
linear  or  quadratic  least  square  fitting.  The  least  square 
fitting  determines  either  a  first  or  second  degree  poly¬ 
nomial  to  represent  the  input  data.  Polynomial  inter¬ 
polation  then  is  used  to  determine  dl/dV  at  V^,V2,V3, . . . , 
etc.  The  advantage  of  the  linear  or  quadratic  fitting 
option  is  that  it  tends  to  eliminate  the  random  noise 
found  in  the  input. 

Stage  I  was  designed  modularly  so  that  it  can  be 
easily  changed  if  different  experimenters  should  have 
slightly  different  requirements.  In  addition  to  this, 
an  effort  was  made  to  put  in  a  variety  of  features  for 
extra  flexibility.  Among  the  features  included  are  a 
choice  of  either  polynomial  interpolation  using  from  1 
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to  10  points  or  linear  or  quadratic  fitting  using  from 
1  to  100  points;  the  ability  to  measure  dl/dV  with  V 
positive  or  negative  I(dI/dV)  is  symmetric  around  V=0] ; 
and  the  ability  to  find  the  normalized  first  derivative 
even  though  the  second  derivative  curve  is  unavailable. 


Stage  II 


r; 


Stage  II  is  a  generalized  deconvolution  routine. 

As  mentioned  in  Chapter  I  if  both  metals  are  supercon¬ 
ducting  the  Formalized  first  derivative  equals 

( a  r 
Hrf  ~  77  j 
V  dv  )N  .  *** 

WP 

where  is  the  normalized  density  of  states  of  metal  I 
and  A01  is  its  energy  gap  value  and  D2  and  A02  are  the 
analogous  quantities  for  metal  II. 

.  When  the  phonon  spectrum  of  some  material  is  desired, 
it  is  not  always  possible  to  fabricate  a  tunnel  junction 
where  the  counter-electrode  remains  normal  at  the  tem¬ 
perature  where  the  experiment  must  be  done.  Instead, 
the  counter-electrode  becomes  superconducting;  if  it 
is  a  BCS-like  super  conductor  then  Eqn. (11)  can  be  used, 
With  little  difficulty,  to  obtain  the  unknown  electronic 
density  of  states.  If,  on  the  other  hand,  the  counter¬ 
electrode  is  non-BCS,  it  introduces  its  own  devia¬ 
tions  into  the  normalized  experimental  conductance 
data.  Therefore,  the  effect  of  the  counter-electrode 
must  be  eliminated  from  the  experimental  data 

I26< 


-43- 


/ 

at  the  same  tiro'-!  that  the  deconvolution  of  Eqn.  (11) 
is  being  accomplished.  To  deconvolute  the  experimental 
data  Stage  II  requires  two  sets  of  inputs  s  the  smoothed 
experimental  data  from  Stage  I  and  the  normalized  first 
derivative  of  the  "given"  counter-electrode. 

The  exact  workings  of  the  program  are  explained  in 
the  appendix.  Only  one  problem  involved  in  handling  the 
input  data  will  be  discussed  here:  the  problem  revolves 
about  the  experimental  fact  that  it  is  difficult  to 
determine  accurately  junction  conductances  very  close  to 
the  gap  edge  (13) .  It  is  usual  to  assume  that  the  junc¬ 
tion  conductance  equals  the  value  predicted  by  the  BCS 
theory  at  these  biases.  This  is  not  really  too  gross 
an  assumption  since  the  BCS  function  is  rising  rapidly 
towards  its  singularity  for  V=A0,  and  the  deviations  due 
to  strong-coupling  effects  are  far  smaller  than  those 
due  to  the  finite  rounding  off  of  the  BCS  function.  This 
means,  however,  that  the  first  non-zero  experimental 
value  in.Ainput  data  suddenly  has  a  deviation  from  the 
BCS  value  while  the  preceding  point  does  not.  Theoreti¬ 
cally  the  deviation  from  BCS  should  start  at  zero  and 
continuously  increase  for  very  low  offset  voltages. 

The  lack  of  experimental  information  at  low  offset 
voltages  introduces  a  large  error  into  the  deconvolution 
routine  especially  for  the  first  data  points.  The  computer 
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sub-routine  is  written  so  that  the  first  few  "input" 
data  slots  are  zero.  The  deviation  of  the  first  non¬ 
trivial  experimental  input  should  be  looked  at  by  the 
program  operator  and  "data"  inserted  in  the  "empty 
slots"  which  takes  the  deviations  from  0  to  the  first 
experimental  deviation  in  a  smooth  way.  So,  for  example, 
if  the  first  deviation  is  .0035,  the  "data  points"  .001, 
.002,  .003  migi.t  be  inserted. 

The  normalized  density  of  states  (NDS)  for  the  known 
counter-electrode  ("given")  can  be  determined  in  one  of 
two  ways.  The  first  method  is  to  make  two  tunneling 
junctions  simultaneously  with  films  of  the  "given":  one 
with  the  metal  of  interest  and  the  second  with  a  BCS  or 
normal  metal.  The  NDS  of  the  given  material  can  then  be 
determined  using  Stages  I  and  II  with  tunneling  data  from 
the  BCS  or  normal  metal  junction.  This  computed  NDS  of 
the  given  is  then  simply  "read-in"  Stage  II  when  the 
tunneling  data  from  the  junction  containing  the  material 
Of  interest  is  processed  by  the  inversion  program.  Since 
the  two  junctions  were  processed  simultaneously,  and  thus 
under  the  exact  same  conditions,  one  would  hope  that  the 
NDS  is  the  same  for  both  junctions.  Needless  to  say, 
this  may  not  be  true.  However,  by  running  several  dif¬ 
ferent  samples,  or  junctions,  with  the  material  of  interest 
(all  either  processed  at  the  same  time  as  the  "reference 
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junction"  with  the  non-BCS  counter-electrode,  or  each 
made  separately  and  each  having  its  own  "reference 
junction")  a  close  check  can  be  made  on  the  validity 
of  the  procedure  by  comparing  the  NDS  for  the  "unknown" 
from  the  different  runs. 

For  the  case  where  it  is  not  possible  to  prepare 
simultaneous  "reference  junctions"  another  option  is 
provided  in  the  program.  It  consists  of  reading  in  a 
phonon  spectrum  for  the  given  material  and  from  that 
phonon  spectrum  computing  the  NDS  of  the  given.  This 
routine  has  been  taken  directly  from  Stage  III  and  pre¬ 
sumes  that  the  "given"  material  is  sufficiently  well- 
characterized  that  its  phonon  spectrum  is  easily  obtained 
and  is  not  sensitive  to  material  changes  inherent  in  the 
present  method  of  making  the  tunnel  junction.  The  only 
check  on  the  accuracy  of  the  unknown’s  NDS  and  phonon 
spectrum  (calculated  in  Stage  III)  is  the  ability  to 
reproduce  the  results  from  sample  to  sample.  This  can, 
of  course,  be  made  quite  rigorous  by  using  different 
counter-electrodes  for  the  given  material.  Agreement  of 
the  calculated  NDS  and  phonon  spectrum  for  the  unknown 
material  under  such  circumstances ,  makes  it  almost  certain 
that  the  basic  material  properties  have  been  calculated. 
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Stage  III 

Stage  III  is  the  original  McMillan- Rowe 11  Inversion 
Program  Stage  III  (13) .  The  main  idea  in  this  program  is 
to  calculate  the  normalized  electronic  density  of  states 
for  the  unknown  and  compare  these  values  with  those 
determined  in  Stage  II.  If  these  two  data  arrays  agree,, 
then  it  is  assumed  that  the  phonon  spectrum  used  to  cal¬ 
culate  the  data  of  Stage  III  is  that  characterizing  the 
unknown.  Although  this  may  seem  like  a  simple  enough 
procedure,  it  is  not... it  hinges  on  knowing  the  approxi¬ 
mate  phonon  spectrum  and/or  knowing  "how"  to  ctenge  the 
assumed  phonon  spectrum  to  bring  the  calculate*  i  electronic 
HDS  of  Stage  III  into  agreement  with  that  derived  from 
Jin  -raw  tunneling  data;  i.e.,  the  output  of  Stage  II. 

The  functions  that  are  calculated  and  plotted  in 
this  stage  were  described  previously  in  Chapter  II.  In 
this  section,  only  the  general  mechanics  of  this  stage 
will  be  discussed,  '.he  program  is  written  to  operate  in 
either  of  two  modes :  (a)  the  approximate  form  of  the 

phcnon  spectrum  is  known  and  (b)  it  is  not  known.  Only 
mode  (b)  will  be  described,  since  during  the  iterations 
of  the  program,  in  this  instance,  the  phonon  spectrum 
that  is  directly  "read-into"  the  program  in  mode  (a)  is 
calculated  and  subsequently  "improved"  or  iterated  to 
obtain  the  final  phonon  spectrum.  Clearly,  if  a  good 
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approximation  to  the  phonon  spectrum  is  already  known, 
the  cost  of  running  this  program  is  substantially  re¬ 
duced.  It  should  be  noted,  therefore,  that  in  an  experi¬ 
mental  investigation  once  an  <x2(e)F(E)  spectrum  is  obtained 
fcr  one  sample  (using  mode  (b) ,  say)  then  this  spectrum 
should  be  used  in  Stage  III  operated  in  mode  (a)  for 
other  samples.  If,  in  fact,  the  a2F  spectrum  calculated 
for  a  single  sample . is  characteristic  of  the  material, 
then  the  changes  of  spectrum  from  sample  to  sample  should 
be  within  the  "error-bar? . "  of  the  experimental  data  set, 
itself.  Given  that  this  occurs,  the  investigation  is 
completed. 

Mode  (b)  of  Stage  III  operates  in  the  following  way.* 
An  iteration  process  is  used  by  McMillan's  inversion  pro¬ 
gram  to  solve  the  gap  equations .  Starting  zeroth  order 
values  must  be  chosen  for  y*,  A (E) ,  and  a2F(E).  y*(0) 

the  first  value  chosen  for  y*  is  usually  set  equal  to 
.11  and  a2F(E) can  be  an  almost  arbitrary  function 
defined  up  to  E-Eq,  the  maximum  energy  of  the  phonon 
spectrum.  A(0)(E),  the  2eroth  order  value,  is  taken 
to  be  equal  to  Aq  for  E<  Eg,  and  equal  to  zero  for 
E  >  Eg,  where  Ec~3Eq  and  Aq=A(E=Aq),  the  BCS  gap  value. 


♦The  following  discussion  has  been  adapted  from  Ref.  C13)-, 
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Inserting  these  starting  values  in  the  right  hand 
side  of  Eqn. (14)  and  Eqn.  (16)  will  then  give  a  calcu¬ 
lated  value  for  the  gap  function,  AlD(E).  Reinsert¬ 
ing  A (E)  in  the  equations,  with  the  same  y*  and 
<*2F,  will  give  a  second  calculated  value  A  ^  (E) .  This 
process  is  continued  until  A ^ (E)  converges.  Then  the 
density  of  states,  Nq(E) ,  can  be  calculated  using 
Eqn« (12) * 

(  It  is  now  necessary  to  correct  the  starting  a2F 
and  y*  based  on  the  difference  between  Nc,  the  calcu¬ 
lated  density  of  states,  and  the  input  data  from  Stage  II, 
Ne  tha  experimental  density  of  states.  17a  use  linear 
feedback  techniques  (17)  to  correct  a^F  for  the  error 
lNe(E)-Nc(E)] .  This  is  done  by  first  calculating  the 
effect  on  bi(E),  for  all  data  points,  of  a  small  change  in 
a2F;  i.e., 

SL nU'H  / CO]  . 

Then 

a*af(^)  ---<X  F  (  , 

where 


£L*<xF(fc<*V]  = 


I 

I 


132< 


19 


i 

i 

I 

where  N^®)  is  the  first  value  calculated  for  Nq  and 

the  functional  derivative  is  calculated  from  the  ana- 

2 

lytic  form  of  the  gap  equation.  The  iteration  of  a  F, 

namely,  solving  for  A(E)  and  NDS  from  a^F  and  changing 

a2F  depending  on  how  the  calculated  NDS  compares  to 

2 

the  experimental  NDS  is  continued  until  a  F  converges. 
Thrt  is,  since  the  gap  equation  is  non-linear  the  total 
process  must  be  iterated  until  a^F  converges.  y*(n)  is 
'determined  by  the  requirement  that  A(E»Ag)  always  equals 
the  BCS  value,  AQ. 

Note,  however,  that  because  the  data  does  not  ex¬ 
tend  to  the  gap  edge,  it  is  necessary  to  make  an  assump¬ 
tion  regarding  the  energy  dependence  of  a2F  at  low  ener¬ 
gies.  The  program  assumes  that  a2F  i r  : raportional  to 
the  square  of  the  energy  up  to  some  prespecified  value 
of  energy. 

It  is  appropriate  here  to  give  a  few  guidelines  on 
the  operation  of  this  stage  of  the  inversion  program. 

A  "normal"  convergence  run  will  have  A  iterated  approx¬ 
imately  5  times. and  ct2F  approximately  8  times.  If  the 
form  of  a^F  is  known  from  previous  samples,  the  "new- 
sample"  might  r'id  to  i cerate  as  many  as  4  times. 

Notice  that  it  is  impossible  for  a2F  to  converge  while 
the  experimental  density  of  states  in  not  close  to  the 
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calculated  density  of  states.  This  is  because  it  is  the 
difference  betv:een  the  two  that  determines  how  much 
a  F  will  change;  a  big  difference,  Ne“Nc»  implies  a 
large,  or  significant  change  in  a?F. 

(d2I/dV2)  data  is  used  to  resolve  the  (dl/dV) 
data  because  the  critical  points  in  a^F  should  appear 
as  singularities  in  (d2I/dV2) .  On  the  other  hand,  the 
singularities  in  this  curve  are  smeared  due  to  thermal 
excitations  and  the  most  sensitive  way  to  determine  if 
the  critical  points  of  a2F  agree  with  those  in  (d2I/dV2) 
is  the  human  eye.  Although  the  computer  could  be  used 
for  this  cross-checking  procedure,  it  would  not  be  very 
efficient. 

It  is  also  important  to  know  if  changes  in  a2F(E) 
correspond  to  experimental  uncertainty  or  to  real  differ¬ 
ences  in  the  experimental  density  of  states.  A  measure 
of  precision  of  a  F  can  be  obtained  by  running  the  in¬ 
version  program  with  the  experimental  density  of  states 
at  their  greatest  and  least  values  within  the  experi¬ 
mental  error  and  thus  determining  the  maximum  change  in 
a2F(E) .  Any  change  less  than  this  is  certainly  inherent 
in  the  experi  *»nt,  itself.  'This  assumes,  of  course, 
that  both  sets  of  data  are  such  as  to  produce  a  converging 
o2F)  . 
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COSTS 

The  three  stages  were  run  on  the  IBM  170-155  at 
the  MIT  Information  Processing  Center.  The  costs  of 
running  each  of  the  three  programs  can  be  kept  low  by 
running  them  under  low  priority  and  by  using  object 
decks.  The  range  of  costs  are: 

Stage  I  $  6  -  $10 

Stage  II  $  2  -  $  4 

Stage  III 

(McMillan  program  alone)  $  8  -  $12 
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CHAPTER  IV 


"Future  Work" 

The  previously  mentioned  programs  can  certainly 
be  improved.  The  first  improvement  is  to  put  temper¬ 
ature  dependence  into  the  convolution  program.  The 
theory  of  strong-coupling  superconductivity  at  finite 
temperatures  has  been  completely  worked  out  and  equa¬ 
tions  analogous  to  those  used  in  the  zero  temperature 
case  are  available  (22)  .  With  these  equations  in  hand, 
it  does  not  seem  as  if  in  principle,  there  will  be 
problems  in  addition  to  those  already  treated,  to  be 
solved.  However,  until  such  a  project  is  undertaken 
it  is  impossible  to  estimate  what  difficulties  could 
be  encountered. 

Secondly,  it  would  be  interesting  to  consider  the 
effects  of  various  choices  of  the  parameter  y*  on  the 
final  a^(E)F(E)  spectrum  selected.  In  some  recent 
tunneling  studies  (29,14)  it  has  been  found  that  once 
y*  is  chosen  to  give  a  phonon  spectrum  which  duplicates 
the  experimental  input  data,  just  by  scaling  up  (in 
amplitude)  a^F,  a  different  y*  which  is  in  better 
agreement  with  the  "expected"  y*  can  be  chosen  that 
will  also  duplicate  the  experimental  input  data.  This 
behavior  is  NOT  understood  and  the  reason  for  this 
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behavior  should  be  determined.  Understanding  the 
mechanism  behind  it  might  enable  the  theoretician 
to  determine  if,  in  these  cases,  it  really  is  totally 
the  electron-phcnon  mechanism  which  produces  super¬ 
conductivity  (30).  Or,  more  pessimistically,  if  the 
parameter  p*  really  is  a  sensitive  superconducting 
parameter  (6) . 

Another  improvement  would  be  to  arrange  Stage  II 
and  Stage  III  in  such  a  way  that  a  phonon  spectrum 
for  the  "given"  electrode  would  be  used  to  calculate 
its  super conducting  electronic  density  of  states 
(this  ii  fed  into  Stage  II) .  Then  Stage  III  would 
iterate  the  phonon  spectrum  of  the  "unknown"  until 
the  "best  possible  value"  for  its  calculated  density 
of  states  is  obtained.  If  total  convergence  with 
experimental  data  has  not  been  achieved  with  this 
"best  possible  value",  small  changes  would  be  made 
in  the  phonon  spectrum  of  the  given,  and  a  new  it¬ 
eration  on  the  "unknown"  phonon  spectrum  carried  out. 
This  is  done  repetitively  until  McMillan’s  program 
converges.  Unless  the  "given"  counter-electrode  is  a 
well-characterized  material,  however,  such  a  pro¬ 
cedure  might  be  prohibitive  from  a  cost-efficiency 
point  of  view.  In  addition,  the  problem  of  the 
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accuracy  of  the  given's  normalized  density  of  states 
(NDS)  remains.  That  is,  if  this  parameter  does  NOT 
truly  reflect  the  density  of  states  of  the  counter- 
alectrode,  even  though  the  inversion  program  con¬ 
verges,  there  is  no  assurance  that  the  phonon  spec¬ 
trum  or  electron-phonon  coupling  constant  obtained 
for  the  "unknown"  material  accurately  reflect  its 
basic  properties. 

This  brings  us  to  the  final  recommendation  of 
this  thesis.  A  careful  and  quantitative  study  of 
the  effects  of  errors  in  the  "given"  NDS  on  the 
final  phonon  spectrum  for  the  "unknown"  material 
should  be  carried  out.  In  particular,  it  should 
be  determined  if  Stage  III  will  converge  when  this 
NDS  in  not  accurate.  In  the  event  that  it  would, 
such  that  several  tunnel  junctions  composed  of  the 
same  two  materials  all  give  identically  the  same 
"unknown"  a  (E)F(E)  (with  the  wrong  given),  the  pres¬ 
ent  method  of  deconvoluting  the  tunneling  data  could 
not  be  used  in  the  case  that  both  junction  electrodes 
are  strong-couplers.  The  best  theoretical  estimates* 

*  J.  Bostock,  private  communication 
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at  the  present  time  indicate  that  this  will  not  be 
true.  Hence,  we  expect  that  implementing  a  "double 
deconvolution"  is,  in  principal,  quite  possible  and 
should  be  done  in  the  near  future. 
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STAGE  I 

STAGE  I  takes  experimental  data  from  the  laboratory 
setup  and  converts  it  to  the  following  form: 


1  -  cr 

Details  will  be  given  in  comment  cards  within  the  pro¬ 
gram  and  in  the  introductions  to  the  various  subroutines. 

INPUT 

TITLE  an  array  of  80  characters  to  be  printed  which 
gives  a  title  or  contains  information  such  as 
date,  metals,  junction  number,  and  any  comments 
on  the  data  set. 

H  is  the  interval  in  millivolts  at  which  values 

for  the  normalized  density  of  states  is  out¬ 
put.  The  default  value  for  H  =  .2  millivolts. 

MAXX  the  maximum  desired  energy  at  which  the  nor¬ 

malized  density  of  states  (NDS)  and  a  is  to 
be  computed.  If  MAXX  is  greater  than  the 
maximum  energy,  E,  for  which  there  exists  in¬ 
put,  MAXX  will  be  set  equal  to  E. 
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IROUT  If  IROUT  =  1  the  data  will  be  interpolated. 

If  IROUT  =  2  the  data  will  be  linearly  fitted. 
Default  value  for  IROUT  ■  1.  A  further  explana¬ 
tion  of  the  fitting  and  interpolation  routines 
will  be  given  within  the  program. 

GAP  1  the  gap  value  for  the  metal  whose  phonon  spec¬ 
trum  is  desired. 

GAP  2  the  gap  value  for  the  counter-electrode. 

IPLT  If  IPLT  *  1  the  original  graphical  input  will 

be  plotted  as  it  appears  on  the  laboratory 
X-Y  recorder.  If  IPLT  =  0  no  plotting  will 
occur.  Default  value  for  IPLT  =  0. 

IFIX  The  value  of  IFIX  will  determine  the  routines 

used  to  convert  the  input  values  to  dl/dV  values. 
In  the  present  program  only  one  routine  exists 
and  IFIX  =  1  as  the  default. 

ISS  If  ISS  =  0  the  first  derivative  data  is  not 

fitted  to  the  second  derivative  data  but  the 
NDS  and  o  are  found.  If  ISS  is  not  specified, 
t>econd  derivative  data  is  used, 

ISKIP  If  ISKIP  =  0  then  only  the  NDS  is  found.  Again, 
the  first  derivative  data  is  not  fitted  to  the 
second  derivative  data. 
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M 


INEGl 

INEG2 

INEG3 


IPLT1 

IPLT2 

IDEBG 


ISM 

ISMAX 


E3M 


M  is  the  number  of  points  around  the  point  of 
interest  used  in  both  the  interpolation  and 
fitting  routines. 

If  INEG  =  1  the  junction  was  reversed  biased 
when  taking  dl/dV  data.  If  INEG  =  0  the  junc 
tion  was  forward  biased.  INEGl  is  for  the 
superconducting  data;  INEG2,  the  normal  data; 
INEG3,  the  second  derivative  data.  Default: 
INEGl  -  1;  INEG2  and  INEG3  equal  by  default 
the  value  read  in  or  assigned  by  default  to 
INEGl . 

If  IPLT1  =  0  the  NDS  will  not  be  plotted. 
Default  value ,  IPLT1  *  1. 

If  IPLT2  =  .0  a  will  not  be  plotted.  Default, 
IPLT2  =1. 

If  IDEBG  =  1  the  input  data  points  with  anoma 
lous  points  removed  is  printed.  Default, 
IDEBG  =  0,  means  nothing  is  printed, 
the  number  of  points  greater  than  the  first 
relevant  point  at  which  smoothing  occurs, 
the  number  of  points  less  than  the  last  rele¬ 
vant  point  at  which  smoothing  occurs, 
a  smoothing  parameter  to  be  explained  within 
the  program.  Default,  ESM  =  .25H. 
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IFLAG  If  fitting  is  used:  If  IFLAG  =  1,  linear 
fitting  is  used;  If  IFLAG  =  2,  quadratic 
fitting  is  used. 

XSCLE1,  XSCLE2 ,  XSCLE3,  DXl,  DX2 ,  DX3,  EXl,  EX2 ,  EX3f 
DY1,  DY2 ,  DY3 ,  EYl ,  EY2 ,  EY3 ,  XI,  X2 ,  X3,  AX1,  AX2 ,  AX3 , 
YI,  Y2-,  Y3 ,  AY1,  AY2,  AY3 ,  ZEROl,  ZER02 ,  ZER03 ,  DXOl , 

DX02 ,  DX03 ,  DYOl,  DY02 ,  DY03,  BEROl ,  BER02 ,  BER03 ,  are 
parameters  for  subroutine  IFIX.  If  the  parameter  name 
ends  with  a  ^one'  it  is  associated  with  the  superconducting 
data.  If  its  name  ends  with  a  two,  the  normal  data;  if 
its  name  ends  with  a' three J  the  second  derivative  data. 

The  superconducting  parameters  have  no  default  values. 

Each  parameter  of  the  normal  and  second  derivative  data 
is  defaulted  to  be  equal  to  the  value  of  the  same  para¬ 
meter  for  the  superconducting  data.  For  example,  X2  =  XI 
and  X3  =  XI  are  the  default  values  for  X2  and  X3. 

XFS,YFS  the  experimental  superconducting  "first  deriva¬ 
tive"  data. 

XFN,YFN  the  experimental  normal  "first  derivative"  data. 
XSS,YSS  the  experimental  superconducting  "second 
derivative"  data. 
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FORMATS  FCR  INPUT  PARAMETERS 
INPUT  FORMAT 

TITLE  20A4 

The  elements  of  the  PAREM  NAMELIST  are:  H,  MAXX,  IROUT, 
GAP1 ,  GAP2,  IPLT,  IFIX,  ISS ,  M,  INEG1 ,  IPLTl,  IPLT2 , 
ISKIP,  Nl,  XSCLEl ,  YSCLEl ,  DX1 ,  EX1 ,  DYl,  EYl ,  XI,  AX1, 
Yl,  AY1,  ZEROl ,  DXOl ,  IDEBG ,  DYOl,  BEROl ,  I FLAG. 

The  elements  of  the  INPUT  NAMELIST  are:  ISM,  ISMAX, 

ESM,  DX2 ,  DX3 ,  EX2 ,  EX3 ,  DY2 ,  DY3 ,  EY2 ,  EY3 ,  X2 ,  X3, 
INEG2 ,  INEG3,  AX2 ,  AX3,  Y2 ,  Y3 ,  AY2 ,  AY3 ,  ZER02 ,  ZER03 , 
BER02,  BER03 ,  DX02 ,  DX03 ,  DY02 ,  DY03 ,  N2 ,  N3,  XSCLE2 , 
XSCLE3 ,  YSCLE2 ,  YSCLE3. 

INPUT  FORMAT 

XFS , YFS  10(1X,I1,  F6.0) 

XFN,YFN  10 (IX, II ,  F6.0) 

XSS ,YSS  10(1X,I1,  F6.0) 

OUTPUT  PARAMETERS 

TITLE,  NAMELIST  PAREM,  and  NAMELIST  INPUT  are  the  same 
as  their  definition  in  the  input  parameters  list. 

XFS , YFS  arrays  of  the . superconducting  dl/dV  x  and  y 

values,  respectively. 
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XFS ,YFS 


XFS,YFS 


XFS ,YFS 


If  IDEBG  =  1  then  the  previously  mentioned 
arrays  with  anomalous  points  removed  are 
printed . 

arrays  containing  the  x  and  y  coordinates , 
respectively,  of  the  input  points  as  they 
appeared  on  the  xy-recorder  expressed  in  terms 
of  inches  away  from  the  xy-recorders  midpoint, 
arrays  containing  the  voltage  and  corresponding 


«| 

I  , 


values . 


dV' SUPER 

The  exact  same  information  is  printed  for  dl/dV  N0RMAL 
data  and  the  d2I/dV2  data, 

FSI  the  normalized  first  derivative. 

Information  to  track  the  fitting  and  smoothing  section  is 
printed  out.  The  titles  printed  are  self-explanatory. 

ECS  the  values  of 


<L_ 

dV 


Cv-O 


Y-c£ 


c\E 


o\ 


FRAT 


a. 


146< 


63  - 


PUNCHED  OUTPUT 

TITLE  same  as  input 

IMAX  number  of  output  a.  FORMAT  -  (6X,I4) 

FRAT  a  FORMAT  -  7F10.6 

*  Notice  that  all  plotting  routines  listed  in  these 
programs  are  peculiar  to  the  MIT  Computation  Center. 
They  must  be  replaced  by  equivalent  programs. 

**  A  description  of  the  second  derivative  fitting 
routine  can  be  found  in  Ref.  11  and  Ref.  13. 
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SUBROUTINE  READR 


Subroutine  READR  reads  in  raw  input  data  and  puts  the 
x  values  sequentially  into  array  X,  and  the  corresponding 
y  values  sequentially  into  array  Y.  The  input  of  a  x  value 
is  of  the  form  Xa±nnnnn  and  its  numerical  values  is 
(±i.nnnn)10  a.  The  input  of  a  y  value  is  of  the  form 
Yainnnnn  and  its  numerical  value  is  computed  in  an 
analogous  manner.  As  an  example,  if  the  first  input  card 
is  as  follows: 

X4+01299Y4+0329 8X4+013 02Y4+03 32 0X4+01 3 08Y4+03332X4+01312Y4 
+03348X4+01318Y4+03356 


then  X(l)  =  .1299 

Y (1)  =  .3298 

X (2)  =  .1302 

Y ( 2 )  =  .3320 

X<3)  =  .1308 

Y (3)  =  .3332 

X(4)  =  .1312 

Y (4)  =  .3348 

X(5)  =  .1318 

Y (5)  =  .3356 

INPUT  PARAMETERS  (passed) 

N  =  the  number  of  x(y)  values  to  be  read. 


1G2<- 


-79- 


input  PARAMETERS  (from  I/O  device) 

X  =  the  input  x  values  in  untranslated  form. 

Y  -  the  input  y  values  in  untranslated  form. 

OUTPUT  PARAMETERS 

X  =  a::  array  of  the  translated  input  x  values 

Y  =  an  array  of  the  translated  input  y  values 

THE  INPUT  RECORDS  ARE  OF  THE  FOLLOWING  FORM 
5 (Xalnnnnn  ,Ya±nnnnn) . 
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SUBROUTINE  EVENO 


If  a  routine  requires  argument  and  function  arrays  as  input, 
but  the  given  input  data  has  certain  anomalies  where  the 
argument  values  are  not  monotouically  increasing,  EVENO 
when  called  will  drop  all  argument  and  their  corresponding 
function  values  from  the  input  arrays  that  cause  the  input 
argument  array  to  be  non-monotonically  increasing.  If 
several  arguments  which  follow  successively  are  equal, 
one  argument-function-value-pair  is  generated  whose  function 
value  equals  the  average  of  the  function  values;  e.g.,  if 
X  is  an  array  of  argument  values,  and  Y  is  an  array  of 
function  values  if 

X  =  1,1, 1,1, 2, 3, 4, 2, 3, 1,5, 5, 4, 6, 6 
and 

Y  =  0,1,2,3,4,5,6,7,8,9,10,11,12,13,14 

and  X  and  Y  are  input  to  EVENO,  they  will  be  output  as 
follows : 

X  =  1,2, 3, 4, 5, 6 

Y  =  3/2,4,5,6,10  1/2,  13  1/2 

(referring  to  the  original  arrays  the  pairs  with  y  =  0.1,2, 


lGo< 
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and  3  are  averaged,  y  =  7,8,9,  and  12  are  dropped,  y  =  i0  and 
11  are  averaged,  and  y  =  13  and  *4  are  averaged). 


INPUT  PARAMETERS 

X  -  an  array  of  argument  values. 

Y  an  array  of  corresponding  function  values. 

N  -  the  first  N  points  of  X  and  Y  are  argument  and  function 
values  respectively. 

OUTPUT  PARAMETERS 

X  -  an  array  of  monotonically  increasing  argument  values. 

Y  -  an  array  of  their  corresponding  function  values. 

ITC  -  the  first  ITC  points  of  X  and  Y  are  argument  and 

function  values  respectively. 
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IFIX* 

This  subroutine  was  written  for  a  particular 
experimental  investigation.  In  these  experiments 
it  became  convenient  (as  it  usually  is)  to  off-set 
the  absolute  magnitude  of  the  conductance  signal  so 
that  a  small  portion  of  the  curve  could  be  centered 
on  the  xy-recorder  and  measured  on  a  highly  magnified 
scale.  Hence,  the  purpose  of  the  subroutine  is  to 
convert  the  xy-recorder  data  to  its  absolute  values. 
Schematically  one  could  say,  we  seek  the  nor¬ 
malized  conductance  ratio: 

FS  Rq  +  (data  point) s 

W  ko  V  (dTata*  point  )n 

where  "(data  point)"  values  are  read  into  the  program 

and  R  is  a  constant  determined  from  the  dial  settings 
o 

of  the  xy-recorder  and  the  Doric  voltmeter.  In  a  crude 
way  we  may  also  represent  the  xy-recorder  traces  as: 


*  Only  the  general  features  of  this  subroutines  will 
be  considered. 


17  0< 
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Physic  ally .  the  signal  on  the  xy-recorder  as 
measured  from  ground  is  too  "flat"  to  observe  fine 
structure,  so  the  conductance  curve  is  shifted  down 
to  a  new  xy-recorder  zero  and  magnified  by  a  factor 
tf  (in  chis  case  10)  for  sensitivity.  With  respect 
to  the  new  axis  the  point  (z)  becomes  BERO  with 
an  amplitude  1C  times  the  original  signal  on 

the  xy-recorder.  Thus,  the  shift/point  from  the 
(0,0)  — -»  (0,0)/  recorder  zero  would  be: 

A=(ZER0-DY0)  -  BERO 

~vr 

This  gives  an  absolute  signal  magnitude  of* 


3=  10 ( ZERO- DYO ) - 3EP 0  1 

TO  T  15 

for  any  arbitrary  point  s(I),  measured  with  respect 
to  the  recorder  zero  (0,0). 


All  "values  for  this  experiment  are  multiplied  by 
the  factor  H.  Since  we  desire  a  ra^io  of  super 
conducting  to  normal  state  conduct;nces .  this 
causes  no  problem. 
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SUBROUTINE  LINFIT 

Subroutine  LINFIT  finds  the  values  of  a  function  F  at  an  equally 

spaced  interval  h.  The  inputs  are  (x,y)  pairs  such  that 

y=f  ix) .  f  ‘^ntains  random  noise  and  F  is  to  be  a  smoothed  f. 

To  determine  F(x0)  the  points  in  the  vicinity  of  x„  are  least 
square  fitted  to  either  a  straight  line,  P(x),  to  a  parabola, 
P(x),  and  F(x0)  =  P(x0). 


INPUT  PARAMETERS 

H  -  the  values  of  F(x)  are  generated  at  0,H,2H,... 

X  -  an  array  of  input  argument  values,  this  array  must  be 
mcnotonically  increasing. 

N  r  only  the  first  N  points  of  arrays  X  and  Y  are  valid 
argument  and  function  values  respectively. 

M  -  the  number  of  points  around  xg  which  are  used  to  determine 
P(x)  by  fitting.  If  M  is  even,  M/2  points  whose  arguments 
is  less  than  x0  and  M/2  points  whose  argument  is  greater 
than  xQ  are  used.  If  M  is  odd,  M^2  -  1/2  points  whose 
argument  is  less  than  x0  and  M/2  +  1/2  points  who'-e 
argument  is  greater  than  x0  are  used. 


4 


IFLAG 

MAXA 

i 

OOTPU , 

XT 

YT 

MAXO 


-  if  IFLAG  =  1,  linear  least  square  fitting  is  used. 

If  IFLAG  =  2,  quadratic  least  square  fitting  is  used. 

-  values  of  F(x)  are  generated  for  values  up  to 
(MAXA-l)H.  If  the  greatest  input  argument,  c,  is 
less  than  (MAXA-l)H,  the  values  of  F(x)  are  generated 
only  for  x  <  c. 

PARAMETERS 

-  an  array  of  output  arguments,  namely  0,H,2H,... 

-  an  array  of  output  function  values,  namely,  l'(O), 

F (H) , .  . . 

-  the  index  in  arrays  XT  and  YT  of  the  first  meaningful 
argument  and  c-~iction  value  respectively.  [F(x)  =  0 
for  all  x  which  are  less  than  the  minimum  input 
argument. ] 


MAXA  -  the  index  in  arrays  XT  and  YT  of  *-.he  last  meaningful 
argument  and  value  of  F(x)  respectively. 


TEST  IS  XINT  IS  >  MAXIMUM  INPUT  ARGUMENT 


!F  XINT  <  THE  MINIMUM  INPUT  ARGUMENT  THE  FULLt'K  ING  IS  OONE: 


AND  V  4RF  LOADED  WITH  THE  LAST  M  INPUT  POINTS 
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P  < 1 )  =  TY X2 
R(?)=TYXl 
R  (  3  )  =  TY 

CALL  THRFR  {  A,fl,  AA, BB,CC  ) 
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SUBROUTINE  INTER 

Subroutine  INTER  finds  the  values  of  a  function  F  at  an  equally 
spaced  interval  h.  The  inputs  are  (x,y)  pairs  such  that 
y  =  F  (x) ,  To  determine  F(xq)  the  n  points  in  the  vicinity  of 
x0  are  fitted  to  a  n-1  degree  polynomial,  P(x),  and 
F(x0)  =  P(xo), 

INPUT  PARAMETERS 

M,H,X,Y,N,  and  MAXA  are  the  same  as  in  LINFIT.  A  and  B  are 
dynamic  storage  arrays  which  are  used  in  subprograms  but  are 
-  not  declared  in  the  input  lists. 

OUTPUT  PARAMETERS 

XT,  YT,  MAXO,  and  MAXA  are  the  same  as  in  LINFIT. 
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SUBROUTINE  THREE 


Subroutine  THREE  solves  a  system  of  simultaneous  linear 
equations  by  computing  the  solution  in  terms  of  the  coefficients 
of  the  variables  and  the  constant  terms  in  the  equations. 


INPUT  PARAMETERS 

If  the  equations  -e  of  the  following  form 

anx  +  ai2y  +  ai3z  =  bi 
a21X  +  a22y  +  a23Z  85  b 2 

a,.x  +  a.  y  +  a  z=b. 

31  3  2-1  33  3 

A  is  the  array  of  the  input  coefficients  of  the  variables. 
B  is  an  array  such  that  B(l)  =  bj,  B(2)  ~  b2,  B(3)  =  b3. 

OUTPUT  PARAMETER 

X  -  the  value  of  x. 

Y  -  the  value  of  y. 

Z  -  the  value  of  z. 
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SUBROUTINE  FBCS* 


The  BCS  theory  predicts  that  the  tunneling 
current  for  two  weak-coupling  superconductors  is: 


V-  A 


Oi 


E  -  <L\f 


■>/e» - 


d  E  (O 


Therefore,  the  normalized  tunneling  conductance  is 


(  slL) 

'  <*v4 

cA~r  \ 

avL 


( 


A. 


N»-  A, 


01. 

E 


E  -tV 


cA^- 


ot 


Now  in  general, 
V>(^ 


dx  \  ?(x,)dx  -  Ux)4x  ViCx)  I  -  9(x)  •ak  ‘'U') 

aCx)  « 


X*fc(>0 


x=c»(x) 


r  *CO 

♦  )  fettxMx 


Cb) 


*  This  section  was  kindly  provided  by  J.  Bostock. 
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Because  of  the  singularities  in  the  two  densities  of 
states  in  Eqn. (2) ,  the  term  analogous  to  the  first 
two  terms  in  Eqn. (3)  will  be  non- analytic.  Hence, 
the  general  form  of  Eqn. (3)  applied  to  Eqn. (1)  is  not 
suitable  for  programming.  To  avoid  this  difficulty  a 
different  form  of  Eqn.(l)  must  be  used.  The  current 
Iss  can  also  be  written  as 


This  is  most  easily  seen  by  writing 

&OX 

and  evaluating  the  integrands  ^  and  Jl/  at  equivalent 
points  in  their  ranges  of  integration,  that  is  eval¬ 
uating  <S  atfcthe  lower  limit,  AQ1,  then  at  AQ1  +  H,..., 
(V  -  Aq2)  -  H,s(V  “  AQ2);  and  c?'  at*(V  -  A  )  , 

(V  -  A^)  +  H „ . . . ,  Aq2  “  h/,'a02*  eas*lY  seen 

then  that: 

XuivuJb  C  .  t 

e^UL  (,«•  (upper  limit)  =  <5.  (upper  limit)j 

d  (raid-point)  =  <!'  (mid-point) 

^  ^lower  =  rft *  (lower  limit)  j 

for  example;  and  this  is  true  for  all  "H-intervals" . 
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Thus, 


/-^O* 

0«  * ) 


where  VI  is  the  mid-point  of  the  first  integral: 


(O 


and  V2  is  the  mid-point  of  the  second  integral: 

va=(v-£>01V^(lW^-kA>- <  iW-t  -  ^  . 

Thus,  the  analogous  aquation  to  Eqn. (3)  is  the 


derivative  of  the  Eqn. (6): 

♦  +  S  ^(Is^OdE  (T> 

The  final  form  of  Eqn. (7)  then  becomes  (realizing 


that 


jpHEISi 


Tfa 


.  a 
koa. 


J07. 


U£-V>*-  aTI 


where  VI  =  1/2  (V-  AQ2  +  AQ1) 
and  V2  =  1/2  (V  -  AQ1  +  AQ2)  . 
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This  functior  was  not  that  used  in  the  original 
inversion  program.  The  accuracy  with  which  it  eval¬ 
uates  the  normalized  conductance  can,  however,  be 
seen  by  using  it  to  calculate  the  BCS  density  of  states 
for  a  normal-superconductor  tunnel  junction.  That  is, 
the  calculated  values  from  this  function  should  be  com¬ 
pared,  in  this  case,  to  (Aq2  =  0): 

E 

The  comparison  shows  that  FBCS  is  correct  to  5  signifi¬ 
cant  figures. 
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STA^E  II 


Eqn. (11)  gives  the  expression  for  the  normalized  first 
derivative  of  a  tunnel  junction  assuming  that  both  metals 
of  the  junction  have  BCS  NDS's.  If  the  metals  don't  have 
BCS  NDS's  the  analysis  of  Chapter  I  in  deriving  Eqn. (11) 
still  applies  and 


where  Pj  is  the  NDS  of  one  electrode,  D2  the  NDS  of  the 
other,  and  C  the  experimental  normalized  first  derivative. 
Aqi  and  Aq2  are  the  BCS  energy  gap  values  of  the  respec¬ 
tive  metals.  D2,  Aoi*  Aq2,  and  C  are  assumed  to  be  known 
so  that  STAGE  II  solves  for  D1.  * 

INPUT  PARAMETERS 
TITLE  same  as  in  STAGE  I 

GAPl  A01 

GAP 2  A02 

the  interval  at  which  values  of  C  and  D2  are 
known  and  at  whicv  the  values  of  Dy  are  desired. 


#~VO.  | — 


H 
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The  following  6  integer  variables  are  indexing 
parameters. 

IEXPF  (IEXPF)H  +  GAPl  +  GAP 2  is  the  first  value  of 
V  for  which  C  is  known.  For  V  <  (IEXPF) H  + 
GAP1  +  GAP2  C  is  assumed  to  be  equal  to 


IEj&Jj  (IEXPL)H  +  GAP1  GAP2  is  the  last  value  of  V 
for  which  C  is  known. 

IGIVF  (IGIVF-1)  H  +  GAP 2  is  the  first  value  of  o>  for 
which  D2  is  known.  For  E  <  (IEXPF) H  +  GAP 2 
D2  is  assumed  to  be  equal  to 


E 


IGIVL  (IGIVL-l)H  +  GAP2  is  the  last  value  for  which 
D2  is  known. 

LF  (LF)H  +  GAP1  is  the  first  value  of  V  at  which 

Di  is  desired.  For  £  <  (LF)H  +  GAP1 

■  E  _ 

*  o> 
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LL  (LL)  H  +  GAP1  is  the  last  value  of  V  i.t  which 

D2  is  desired. 

EXP  C  in  a  form;  namely,  divided  by  the  BCS  pre¬ 

dicted  values  and  decreased  by  one. 

GIVEN  in  a  form. 

INPUT  FORMATS 

INPUT  FORMAT 

TITLE  20A4 

The  elements  of  the  INPU1  NAMELIST  are  LF,LL,IEXPF,  IEXPL, 
IGIVF,  and  IGIVL. 

GAPl,  ,H  3(10X,F10.5) 

GIVEN  (I).  1=1,  IGIVL  7F10.6 

EXP (I) ,  1=1,  IEXPL  7F10.6 

OUTPUT  PARAMETERS 

TITLE,  NAMELIST  INPUT,  GAPl,  GAP2 ,  and  H  are  the  same  as 
their  definition  in  the  input  lis^. 

EXP 


the  input  from  STAGE  I,  the  first  time  EXP 
is  printed  out. 


n 


GIVEN 


UNKWN 


ERR 

EXP 

GIVEN 

UNKWN 

UNKWN 

UNKWN 


-125- 

E 


the  BCS  density  of  states,  the  first  time 
GIVEN  is  printed  out. 


the  first  time  UNKWN  is  printed  out. 

is  the  change  in  value  for  the  equivalent 

BCS  density  of  states  in  the  interval  H. 

the  convoluted  input  in  a  form. 

the  values  c f  D2  in  a  form. 

the  values  of  which  are  assumed  equal  to 

the  BCS  predicted  values 

the  values  of  D^. 

the  values  of  in  a  form. 


2C8< 
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DINT (GAP, H) 


SQ(X,A)  = 


END  = 


ENC  = 


FUNCTIONS 


r  GAP  +  H 

_ E 


J. 


■J 


dt 


E  -GAP 


\ 


J  '\it1-&/\P2 

H 


dt 


GAPl  +  H 
Ea 


dt 


2G9< 


n-E  FIRST  TIME  THROUGH  THIS  LOOP  BCS  VALUE S  ARE  PRINTED  OUT  FOR  BOTH  METALS 
AS  WFLL  AS  THFIP  BCS  CONVOLUTION. 

THE  SEC JNO  TIME  THROUGH,  THE  ACTUAL  NORMALIZED 

CCNDUCT ANCF  FOR  THF  JUNCTION  IS  PRINTED  OUT,  AS  WELL  AS  THE  DENSITY  OF 
STATES  FUR  THE  KNOWN  (CR  GIVEN)  MFTAL. 
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STAGE  III* 

The  object  of  this  stage  is  to  find  a  complex  energy 
gap  function  A(E),  a  complex  renormalization  function  Z(E), 
an  effective  phonon  spectrum  a2(E)F(E)  ,  afid  a  Coulomb  pseudo¬ 
potential  y*  which  will  reproduce  the  tunneling  electron  den¬ 
sity  of  states  and  also  gives  the  measured  value  for  the  gap 
function  at  the  gap  edge,  Aq  =  A(e=Aq).  The  theory  behind 
these  calculations  has  been  discussed  in  Chapters  2  and  3. 

FUNCTION  NAMES  AND  INPUT  PARAMETERS 
Function  Name 

SAMPLE  Identifies  the  input  data 
DM  Initially  is  used  to  label  the  input  data 

TZPE  Numerically  computed  derivative  of  the  input  data 
TZ  Calculated  density  of  states,  N  (E) 

V# 

El  Energy  in  10” 3  electron  volts 

G  (a)  a2F (Eq) 

(b)  The  above  divided  by  a  sort  of  average  value 
of  a2;  equal  to  la2F(Eq)  r/a2F(Eq)dEq] . 

DG  The  incremental  correction  to  o2F (Eq) . 

DR  A^ (E) ,  the  real  part  of  the  gap  function 

DI  Aj (E ) ,  The  imaginary  part  of  the  gap  function 

ZP.  Z^(E)  ,  the  real  part  of  the  renormalization 

function 

*  The  ~ l lowing  has  been  adapted  from  Ref.  (13)  . 


-TTTrg^T  tat  , 


ZI 


lb<t 

ZjtE),  the  imaginary  part  of  the  renormal¬ 
ization  function 
FR  Real  { [1  -  Z(E)]E} 

FI  Imag  {U  -  Z(E)]E} 

PAR  Real  (A(E)Z(Z)  } 

PAI-.  Imag  {A(E)Z(E)> 

Input  Parameters 

JO  The  index  specifying  the  maximum  phonon  energy, 

i 

Eq.  For  an  energy  spacing  of  H=.lxlO-^  electron 
volts,  for  example,  J0=161  corresponds  to 

E  =16.0x10""^  electron  volts. 

o 

JC  The  index  specifying  E  .  JC=451  corresponds  to 

C 

E  *45, 0x10""-*  electron  volts  for  H=.lxl0“^  electron 
c 

volts,  for  example. 

MAXG  The  index  specifying  the  maximum  energy  to  which 

2 

corrections  to  a  F  will  be  made  in  the  iteration 
procedure.  Usually  MAXG  is  set  equal  to  JO,,  but 
lower  values  may  be  useful  in  refinements  of  a 
previously  calculated  spectrum. 

NITER  Number  of  iterations  of  the  gap  equation. 

Number  of  iterations  in  which  a^F(Eq)  is  refined. 

I 

I 

ZZ2< 


NREF 


1**0. 

-• 

NG  This  is  the  ratio  of  the  energy  spacing  of  the 

input  a^F  (Eq)  to  the  energy  spacing  given  by  the 
parameter  H.  For  example,  if  the  input  <j^F(Eq) 
has  an  energy  interval  of  0*l3tlO“^eV  while  H=.05, 
then  NG=2.  The  program  will  interpolate  to  obtain 
the  same  energy  interval  for  all  the  functions. 
However,  if  NG  <_  0,  no  data  will  be  read  in  and 
the  program  will  use  a  starting  a2F(Eq)  of 
tt2F(Eq)  =  (Eq/H)  2  (E^-Eq/H)  2 . 

ND  This  is  the  ratio  of  the  energy  spacing  of  the 

input  A^(E)  and  A2(E)  to  the  energy  spacing  given 
by  the  parameter  H,  in  the  same  manner  as  NG 
served  for  the  input  a2F(Eq)  data.  For  ND  _>  1, 
the  program  will  interpolate  to  obtain  the  same 
energy  interval  as  for  the  other  functions .  For 
ND  <  0,  no  data  will  be  read  in  and  the  program 
will  use  a  starting  function  of  A^=Aq  for  E  <  Ec, 
A-j=0  for  E  >  Ec,  and  Aj  =  0. 

NH  This  is  the  ratio  of  the  energy  spacing  given  by 

the  parameter  H  to  the  energy  spacing  of  the  in¬ 
put  density  of  states  data,  DM.  If  NH=2,  for  ex¬ 
ample  ,  then  the  program  will  use  only  every  other 
data  point. 


Z23< 


1*41 


MING  The  index  specifying  the  energy  up  to  which 

a^F(Eq)  is  to  be  taken  as  proportional  to  Eq2. 
MING=9,  NH=1 ,  and  H=.l  mean  that  a2F(Eq) 
is  taken  to  be  proportional  to  Eq2  up  to 
I (MING/UH)  +1]  x  H,  or  lxlO"3  eV. 

N  The  starting  value  for  the  Coulomb  pseudo¬ 

potential  y*.  Usually  set  equal  to  0.11. 

DO  This  is  the.  GAP1=A^  of  the  two  earlier  stages. 

H  The  energy  spacing,  in  eV  x  10"3,  for  the  cal¬ 

culations  . 

NR  The  number  of  data  points  in  the  input  density 

of  states.  If  the  second  stage  is  used,  these 
last  three  parameters  will  be  part  of  the  punched 
output  of  that  stage. 


INPUT  DATA  FORMAT 
Data  Card  Number 

1 

2 

3 

4 


Information  on  Card  Format 

JO, JC,MAXG, NITER, NREF,  7(6X,I4) 

NG, ND 

NH, MING,W  2(7X,I3) ,10X, 

F10.2 

SAMPLE  20A4 

DO,H,NR  (6X,F10 . 6 , 3X, 

F10.  6 ,4X,I5) 


22%< 


-^6- 

5ff 

DM 

7P10.6 

6ff (only  if  NG>1) 

G  ‘ 

1P7E11. 3 

7f f  (only  if  ND>1) 

DR 

7E11.3 

8ff(only  if  ND>1) 

DI 

7E11.3 

nor*-- 


IF  PROGRAM  IS  BEING  USED  TO  CALCULATE  IOI/DVI  GIVFNt  FOR  STAGE 
INSEPT  INSTRUCTION  N ( EXP  )  =  N ( CALCULATED ) . 
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Chemical  Synthesis  Using  High  Temperature  Lithium  Vapor  Species 


Abstract  and  Objectives 


Recent  work  in  our  laboratory  has  resulted  in  the  discovery 
of  several  new  routes  to  polylithioc?  rbons,  perlithiocarbons,  and 
inorganic  polylithium  species.  Previously  for  these  classes  of  com¬ 
pounds,  there  were  very  few  known  examples,  and  there  were  no 
general  synthetic  routes  in  the  literature.  The  new  synthetic  methods 
are  all  in  the  early  stages  of  development  and  all  involve  the  reaction 
of  high  temperature  lithium  vapor  in  the  range  of  800  -  1000°C  with 
various  organic,  inorganic  and  polymeric  species.  During  the  past 
year  we  ha*re  succeded  in  the  synthesis  of  the  first  perlithioalkanes, 
tetralithiomethane^  C(Li)4  and  hexalithioethane,  C2Li6  .  We  have  also 
prepared  hexalithiobenzene  C6Li6  .  An  intensive  study  of  these  species 
is  now  underway  and  efforts  are  underway  to  prepare  other  polylithiated 
species.  Support  for  this  work  is  urgently  needed  to  expand  the  small 
effort  in  the  exciting  new  area.  It  appears  that  polylithio:arbons  will 
be  very  useful  reagents  in  organic  synthesis  as  well  as  important  mon¬ 
omers  and  crosslinking  agents  in  polymer  synthesis.  The  synthesis 
of  jeveral  type3  of  th.re'*  dimensional  polymers  and  igh  temperature 
materials  is  proposed  based  on  routes  involving  these  new  polylithium 
species.  The  preparation  of  a  series  of  diamond-like  polymeric  mater¬ 
ials  is  proposed  and  would  result  in  a  new  class  of  high  performance 
materials.  Polylithium  compounds  may  also  be  useful  as  catalysts  and 
a  new  synthesis  is  proposed  for  diamond  based  on  these  compounds. 
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Summary  of  Research  Progress  in  the  Synthesis  Characterization 
and  Reactions  of  Polylithium  Compounds 
6/1/73  -  U/ 1/73 

Initially  a  study  of  the  reaction  of  carbon  tetrachloride,  CC14 
and  the  reaction  of  hexachloroethane  CjClt  was  studied  with  lithium 
vapor  at  a  1000°C.  *  The  products  of  the  initial  reaction  were  tetra- 
Uthiomethane  CLi4  and  tetralithioethylene  Li2C-CLi2.  These  compounds 
were  characterized  (see  the  attached  reprint)  and  it  has  been  found 
that  they  react  with  various  organic  and  inorganic  substrates  in  the 
manner  characteristic  of  other  lithium  compounds  such  as  the  con¬ 
ventionally  prepared  monolithium  compounds.  This  experimental 
result  provides  evidence  of  the  suitability  of  such  polylithium  com¬ 
pounds  for  reactions  as  monomers  in  the  proposed  diamond-like 
polymer  systems.  Subsequently,  the  reaction  of  hexachloroethane 
with  lithium  produced  hexolithioethane  in  approximately  80*7«  yield.  • 
Again  this  material  has  been  shown  to  be  a  potentially  reactive  monomer. 

A  study  of  the  reaction  of  carbon  vapor  generated  frovn  a  carbon 

arc  apparatus  at  2500°C  was  undertaken  to  establish  an  alternative 

2 

method,  for  synthesizing  polylithium  compounds.  The  principal  product 
was  CjLi4  which  appears  to  have  the  allene  structure  and  might  be  an 
.  attractive  monomer.  Varying  amounts  of  tetralithiomethane  and 


tetralithioethylene  were  also  produced  in  this  study. 

A  third  publication  ftas  resulted  from  the  reaction  of  lithium 
vapor  with  benzene.  ^  While  for  most  of  our  systems  the  reactants 
are  admitted  to  the  reactor  at  room  temperature  and  reacted  with  the 
high  temperature  lithium  vapor,  initial  studies  giving  a  small  yield  of 
hexalithio  benzene  prompted  us  to  do  a  more  extensive  study  based  on 
the  reaction  of  lithium  vapor  plus  benzene  at  various  temperatures. 

A  new  type  of  reactor  was  designed  for  this  purpose.  At  room  temper¬ 
ature  one  gets  a  distribution  of  lithiated  benzene  ranging  from  over  the 
range  of  C^LU_nHn  where  n  =  1  to  6  .  At  room  temperature  the  species 
of  highest  yield  is  dilithiobenzene.  As  one  raises  the  temperature 
progressively  to  800° C  a  shift  toward  the  higher  benzene  polylithium 
compounds  is  noted.  The  temperature  of  the  benzene  is  raised  to  900°C 
and  a  striking  new  result  is  attained.  One  begins  to  obtain  polyli  thiated 
cyclohexanes  from  the  benzene  starting  material.  This  results  from  a 
tendency  for  the  lithium  to  add  to  the  double  bonds  as  well  as  replace 
the  protons  on  the  ring.  Studies  were  made  up  to  temperatures  of 
J.200°C  producing  such  striking  compounds  as  perlithiocyclohexane  and 
CiLUH*.  It  has  been  found  subsequently  that  the  reaction  of  hexachloro- 
.  benzene  with  lithium  provides  a  much  cleaner  route  to  h„xalithiobenzene. 

The  reaction  of  methane  gas  with  lithium  has  also  been  extensive];, 
investigated.  ^  We  predicted  and  found  that  lithium  at  1000°  will  not 
react  with  methane  gas  at  room  temperature.  The  activation  energy  for 


this  reaction  appears  to  be  in  the  30  kcal  mole  range.  Lithium  at 
1000°C  usually  has  about  3  kcal  mole  of  translational  or  kinetic  energy 
and  only  a  very  small  percentage  of  lithium  in  the  first  excited  elec¬ 
tronic  state.  Therefore,  the  activation  energy  for  the  lithium- methane 
reaction  is  not  likely  to  be  r.vailable  under  these  reaction  conditions. 
The  reaction  between  methane  and  lithium  would  be  a  very  important 
one  for  the  production  of  tetralithiomethane  which  we  consider  to  be 
one  of  the  important  materials  in  subsequent  polymer  synthesis.  We 
have  used  several  approaches  to  provide  the  activation  energy  for 
this  reaction.  We  nave  found  that  if  we  either  electronically  excite 
the  lithium  or  vibrationally  excite  the  methane  the  reaction  proceeds 
in  a  satisfactory  manner.  We  have  conducted  an  experiment  in  which 
the  methane  was  heated  to  900°C  and  obtained  satisfactory  results. 
Increased  temperatures  populate  higher  vibrational  strates  in  the 
methane  and  thus  allow  a  very  significant  decrease  in  the  activation 
energy  of  the  reaction.  We  also  used  a  radio  frequency  source  to  excite 
the  lithium  vapor  electronically  and  produced  a  successful  reaction. 

Still  another  piece  of  work  has  been  done  using  a  methane  plasma 
source  which  vibrationally  excites  the  methane  and  allows  the  reaction 
to  proceed. 

In  the  study  of  tetralithiomethane,  a  fourth  synthesis  is  being 
completed  which  involves  intitial  synthesis  of  dilithiomethane  and 


subsequent  pyrolysis  to  give  a  substantial  yield  of  tetralithiomethane. 
Using  this  method  we  have  obtained  substantial  quantities  of  tetra¬ 
lithiomethane.  Extensive  work  has  been  done  on  solvent  systems  for 
tetralithiomethane  so  that  one  may  do  the  polymer  synthesis  in  solution. 
Several  successful  solvent  systems  are  being  developed  and  the  concen¬ 
trations  of  the  polylithium  compounds  in  solution  are  being  determined 
by  NMR  studies  and  the  preparation  of  derivatives. 

An  extensive  study  has  been  recently  completed  on  the  reactions 
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of  olefins  with  lithium  vapor  A  large  number  of  new  polylithium  com¬ 
pounds  such  as  tetralithio  cyclohexadiene,  2,3  dilithio  2  butene, 

2, 3  tetralithio  2  butene,  1, 2  dilithio  isobutane,  1,4  dilithio  2  butene, 
and  l,  2  trilithio  isobutane  have  been  prepared  in  this  study.  We  have 
seen  instances  of  both  substitution  for  hydrogen  on  sites  requiring 
alow  activation  energy  for  reaction  and  addition  to  double  bonds. 

This  study  has  established  another  general  route  to  polylithium  com¬ 
pounds. 

Studies  of  reactions  of  inorganic  species  with  lithium  have  also 
been  undertaken.  The  reaction  of  silicon  tetrachloride,  SiCl4,  with 
lithium  has  been  recently  studied  and  the  polylithium  compound  SiL^ 
was  obtained  in  good  yield.  Hydrolysis  and  deuterolysis  of  the  com¬ 
pound  produced  silanes  and  perdeutero  silanes.  A  subsequent  study 
of  the  reaction  of  tetralithiosilane  with  methyl  chloride  resulted  in  a 
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15%  yield  of  tetramethyl  silane  indicating  that  the  SiLi4  is  also  a 
potential  monomer.  The  reaction  of  BC13  and  PC13  are  currently 
under  study.  A  fair  yield  of  trimethylborane  B(CH3)3  was  obtained 
from  the  reaction  of  BLi3  with  methyl  chloride.  Additionally  the  syn¬ 
thesis  of  polylithiated  carboranes  has  been  attempted  and  a  workup 
is  now  in  progress.  Isolation  of  a  higher  purity  SiLi4  is  being  under¬ 
taken.  It  should  be  pointed  out  that  if  the  diamond  syntheses  from 
tetralithiomethane  are  successful,  a  synthesis  of  a  silicon  analog  for 
a  diamond  will  be  attempted  using  the  SiLi4.  Such  a  diamond-like 
silicon  structure  would  have  many  unusual  electrical  and  physical 
properties.  This  type  of  structure  has  been  the  object  of  numerous 
unsuccessful  synthetic  studies  utilizing  the  phase  diagram  of  silicon. 
This  potentially  valuable  species  has  previously  eluded  all  other  syn¬ 
thetic  approaches. 

Recently  structural  studies  have  been  initiated  on  several  poly- 
lithium  compounds.  We  have  been  interested  in  the  lithium  and  carbon 
13  NMR  spectra  and  the  laser- raman  spectra  of  polylithium  compounds 
both  in  solution  and  in  the  solid  state.  Very  encouraging  results  have 
been  obtained.  Such  studies  will  also  lead  to  accurate  methods  for 
determination  of  concentration  of  reactive  polylithium  intermediates 
In  solution  and  for  identification  of  polylithium  compounds. 


The  synthesis  of  diamonds  which  were  proposed  from  the  poly- 
llthlum  compounds  Is  almost  operational  and  a  first  attempt  should  be 
made  within  a  month.  A  series  of  steel  containers  have  been  prepared 
for  handling  the  polylithium  compounds.  These  containers  are  cylin¬ 
drical  and  have  a  length  of  about  3  inches  and  are  3  inches  in  diameter. 
Arrangements  have  been  made  to  rent  an  ASEA  high  temperature  - 
high  pressure  press  from  a  local  industrial  firm.  This  press  has 
reaction  chamber  dimensions  of  6  inches  in  diameter  by  20  inches  in 
length,  and  is  capable  of  pressures  on  the  order  of  35  kilobars  at 
temperatures  above  1200°C.  The  size  of  this  high  temperature  reactor 
is  indeed  remarkable  and  results  font  a  technological  breakthrough 
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within  the  last  several  years.  Previously,  high  pressure  chemistry 
had  been  limited  to  rather  small  volumes.  These  temperatures  and 
pressures  are  not  adequate  for  the  synthesis  of  diamond  by  brute  force 
'methods  involving  graphite.  We  are  very  hopeful  and  thoroughly 
confident  that  our  methods  will  generate  diamonds  at  a  much  lower 
temperature  and  pressure.  At  this  stage,  of  course,  it  would  be 
impossible  to  predict  the  size  of  the  diamonds.  However,  :he  reactors 
have  been  designed  in  such  a  manner  that  a  successful  experiment 
would  result  in  many  karats  of  diamonds.  The  fee  for  this  press  is 
about  $650  per  run  and  we  are  able  to  run  approximately  10  reactions 
at  once.  This  is  done  by  placing  a  number  of  smaller  steel  containers 
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Inside  the  large  reaction  container .  Our  laboratory  has  worked  very 
hard  in  obtaining  the  samples  for  this  work  and  we  are  awaiting  results 
with  great  anticipation.  Details  of  the  proposed  diamond  syntheses 
are  contained  in  the  body  of  the  original  proposal. 
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ABSTRACT 


INFLUENCE  OF  STRAIN- INDUCED  MARTENSITIC  TRANSFORMATIONS  ON 
FATIGUE  CRACK  GROWTH  RATES  IN  STAINLESS  STEELS 

A.G.  rineau  and  R.M.  Pelloux 

The  fatigue  crack  growth  rates  (FCGR)  of  two  unstable  austenitic  stain¬ 
less  steels  (Fe-1 6Cr-13Ni )  and  (Fe-18Cr-6.5Ni-0.19C)  were  determined  in  the 

M  -M.  temperature  range  where  a  strain  induced  y  -*■  ct'  martensitic  transfor- 
s  d 

mation  occurs  near  the  crack  tip.  These  FCGR  were  compared  to  the  rates 
measured  in  the  stable  austenitic  phase  of  a  Fe-31 . 5Ni  and  a  Fe-34Ni  alloy 
$nd  in  the  martensitic  phase  obtained  by  quenching  the  Fe-31 . 5Ni  alloy  below 
M$.  In  the  Fe-31. 5Ni  alloy,  the  FCGR  are  an  order  of  magnitude  higher  in 
the  martensitic  than  in  the  austenitic  structures  for  A  K  <  40  ksi  /uT  .  The 
FCGR  of  the  stainless  steels  decrease  markedly  when  the  test  temperature  ' 
approaches  in  the  Ms-Md  range.  The  FCGR  for  the  alloy  Fe-18Cr-6.5Ni-0.19C 
In  a  warm-worked  condition  are  consistently  higher  than  for  the  same  alloy  in 
the  annealed  condition  for  AK  <  40  ksi  /Tn.  The  results  are  discussed  in 
terms  of  the  influence  of  phase  structures,  stacking  fault  energy  and 
work  hardening  exponent  on  the  FCGR, 


INFLUENCE  OF  STRAIN-INDUCED  MARTENSITIC  TRANSFORMATIONS  ON 
FATIGUE  CRACK  GROWTH  RATES  IN  STAINLESS  STEELS 

A.G.  Plneau  and  R.M.  Pelloux* 

INTRODUCTION 

Martensitic  transformations  can  be  induced  by  plastic  deformation  in 
the  temperature  range  Mg  to  M^.  M$  is  the  temperature  for  spontaneous 
transformation  and  Md  is  the  maximum  temperature  at  which  the  transformation 
can  be  induced  in  a  tensile  test.  A  comparative  study  of  these  phase  trans¬ 
formations  in  Fe-Ni-Cr-C  alloys  was  recently  completed  bv  Lecroisey  and 
Pineau  (1).  Depending  upon  the  alloy  composition,  the  Mr-Md  temperature 
range  can  vary  from  30  to  200°C.  The  work  hardening  exponent  and  the  UTS  of 
an  alloy  undergoing  strain- induced  phase  transformation  increase  markedly 
as  the  defonnation  temperature  approaches  M  .  The  uniform  elongation,  which 

is  also  strongly  dependent  upon  the  test  temperature, shows  a  well  defined 
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maximum  in  the  Ms-Md  temperature  range. 

Following  the  results  of  Lecroisey  and  Pineau,  it  was  decided  to  use 
some  of  those  alloys  for  which  the  M  and  Md  temperatures  and  the  tensile 
properties  were  well  known,  in  order  to  study  tne  influence  of  the  strain 
induced  martensitic  transforation  on  the  fatigue  crack  growth  rate  (FCGR). 
The  aim  of  this  investigation  was  to  measure  and  compare  tie  FCGR  of  the 
alloys  in  four  different  conditions: 

a)  Stable  annealed  austenite 

b)  Warm-worked  austenite 

c)  Martensite  formed  by  quenching  belov/M,. 

*  A.G.  Pineau  is  Maitre  de  Recherches  at  Centre  des  Materiaux  de  l'Ecola 
des  Mines  de  Paris  -  91  Corbeil  Essonnes,  France.  R.M.  Pelloux  is  Asso¬ 
ciate  Professor,  Massachusetts  Institute  of  Technology,  Cambridge,  Mass. 
02139 
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d}  Strain-induced  marten. ite  at  the  fatigue  crack  tip  in 
an  austenitic  matrix. 

The  alloys  tested  in  condition  (a)  were  Fe-Ni  compositions  and  two 
stainless  steels,  a  1 6-Cr- I 3Ni  and  a  type  301  steel  with  a  high  carbon  content. 
This  last  alloy  was  also  tested  in  condition  (b),  because  it  is  representa¬ 
tive  of  a  new  class  of  high  strength,  high  toughness  steels  referred  to  as 
TRIP  steels.  The  comparison  between  'J..2  FCGR  in  austenite  and  martensite 

(conditions  a  and  c)  was  made  on  a  Fe-3 I . 5Ni  alloy.  The  influence  of  the 

(condition  d) 

strain  induced  martensitic  transformation  on  the  FCGRVwas  determined  for  the 
two  stainless  steels  which  were  tested  at  various  temperatures  in  the  M  -M . 

t,  5  0 

temperature  range. 

Most  of  the  fatigue  crack  growth  rate  tests  were  performed  in  dry  argon 
to  minimize  the  effect  of  humidity.  Microhardness  measurements  were  used 
to  determine  the  plastic  zone  sizes  and  to  assess  the  degree  of  cyclic  harden- 
Ing  or  softening  near  the  crack  tip.  The  microscopic  crack  growth  rates  were 
measured  by  electron  microscopy  and  compared  to  the  overall  growth  rates. 

The  test  results  serve  as  a  basis  for  the  evaluation  of  the  influence 
of  austenite  stacking  fault  energy,  work  hardening  exponent  and  phase 
structure  on  fatigue  crack  growth  rates. 

EXPERIMENTAL  PROCEDURE 

1.  Materials  -  Composition  and  Preparation 

The  chemical  analysis  of  all  the  alloys  are  given  in  Table  I.  The  Fe-Ni 
and  16-13  alloys  were  induction  melted  under  vacuum  and  cast  as  20  lb.  ingots. 
These  ingots  were  cut  length-wise  into  two  pieces,  reduced  50%  by  hoc  rolling 
and  homogenized  for  70  hours  at  1200°C  under  argon.  They  were  firally  hot 
rolled  to  plates  of  9  mm  in  thickness. 

Tensile  specimens  cut  parallel  to  the  rolling  direction  were  sealed  in 
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evacuated  vycor  capsules  and  given  the  austenitizing  treatment  reported  in 
Table  I,  followed  by  water-quenching.  The  austenitizing  conditions  were 
selected  to  give  a  final  grain  size  of  approximately  100  microns. 

A  200  lb.  heat  of  a  301  type  stainless  steel  was  used  to  compare  the 
tensile  and  fatigue  properties  of  anneal edand  warm-worked  austenite.  The 
annealed  301  specimens  were  heat  treated  and  prepared  ir.  the  same  manner  as 
the  Fe-Ni  and  16-13  alloys.  The  warm-worked  specimens  were  obtained  by 
annealing  25  mm  thick  plates  at  1120°C  followed  by  water  quenching  and  by 
rolling  65%  at  450-500°C.  The  reduction  of  thickness  from  25  to  9  mm  was 
achieved  in  about  10  passes.  Betv/een  eaqh  pass,  the  plates  were  reheated  to 
500°C  in  a  salt  bath  furnace.  The  total  warm-working  operation  took  approxi¬ 
mately  30  minutes.  Tensile  specimens  parallel  and  perpendicular  to  the  rol¬ 
ling  direction  were  cut  from  these  warm-worked  plates. 

Double  cantilever  beam  fatigue  specimens  were  cut  from  9  run  thick  plates, 
rough  machined,  annealed  at  the  temperatures  given  in  Table  I  and  water  - 
quenched.  Cutting  of  the  notch  and  final  machining  was  always  done  after  the 
final  quench  in  order  to  minimize  residual  stresses.  In  the  fatigue  tests, 
the  notch  was  perpendicular  to  the  rolling  direction.  The  phase  transforma¬ 
tion  temperatures  (M$  and  M^)  were  determined  by  the  same  methods  which  have 
been  used  previously  (1,2). 

2.  Mechanical  Properties 
2:1  Tensile  tests 

Cylindrical  tensile  specimens  of  3  mm  diameter  and  30  mm  gage  length 

-3  -1 

were  tested  at  a  deformation  rate  of  2.8  x  10  sec  *  The  specimens  were 
immersed  in  liquid  baths  (liquid  nitrogen,  alcohol  +  dry  ice,  oil)  controlled 
to  +  2°C  in  the  range  -196°Cto  150°C.  In  the  two  stainless  steels  and  the 
Fe-31.5  alloy,  the  intensity  of  magnetization  was  measured  continuously  during 
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tenslle  tests.  The  magnetic  measurements  were  carried  out  by  using  a  double¬ 
coil  device  described  elsewhere  (2).  These  measurements  give  a  good  estimate 
of  the  volume  fraction  (f)  of  strain  induced  a'  martensite. 

2:2  Fatigue  tests 

Double  cantilever  beam  specimens  were  used  for  the  fatigue  crack  propa- 
*  gation  tests.  The  specimens  dimensions  were  70  mm  wide  x  67.2  mm  high 
x  7  ran  thick.  The  original  notch  length  was  32.5  mm.  The  tests  were  run  at 
20  Hz  with  sinusoidal  tension  stress  waveform  (R  =  0.05)  in  a  load  controlled 
closed  loop  hydraulic  fatigue  machine.  Crack  lengths  were  measured  to  an 
absolute  accuracy  of  0.02  mm.  The  stress  intensity  range  A  K  was  measured 
from  the  following  equation  (3): 

ap  ,  1/2  .  3/2  a  5/2  7/2  9/2 

4  K  «  [29.6(i)  -  185.5(A)  .  655.7(f)  -1017(f)  .  638.9(f)  ] 

with  AP  the  applied  load  range;  a  =  crack  length;  w  =  56  mm;  B  =  7  mm. 

The  da/dN  measurements  were  taken  only  after  the  initial  fatigue  crack  length 
increment  had  reached  a  value  between  0.5  and  1  ran. 

Most  of  the  tests  were  run  in  an  inert,  argon  environmt.it  (less  than 
200  ppm  of  water  vapor).  For  tests  run  in  air  at  25°C,  the  humidity  was 
between  20  and  60*.  An  air-tight  temperature  chamber  was  built  for  tests  at 
various  temperatures  between  -30°C  and  100°C.  The  low  temperatures  were 
achieved  by  passing  cold  alcohol  in  a  copper  coil  within  the  test  chamber. 
Heating  tapes  were  used  for  the  experiments  above  25°C.  The  temperature  of 
the  specimen  was  measured  with  a  thermocouple  and  was  controlled  to  +  2°C. 

2:3  Metallography  and  microhardness  measurements 

Microhardness  measurements  at  mid-thickness  of  the  fatigue  specimens 
were  used  to  determine  the  plastic  zone  sizes  on  each  side  of  the  fatigue 
fracture  surfaces.  The  measurements  were  made  at  different  crack  lengths  in 
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a  direction  perpendicular  to  the  plane  of  fracture.  The  fracture  edge  was 
protected  by  an  edge  retainer  made  of  a  mixture  of  epoxiresin  and  silicon 
carbide  and/or  iron  powders  (600  mesh).  Mechanical  polishi?ig  was  followed 
by  electropolishing.  The  microhardness  measurements  were  made  with  a  micro  - 
durimet  Leitz  with  a  hardness  load  of  25  g  ,  except  for  the  warm-worked  301 
steel  which  required  o  50  g  load. 

The  fracture  surface-  of  fatigue  specimens  were  studied  both  by  scanning 
and  transmission  electron  microscopy,  Vhe  microscopic  growth  rates  were 
obtained  by  measuring  the  average  spacing  of  fatigue  striations  (i)  in  a 
direction  parallel  to  the  microscopic  crack  growth  direction,  For  striation 
spacings  lower  than  about  0.1  -  0.2  micron  all  the  measurements  were  made  on 
Cr-shadowed  carbon  replicas. 

2:4  X-ray  measurements 

The  amount  of  martensite  on  the  fracture  surface  of  301  fatigue  specimens 
tested  at  -30°C  and  25°C  was  determined  by  X-ray  diffraction.  A  comparison 
of  the  integrated  intensity  of  the  (311)  austenite  peak  with  the  (220) 
martensite  peak  was  made.  CoKa  radiation  was  used  and  the  standard  correc¬ 
tions  ware  employed.  (4) 

RESULTS 

1.  Transformation  Temperatures 

The  M$(y  a')  and  M^y  -*  a‘)  femperatures  are  given  in  Table  1.  The 
martensitic  transformation  y  -»•  e  (hep)  also  occurs  in  the  two  stainless  steels. 
In  alloy  16-13,  the  Ms(y  -*■  e)  temperature  coincides  with  Mg (y  -*•  a').  On  the 
other  hand,  in  an  alloy  very  similar  to  the  301  steel  studied  in  this  work, 
it  was  shown  by  Abrassart  (5)  that  the  y  -*■  e  martensitic  transformation  does 
occur  at  a  temperature  above  M^(y  -*•  a').  After  quenching  the  301  alloy  to 
room  temperature,  only  a  few  e  platelets  were  observed,  while  after  quenching 
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to  -196°C  (below  M$  (y  a')  )  ,  the  structure  of  the  301  steel  consists  of 
a  martensite  laths  imbedded  in  bands  containing  e  platelets  (Fig.  1). 

2.  Tensile  Properties 

Fe-NI  alloys 

Tensile  curves  at  25°C  are  reported  in  Figure  2.  The  high  yield 
strength  of  these  austenitic  alloys  has  been  attributed  to  the  ferromagnetic 
properties  of  these  compositions  (6).  After  quenching  to  -196°C,  the  amount 
of  residual  austenite  in  alloy  Fe-31 . 5Ni  is  low  and  the  tensile  curve  given 
in  Figure  2  is  comparable  to  that  of  a  material  which  is  nearly  completely 
martensitic.  The  phase  transformation  y  ot*  leads  to  an  increase  in  the 
yield  strength  and  a  decrease  in  the  elongation. 

Stainless  steels 

Tensile  curves  a-  e  and  magnetization  curves  J-e  of  301  steel  and 
16-13  alloy  are  given  respectively  in  Figures  3  and  4.  In  both  alloys,  the 
same  qualitative  behaviour  is  observed  when  the  temperature  is  decreased 
within  the  Ms~Md  range.  The  following  results  should  be  kept  in  mind: 

a)  the  elastic  limit  is  not  strongly  temperature  dependent  as  is  usually 
the  case  for  FCC  materials  tested  near  25°C. 

b)  a  strong  increase  in  the  hardening  rate  is  observed  as  the  deformation 
temperature  approaches  Ms>  At  25°C,  the  monotonic  strain  hardening  rate  of 
301  steel  is  higher  than  that  of  16-13  alloy,  mainly  because  the  volume  of 
strain-induced  a'  martensite  at  a  given  plastic  strain  is  higher  at  this 
temperature  in  the  301  steei  than  in  the  16-13  matrix. 

c)  At  low  temperatures,  e.g.  -45°C  for  16-13  alloy  and  -30°C  for  301 
steel,  for  the  same  volume  of  transformed  a'  phase,  the  flow  stress  of  301 
steel  is  higher  than  that  of  16-13  alley.  This  difference  is  mainly  due  to 
the  difference  in  carbon  content  of  the  two  alloys. 
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d)  The  rapid  change  in  the  work  hardening  exponent  with  decreasing 
temperatures  in  the  Md-M  range  results  in  large  variation*  of  the  uniform 
elongation  (Fig.  5).  This  behaviour  is  in  agreement  with  results  on  stain¬ 
less  steels  (7,7). 

The  tensile  properties  of  301  steel,  in  warm-worked  condition,  are 
reported  in  Table  2.  At  25°C.  the  y  -*■  a'  transformation  results  in  a  lower 
work  hardening  exponent  for  the  warm-worked  alloy  than  for  the  annealed  alloy. 
At  25° C,  the  typical  UTS/YS  ratios  are  1.23  for  the  warm-worked  condition 
and  3.78  for  the  annealed  alloy. 

3.  Fatigue  Properties  '• 

3:1  Fatigue  crack  growth  rates  (FCGR) 

FCGR  in  austenite  compared  with  FCGR  in  martensite 
The  crack  growth  rates  da/dN  versus  A  K  of  the  Fe-Ni  alloys  at  20°C  are 
reported  in  Figure  6.  The  FCGR  in  Fe-34  alloy  which  has  a  higher  yield  * 
strength  than  that  of  Fe-31.5  steel  are  slightly  lower.  In  the  Fe-31.5  alloy 
transformed  to  martensite  by  quenching  to  -196°C,  the  FCGR  are  an  order  of 
magnitude  faster  than  for  the  austenite,  for  A  K  <  25  ksi  /fn  .  Moreover, 
the  slope  of  the  da/dN  -  A  K  curve  is  lower  for  the  martensite  than  for  the 
austenite  phase.  The  macroscopic  crack  growth  rates  measured  in  martensite 
and  the  slope  of  the  da/dN  -  A  K  curve  are  in  good  agreement  with  those  rates 
reported  by  Bathias  (8)  and  van  Swam(9)  for  Maraging  steels  which  are  also 
Fe-Ni  martensites 

FCGR  and  strain-induced  martensitic  transformations 
The  FCGR  of  the  16-13  alloy  at  three  different  temperatures  in  the  Ms~Md 
range  are  given  in  Figure  7.  The  propagation  rates  decrease  with  decreasing 
temperature.  No  significant  change  in  the  slope  of  the  da/dN  -  AK  curves 
are  observed  in  the  range  of  temperatures  and  AK  which  were  studied. 
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In  annealed  and  warm-worked  301  steel,  the  results  of  da/dN  measure¬ 
ments  are  reported  in  Figure  8.  A  strong  temperature  dependence  is  observed 
between  95°C  (=Mj)  and  -30°C  for  the  annealed  alloy.  The  FCGR  and  the  slopes 
of  the  da/dN  -AK  curves  decreased  when  the  temperature  approaches  M,..  For 
A  K  >  35  ksi  /Tn,  the  FCGR  at  30CC  is  an  order  of  magnitude  slower  than  at 
95°C.  Tests  in  air  gave  a  propagation  rate  which  was  significantly  higher 
than  in  dry  argon.  This  difference  could  be  due  to  the  influence  of  humidity 
on  crack  growth  rate  in  transformed  a’  martensite  present  at  the  crack  tip. 
However,  it  is  worthwhile  to  note  that  Shahinian  et  al  (10)  have  also 
reported  an  increased  growth  rate  due  te  humidity  in  air  in  a  304  stainless 
steel,  tested  at  25°C,  although  their  alloy  was  mere  stable  than  the  301  alloy 
investigated  here. 

The  fatigue  behaviour  of  the  301  alloy  in  the  warm-worked  and  in  the 
annealed  conditions  is  markedly  different.  Even  without  any  large  amount  - 
of  of  strain-induced  phase  transformation  (6  =  85°C),  da/dN  is  significantly 
higher  for  low  A  K  (AK  <  40  ksi  /in'  )  in  the  previously  warm-worked  struc¬ 
ture.  The  results  obtained  with  warm-worked  301  alloy  are  compared  in  Figure 
8  with  those  results  of  Chanani  et  al  (11)  on  another  TRIP  steel  which  was 
subjected  to  a  thermomechanical  treatment  similar  to  that  applied  to  the  301 
steel.  These  authors  performed  all  their  tests  in  air  on  relatively  thin 
sheet  specimens  (=  1.90  mm  thick). 

3:2  Plastic  zone  sizes  measurements 

Figure  9  shows  two  typical  plots  of  hardness  versus  distance  from  the 
fracture  surface,  fer  the  Fe-34Ni  alloy.  Each  ha.  tiness  curve  shows  two 
plateaus  which  clearly  define  the  boundaries  of  the  two  plastic  zones  usually 
referred  to  as  the  monotonic  (Rm)  and  the  reversed  (Rr)  plastic  zones. 

Microhardness  measurements  were  also  used  for  the  301  steel  both  in 
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annealed  and  warm-worked  conditions  (Table  3).  The  annealed  301  steel  shows 
a  similar  behaviour  to  that  of  the  Fe-34Ni  alloy.  However,  the  heterogeneity 
of  the  phase  transformations  at  the  crack  tip  gives  a  large  scatter  of 
measurements.  The  high  yield  strength  of  the  warm-worked  steel  leads  to 
smaller  plastic  zone  sizes  and  the  reversed  plastic  zones  could  not  be 
.  unambiguously  identified  for  Rm  <  0.2  mm.  Moreover,  when  there  is  no  strain  - 
Induced  y  -*•  <*'  transformation,  the  warm-worked  alloy  cyclically  softens. 

3:3  Metallography 

Optical  microscopy  and  X-ray  observations 

Optical  metallography  showed  a  qualitative  agreement  between  the  struc¬ 
tures  observed  in  tensile  and  fatigue  specimens.  Platelets  present  near  the 
fracture  surfaces  at  95°C  are  twins  and/or  e  martensite  (Fig.  10a).  Testing 
with  a  strong  magnet  shewed  that  the  fracture  surfaces  of  specimens  tested 
at  this  temperature  were  not  magnetic.  On  the  other  hand,  at  -30°C,  the  * 
amount  of  phase  transformation  inside  the  plastic  zone  is  significant  (Fig. 
10b).  X-ray  measurements  s^wed  that  the  amounts  of  a'  martensite  were  re¬ 
spectively  =100%  at  -30°C  and  =50%  at  25°C. 

Fractography 

Fractographic  investigation  was  limited  to  the  stainless  alloys.  These 
steels  exhibited  flat  transgranular  fracture  in  all  specimens  tested.  The 
fracture  surfaces  were  much  rougher  at  low  AK  than  at  high  A  K,  while  the 
fracture  surface  inside  a  given  grain  is  flatter  at  high  A  K.  Microscopic 
fracture  features  which  could  be  associated  with  the  martensitic  transfor¬ 
mations  were  observed  (Fig.  11a,  b). 

The  main  difference  between  the  two  steels  was  the  tendancy  to  form 
secondary  cracks.  In  301  alloy,  secondary  cracks  were  occasionally  observed 
only  at  high  A  K.  On  the  other  hand,  in  16-13  alloy,  even  for  low  AK 
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(*  25  ksl  fin  ),  a  high  density  of  secondary  cracks  was  present  (Fig.  12a,  b). 
In  the  301  steel,  at  high  growth  rates  (e.g.  at  95°C,  AK  *  40  ksi  /fn  ),  the 
fracture  surfaces  exhibited  extensive  tear  dimples. 

Striation  spacings  measurements  for  the  two  stainless  steels  are  reported 
in  figures  13  and  14.  Each  datum  point  corresponds  to  the  mean  value  of  the 
•distance  between  at  least  20  striations.  Measurements  were  taken  with  the 
SEM  for  i  >0.2  micron  and  a  good  agreement  between  measurements  made  by  SEM 
and  electron  microscopy  of  replicas  was  obtained  in  this  range  of  growth  rates. 

In  both  alloys,  it  is  observed  that  the  microscopic  crack  growth  rates 
decrease  with  the  macroscopic  rateswith  decreasing  temperatures.  The  slopes 
of  the  1  -  AK  curves  increase  when  the  amount  of  strain  induced  a'  martensite 
Increases,  and  the  slopes  of  these  curves  are  quite  high  (>  5)  even  when 
there  is  no  y  -*•  a*  phase  transformation. 


DISCUSSION 


Plastic  Zone  Sizes  -  Cyclic  Hardening  and  Softening  at  the  Crack  Til 


Microhardness  measurements  at  different  depths  below  the  fracture  sur¬ 


faces  gave  an  excellent  estimation  of  the  monotonic  and  reversed  plastic  zone 
sizes.  The  results  reported  here  are  similar  to  the  results  obtained  by 
Bathias  (8).  The  monotonic  plastic  zone  sizes  which  are  reported  in  Figure 


15  are  in  good  agreement  with  the  theoretical  sizes  given  by  the  equation  R 

A  K  2  m 

c  A(- — )  .  The  values  of  a  were  obtained  from  low  strain  rate  tensile 

°y  y 

tests.  The  experimental  value  of  0.053  for  A  is  one  half  of  the  value  found 
by  Hahn  et  al  for  a  Fe-3Si  alloy  (12)  and  one  third  of  the  theoretical  value 
of  0.157  given  by  Levy  et  al  (13)  for  the  plane  strain  plastic  zone  size  of 
a  non-hardening  plastic  material.  The  difference  between  the  measured  and 


theoretical  values  of  A  could  be  related  to  the  high  monotonic  workhardening 
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exponent  of  the  stainless  steels.  When  the  reversed  plastic  zone  was  large 
enough  to  be  measured  by  microhardness  indentations,  the  zone  was  approxi¬ 
mately  one  fourth  of  the  size  of  the  monotonic  plastic  zone  which  is  in 
agreement  with  Rice's  modci  of  reversed  yielding  at  the  tip  of  a  fatigue 
crack.  (14) 

The  only  alloy  which  exhibited  cyclic  softening  at  the  crack  ti.  was 
the  warm-worked  301  stainless  steel  at  a  test  temperature  near  Md-  This  soft¬ 
ening  behaviour  could  be  predicted  from  the  value  of  the  ratio  of  UTS  to  Oy. 

At  0  ■  95°C  (Table  2),  this  ratio  is  1.05,  which  is  lower  than  the  critical 
value  of  1.2  for  cyclic  stability  (15).  *At  25°C,  the  same  alloy  undergoes 
cyclic  hardening  because  the  martensitic  transformation  at  the  crack  tip 
offsets  the  cyclic  softening  of  the  warm-worked  austenite.  The  difference 
between  cyclic  softening  at  95°C  and  cyclic  hardening  at  25°C  is  not  reflected 
in  a  marked  change  in  the  slopes  of  the  da/dN  -  A  K  curves  (see  rig.  8).  - 
2.  Crack  Growth  Rate  in  Stable  Austenite 

2:1  Influence  of  SFE  and  deformation  modes  on  FCGR 

By  comparing  the  FCGR  of  different  copper  base  alloys,  some  authors 
(16-19)  have  concluded  that  a  decrease  in  intrinsic  stacking  fault  energy 
(SFE)  will  lower  the  fatigue  crack  growth  rate.  Since  the  SFE  of  the  austen¬ 
itic  steels  used  here  were  known,  it  is  worthwhile  to  compare  the  FCGR  of 
the  alloys  on  the  basis  of  SFE  levels. 

The  austenitic  alloys  studied  here  can  be  conveniently  grouped  in  two 
classes,  depending  upon  their  SFE: 

2 

(i)  the  SFE  of  the  Fe-Ni  alloys  is  greater  than  50  ergs/cm  and  they 
deform  only  by  slip.  The  yield  strength  of  these  ferromagnetic  austenites 
is  relatively  high. 

(ii)  the  SFE  of  the  two  stainless  steels  is  much  lower  than  the  Fe-Ni 
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alloys.  It  Is  of  the  order  of  25  -  30  ergs/cm  at  (1,5).  This  low  value 
for  the  SFE  results  in  plastic  deformation  by  a)  twinning  at  90°C  for  the 
16-13  alloy;  b)  twinning  and  y  e  martensitic  transformation  at  S5°C  for 
the  301  steel. 

The  FC6R  of  the  Fe-34Ni  and  the  annealed  301  alloy  are  comparable, 
although  the  SFE  and  deformation  modes  ere  quite  different.  On  the  other 
hand,  a  comparison  of  the  FCGR  of  the  two  Fe-Ni  alloys  and  the  two  annealed 
stainless  steels  shows  that  the  FCGR  decreases  with  increasing  yield  stress 
(Fig.  16).  It  is  well  known  that  this  correlation  is  not  observed  in 
quenched  and  tempered  high  strength  steels  such  as  maraging  steels  (8,9). 
which  usually  work  soften  at  the  crack  tip.  However,  the  warm-worked  301 
steel  which  work-softens  at  85°C  shows  a  higher  FCGR  than  the  annealed  mater¬ 
ial  although  their  yield  strengths  differ  by  a  factor  of  5. 

It  is  tentatively  concluded  that  the  FCGR  decreases  with  increasing 
tensile  yield  stress  as  long  as  the  alloy  undergoes  cyclic  hardening  or  is 
cyclically  stable. 

2:2  Macroscopic  and  microscopic  FCGR 

A  detailed  comparison  between  da/dN  and  the  striation  spacing  (i)  was 
made  for  the  two  stainless  steels.  The  crack  propagation  rates  (da/dNl 
measured  in  these  two  alloys  near  the  temperature  are  similar  to  those 
reported  for  different  stable  commercial  stainless  steels  tested  at  25°C 
(8,  10,  20,  21  ). 

The  FCGR  behaviour  of  301  steel  is  typical,  that  is  at  low  A  K,  i  <  da/dN 
and  i  >  da/dN  at  high  A  K.  At  low  A  K,  the  difference  between  i  and  da/dN 
shows  that  crack  growth  is  discontinuous  along  the  crack  front.  At  high 
A  K,  static  fracture  modes  such  as  ductile  void  formation  accelerate  the 
growth  rate^and  the  transition  where  the  i  and  da/dN  versus  AK  curves  cross 
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over  corresponds  to  the  appearance  of  voids  on  the  fracture  surfaces.  In 
16-13  alloy,  it  Is  more  difficult  to  compare  the  microscopic  and  macroscopic 
growth  rates  because  of  the  extensive  density  of  secondary  cracks. 

It  is  worthwhile  to  nDte  that  the  same  1  -  AK  dependence  was  found 
for  both  alloys  (Fig.  13,  14  )  and  that  the  slooes  of  the  i  -  AK  curve 
are  greater  than  the  value  of  2  which  is  predicted  by  a  simple  crack  growth 
model  based  on  the  crack  opening  displacement  (COD)  concept.  Slopes  greater 
than  2  have  also  been  measured  by  Weber  and  Hertzberg  (21)  in  a  305  stain¬ 
less  steel.  Moreover,  the  experimental  values  for  i  are  at  least  one  order 
of  magnitude  lower  than  the  values  predicted  by  the  COD  model,  that  is  i  * 

A  K2 

0.25  —  ,  taking  a  cyclic  yield  strength  a  of  115  ksi  for  a  301  steel 

toyc  yc 

(22). 

3.  Influence  of  Strain  -induced  Martensitic  Transformation  on  FCGR 

Since  the  influence  of  the  strain-induced  phase  transformation  on  the 
FCGR  was  studied  by  varying  the  test  temperature,  we  must  first  determine 
the  effect  of  a  temperature  change  on  the  mechanical  and  physical  properties 
which  could  affect  the  FCGR. 

For  stainless  steels,  the  variation  of  the  SFE  (y)  with  temperature  is 
given  by  dy/dT  -  -0.10  ergs/cm  °  K  (5,23,24).  Thus,  in  the  relatively  narrow 
temperature  range  used  here,  the  value  of  y  changes  by  a  factor  less  than  2. 
Furthermore,  as  discussed  above  ,  it  is  felt  that  the  SFE  is  not  a  critical 
factor  in  controlling  the  PCGR.  The  Young  modulus  (E)  will  increase  by  about 
3%  for  a  decrease  in  temperature  of  100  to  150°C.  This  variation  of  E  will 
not  markedly  affect  the  COD  or  the  elastic  strain  distribution  near  the 
crack  tip. 

It  can  therefore  be  concluded  that  the  change  in  FCGR  with  temperature 
is  unambiguously  due  to  the  strain  induced  martensitic  transformation.  The 


263< 


temperature  dependence  appears  to  be  related  to  the  relative  amount  of  marten¬ 
site  formed  near  the  crack  tip.  Metallographic  observations  and  X-ray  diffrac¬ 
tion  analysis  of  the  plastic  zone  show  that  the  amount  of  a'  martensite 
increases  with  decreasing  temperature.  It  remains  then  to  explain  how  the 
phase  transformation  will  lead  to  a  lower  FCGR  and  to  a  decrease  in  the  da/dN 
-  A  K  slope  with  decreasing  temperature.  This  change  of  slope  was  observed 
only  in  the  301  stainless  steel  which  exhibited  the  larger  volume  fraction  of 
transformed  a 4 . 

The  lower  growth  rates  with  a  greater  amount  of  phase  transformation 
are  not  due  to  the  fact  that  the  FCGR  in  partensitt  is  lower  than  the  FCGR 
in  austenite.  On  the  contrary,  the  results  on  Fe-31.5Ni  alloy  showed  that 
the  FCGR  in  the  fully  austenitic  condition  was  lower  than  in  the  fully  marten¬ 
sitic  condition.  However,  In  this  alloy,  the  slope  of  the  da/dN  -  AK  curve 
Is  lower  in  martensite  than  in  austenite.  This  difference  in  slopes  could, 
explain  the  lower  slopes  observed  in  3C1  steel  when  the  test  temperature  is 
decreased. 

The  influence  of  the  strain-induced  phase  transformation  on  the  FCGR  ran 
best  be  explained  by  considering  the  concomittant  changes  of  the  monotonic 
and  cyclic  work  hardening  coefficients.  The  marked  increase  in  the  cyclic 
work  hardening  coefficient  due  to  y  -*-a‘  phase  transformation  was  reported 
previously  (22).  Currently,  there  is  no  complete  FCGR  theory  which  takes  into 
account  the  influence  of  both  the  monotonic  (n)  and  the  cyclic  (n1)  work 
hardening  coefficients.  The  theoretical  models  of  the  da/dN  -  AK  laws  can 
be  conveniently  grouped  in  two  broad  classes;  a)  the  work  hardening  or  inte¬ 
grated  damage  theories,  and  b)  the  COD  theories  based  on  a  plastic  deforma¬ 
tion  and  fracture  criteria  at  the  crack  tip.  In  all  these  theories,  higher 
monotonic  nnd  cyclic  hardening  coefficients  will  lower  the  FCGR. 
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The  theory  of  Head  (25),  which  represents  the  plastic  2one  at  the  crack 
tip  by  rl <j id  work  hardening  elements, shows  da/dN  varylny  as  nv~  —  .  In  the 

„  '  "y 

301  steel,  this  ratio  Increases  by  a  factor  of  3  between  95°C  and  -30°C, 

while  the  growth  rate  decreases  by  a  factor  of  20  for  A  K  >  40  ksi. 

It  can  also  be  shown  that  an  increase  of  the  monotonic  and  cyclic  work 

hardening  coefficients  will  lead  to  a  decrease  of  the  strain  gradient  in  the 

•  plastic  zone  near  the  crack  tip.  For  instance,  for  a  work  hardening  material 

a  c  Ae_r,  the  monotonic  strain  distribution  ir.  the  plastic  zone  will  be 
P  R  1/1-n 

given  by  ep  -  ey(-  -  1)  #  where  R  is  the  plastic  zone  size  and  r  the 

distance  to  the  crack  tip  (26).  Due  to  the  fact  that  n*  increases  as  well 

as  n,  the  cyclic  plastic  strains  distribution  within  the  plastic  zone  will 

also  be  more  uniform.  Thus,  the  integrated  damage  criteria  proposed  by 

McClintcck  (27)  should  lead  to  a  decreasing  FCGk  when  n  and  n’  increase. 

Although  the  CGD  model  cannot  accurately  account  for  the  observed  propa- 

« 

gation  rates  (cf.  2:2),  it  is  worthwhile  to  use  this  concept  to  estimate 

qualitatively  the  influence  of  the  martensitic  transformation  on  the  cyclic 

yield  stress  -od  on  the  work  hardening  exponent.  The  phase  transformation  at 

the  crack  tip  results  in  a  marked  increase  in  the  cyclic  elastic  limit  which  lowers 

the  COD  at  a  given  AK  level.  For  instance,  the  cyclic  stress  strain  curves 

of  a  301  type  steel  showed  that  ayc  at  -30°C  is  3  to  6  times  higher  at  -30°C 

than  at  90°C.  a  depends  strongly  upon  the  value  of  the  perma- 

yc 

nent  offset  plastic  strain  at  which  it  is  measured.  McClintock  also  showed 
that, theoretically,  the  COD  decreases  when  the  monotonic  work  hardening 
coefficient  increases  (28). 

However,  these  predicted  changes  of  the  FCGR  due*  to  the  variation  of  a 
and  of  the  work  hardening  coefficients  are  not  as  large  as  the  obser.ed 
variations  of  FCGR  with  decreasing  temperatures,  particularly  at  high  AK 
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levels  It  Is  felt  that  the  residual  stresses  near  the  crack  tip  and  along 
the  flanks  of  the  crack  behind  the  tip  may  play  an  important  role  In  control¬ 
ling  the  FCGR  when  strain  induced  y  a'  phase  transformation  is  taking  place. 
The  role  of  residual  stresses  has  been  demonstrated  by  the  work  of  Elber,  who 
Introduced  the  concept  of  effective  stress  intensity  factor  (29,30).  At 
this  time,  we  do  not  have  any  theoretical  approach  to  predict  the  distri¬ 
bute  of  residual  stresses.  We  can  only  say  qualitatively  that  the  high 
cyclic  yield  and  ultimate  tensile  stresses  -lilting  from  the  a'  phase  trans¬ 
formation  could  lead  to  very  high  residual  stresses  which  would  lowe*"  the 
value  of  A  K  effective  for  a  given  AK  applied.  These  residui.l  stresses 
appear  to  be  very  stable  and  to  resist  any  cyclic  shakedown  due  to  reversed 
plastic  deformation  within  the  plastic  zone. 

CONCLUSIONS 

1.  The  FCGR  of  metastable  stainless  steels  were  found  to  be -strongly  temper¬ 
ature  depsndert  in  the  M$ -Mq  temperature  range.  At  a  given  AK  level,  the 

,  FCGR  decrease  as  t  e  volume  fraction  of  strain-induced  y  a'  phase  transfor¬ 
mation  at  the  crack  tip  increases  with  decreasing  temperature. 

2.  In  the  stable  austenites,  the  FCGR  do  not  depend  upon  the  intrinsic 
stacking  fauH  energy  or  upon  the  deformation  modes  of  the  austenites.  It  is 
found  that  in  the  annealed  austenites,  the  FCGR  will  decrease  wnen  the  mono¬ 
tonic  tensile  yield  stress  increases. 

3.  In  a  Fe-31 - 5Ni  alloy,  the  FCGR  is  higher  in  the  martensitic  strucsuie 
than  in  the  annealed  austenite  for  A  K  <  25  ksi  vTn. 

4.  A  3C1  type  stainless  steel  has  a  lower  fatigue  crack  growth  resistance 
in  the  high  strength  warm-worked  condition  than  in  the  annealed  austenitic 
condition  for  A  K  <  40  ksi  /Tn. 

5.  The  monotonic  and  reversed  plastic  zone  sizes  which  were  determined  by 
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mlcroliiii dnr<;«i  nipnMiiTiiiPiiti  nrp  In  pond  agreement  with  the  sizes  predicted  from 
thence l  li a  I  model-.. 

6.  The  influence  of  the  strain-induced  martensitic  phase  transformation  on 
the  FCGR  in  the  Ms-Md  range  is  attributed  to  the  strong  concomittant  varia¬ 
tions  of  cyclic  and  monotonic  work  hardening  exponents  and  to  residual  stresses 
left  near  the  crack  tip  and  along  the  flanks  of  the  fatigue  fracture  surfaces. 


ACKNQWirnr.LMI.NT5 

Mil*  work  Wdb  supported  by  UAHC  li»-7JC-G316. 

The  X-ray  measurements  were  performed  by  Professor  Van  Der  Sando, 
Dept,  of  Metallurgy,  M.I.T. 


9* 


4*  4 


REFERENCES 


1.  F.  Lecroisey  and  A.  Pineau,  Met.  Trans,  1972,  Vol.  3,  p.  387. 

2.  F.  Abrassart,  F.  Lecroisey  and  A.  Pintail,  Mem.  Sci.  Rev.  Met.. 1969. 

Vol.  66,  p.  805. 

2.  ASTM  Standard  E  99  -  70T. 

•4.  B.D.  Cullity,  Elements  of  X-ray  Diffraction,  Addison-Wesley  Co.  (1967), 
p.  388. 

5.  F.  Abrassart,  Thesis  (Sc.D),  University  of  Nancy,  1972. 

6.  6.F.  Bolling  and  R.H.  Richman,  Phil.  Mag.,  1969,  Vol.  19,  p.  247. 

7.  J.P.  Brtssanelli  ar.d  A.  Moskowitz,  ASM  Trans  Quart.,  196G,  Vol.  59.  p.  223 

8.  C.  Bathias  and  R.M.  Pelloux,  Met.  Trans,  1973,  Vol.  4,  p.  1265. 

9.  L.F.  van  Swam,  Thesis, (Sc.D),  M.I.T.,  Cambridge,  1973. 

10.  P.  Shahinian,  H.E.  Watson  and  H.H.  Smith,  Ja1  of  Materials,  1972,  p.  527. 

11.  G.R.  Chanani,  D.S.  Antolovich  and  W.W.  Gerberich,  Met.  Trans,  1972,  Vol. 3, 

p.  2661. 

12.  G.T.  Hahn,  M.  Sarrate  and  A.R.  Rosenfield,  Inta^  cf  Fract.  Mech.,  1971 

Vol.  7,  p.  435. 

13.  N.  Levy,  P.V.  Marcal,  W.J.  Ostergren  and  J.R.  Rice,  Inta^  Ja^of  Fract. 

Mech. ,  1971 „  Vol.  i,  p.  142. 

14.  J.R.  Rice,  ASTM  STP  415,  1966,  p.  247. 

15.  R.W.  Smith,  M.H.  Hirschberg,  S.S.  Manson,  NASA  TND  -1574,  1963. 

16.  A.J.  McEvily  and  R.C.  Boettner,  Acta  Met.,  1963,  Vol.  11,  p.  725. 

17.  G.A  Miller,  D.H.  Avery  and  W.A.  Backofen,  T^ans,  AIME,  1966,  Vol.  236, 

p.  1667. 

18.  H.  T shi i  and  J.  Weertman,  Met.  Trans.,  1971,  Vol.  2,  p.  3441. 

19.  C.M.  Kan  and  J.G.  Byrne,  "Fracture",  Brighton  Conf.  1969,  p.  598. 

20.  L.E.  James  and  E.B.  Schwenk,  Met.  Trans.,  1971,  Vol.  2,  p.  491. 


21.  J.H.  Weber  and  R.W.  Hertzberg,  Met.  Trans. ,1973,  Vol.  4,  p.  595. 


22.  A.  Pineau,  L.F.  van  Swam  and  R.M.  Pelloux,  Scrlpta  Met..  1973,  Vol.  7, 

p.  657. 

23.  F.  Lecroisey  and  B.  Thomas,  Phys.  Stat.  Sol,  (a),  1970,  Vol.  2j  K.217 

24.  R.M.Latanisonand  A.W.  Ruff,  Met.  Trans.,  1971,  Vol.  2,  p.  50E 

25.  A.K.  Head,  Phil.  Mag.,  195:,  Vol.  44,  Ser.  7,  p.  925. 

26.  J.R.  Rice  and  G.F.  Rosengren,  Ja^  Mech.  Phys.  Solids,  1968,  Vol.  16,  p.l. 

27.  F.A.  McClintock,  Fractjre  of  Solids.  John  Wiley,  New  York,  1963,  p.  65. 

28.  F.A,  McClintock,  Fracture,  Academic  Press,  New  York  and  London,  1971, 

Vol.  3,  p.  47.  * 

29.  W.  Elber,  ASTM  STP  486,  p.  230. 

30.  W.  Elber,  Eng.  Fracture  Mechanics,  1970,  Vol.  2,  p.  37. 


<(-  t  o< 


Alloy  compositions  (weight  %)  -  Heat  treatments  - 
Austenite  grain  sizes  -  Transformation  temperatures 


Type 


Fe-34Ni  34.09 
Fe-31 .5Ni  31.41 


Composition 


0.004 


6.50  18.25  0.19 


M  (°C)  M.  CC)  Water-  grain  size 

"  a  quenched  (microns) 


0.005  Bal  <-196 


16-13  12.50  15.8  0.010 


1100°C 
(  1  hr) 

1120°C 
(  1  hr) 


ansile  properties  of  alloy  301  warm-worked  65%  at  450°-500°C 


Test 

temperature  (°C) 


Yield 

irengtf 


148 


Ultimate  tensile 
strength  (ksi) 


154.6 


Total 

Elongation 

Test  direction 

15.1 

longitudinal 

45.5 

40.3 

longitudinal 

transverse 

Figures 

Fig.  1  -  Alloy  301:  a'  martensite  imbedded  in  e  platelets  bands  following 
quenching  to  -196°C. 

Fig.  2  -  Fe-Ni  alloys:  tensile  stress-strain  curves  at  25°C. 

Fig.  3  -  Alloy  301:  tensile  stress-strain  curves  and  magnetization  (J)  - 
strain  curves  at  different  test  temperatures  -  (100%  martensite  corresponds 
to  J  -  14,000  gauss  ). 

Fig.  4  -  Alloy  16-13:  tensile  stress-strain  curves  and  magnetization  (J)  - 
strain  curves  at  different  test  temperatures  -  (100% martensite  corresponds 
to  J  =  14,000  gauss). 

Fig.  5  -  Elongation  <5  versus  test  temperature:  uniform  tensile  elongation  ver 
sus  test  temperature  for  alloys  16-13  and  301. 

Fig.  6  -  Fe-Ni  alloys  :  fatigue  crack  growth  rates  as  a  function  of  A  K  at 
20°C. 

Fig.  7  -  Alloy  16-13  :  fatigue  crack  growth  rates  as  a  function  of  A  K  at 
different  temperatures. 

Fig.  8  -  Alloy  301  :  fatigue  crack  growth  rates  of  annealed  and  warm-worked 
301  steel  at  various  temperatures.  At  0  =  20°C,  the  warm-worked  301  alloy 
was  tested  in  the  transverse  direction  (TD)  and  longitudinal  direction  (ID). 
Fig.  9  -  Alloy  Fe-34Ni:  plastic  zone  size  determination  by  microhardness 
measurements  at  two  AK  levels. 

Fig.  ]0  -  Alloy  301  :  metal lographic  evidence  of  the  reversed  plastic  zone 

size  below  the  fracture  surface.  (a)0  =  95°C,  AK  =  26.7  ksi  /Tn  ;  (b)  AK  = 

_  e  =  -30°c. 

42.6  ksi  /TrijV  The  edges  of  the  plastic  zones  are  indicated  by  an  arrow. 

Fig.  11  -  Alloy  301  :  typical  fracture  lamellae  which  are  associated  with  the 
transformed  martensitic  structures  (6  =  25°C).  (a)  AK  =  43  ksi  /in;  (b)  AK 
26.5  ksi  /Tn. 


Z73< 


Fig.  12  -  Alloy  16-13  :  typical  fracture  features  at  90°C.  The  density  and 
residual  opening  of  the  secondary  cracks  are  higher  in  fig.  12b  (AK  =  31.6 
ksi  /fn  )  than  in  fig.  12a  ( A  K  ®  23.3  ksi  /Tn  ). 

Fig.  13  -  Alloy  16-13  :  comparison  of  the  macroscopic  and  microscopic  (i) 
crack  growth  rates  at  90°  and  20°C. 

Fig.  14  -  Alloy  301  :  comparison  of  the  macroscopic  and  microscopic  (i)  crack 
growth  rates  at  95°  and  20°C. 

Fig.  15  -  Measurements  of  monotonic  plastic  zone  sizes  (Rm)  versus  ( A  K/ay)  . 

Fig.  16  -  Fatigue  crack  growth  rates  of  different  stable  austenites  as  a  function  of  ZK. 
The  tensile  yield  strength  of  each  alloy  is  given  for  reference. 
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+  Striation  spacing  9  =  90°C 
•  Striation  spacing  9  -  20°C 
Macroscopic  Crack  Growth  Rates 


Stria*ion  Spacing  8 
Striation  Spacing  6 
Macroscopic  Crack 
Growth  Rate 
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SEMI-ANNUAL  TECHNICAL  REPORT 


July  1  -  Dec.  31,  1973 
George  B.  Benedek 

During  the  period  July  1  -  Dec.  31  our  work  has  been  devoted  to  three  major 
areas. 

a)  Characterization  o£  the  viscoelastic  properties  of  gels  froe  the  quasi- 
elastic  spectrum  of  the  scattered  laser  light. 

b)  The  measurement  of  tb”  equation  of  state  of  a  two-dimensional  monolayer 
"gas". 

c)  The  study  of  the  onset  and  structure  of  turbulence  from  the  spectrum 
of  light,  microwaves  or  sound  waves  scattered  from  the  fluid. 

We  shall  report  here  on  the  progress  of  each  of  these  in  turn. 

a)  Dr.  T.  Tanaka  and  D.  L.  Hockar  have  succeeded  in  measuring  the  spec¬ 
trum  of  laser  light  scattered  from  poly-acrylaaide  gels  as  a  function  of  the 
gel  density,  scattering  angle,  and  the  temperature.  They  have  been  able  to 
interpret  their  data  in  terms  of  a  theory  whose  principal  parameters  are  else 
elastic  constant  of  the  gel  fiber  network,  and  the  friction  factor  between  the 
fiber  network  and  ths  ^el  liquid.  In  fact  from  their  measurements  they  mere 
able  to  measure  the  ratio  of  these  quantities,  and  the  temperature  dependence 
of  this  ratio.  By  means  of  painstaking  mechanical  measurements  they  were 
able  to  measure  macroscopically  both  the  elastic  constant  and  friction  fae'er 
and  were  able  to  demonstrate  that  the  values  so  determined  were  in  excellent 
agreement  with  those  found  from  the  light  scattering  spectra.  They  also 
showed  that  the  temperature  dependence  of  the  elastic  constant ,  as  obtained 
from  the  light  scattering  measurements  was  not  quite  linearly  proportional 
to  the  temperature.  This  measurement  of  the  departure  from  linearity  Km 
important  consequences  for  the  theory  of  randos  polymer  chains.  This  wore 
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has  keen  written  up  2nd  has  appeared  in  the  Journal  of  Chemical  Physics,  59, 
5151,1573.  Dr.  Hocker  is  cow  at  the  Lincoln  Laboratories  of  M.I.T. 

b)  Dr.  Gilbert  Hawkins  has  measured  the  pressure-area  isotherms  of  a  two- 
dimensional  gas  of  pentaaecylic  acid  molec  '«•■?  resting  on  the  surface  of  water. 

Se  designed,  and  brought  into  operation,  an  apparatus  for  this  purpose  whose 
sens!  rvity  exceeded  by  a  factor  10  previou-  studies  of  surface  monolayers. 

Viu  this  apparatus  Dr.  Hawkins  demonstrated  for  the  fiist  time  the  existence 
of  a  critical  point  in  a  gas-liquid  type  of  transition,  for  a  two-dimensional 
system.  Ee  also  measured  tie  shape  of  the  "gas-liquid"  coexistence  curve  and 
the  «p^oent  which  inscribes  the  divergence  oi  the  compressibility.  The  re- 
malts  were  interpreted  in  tens  of  a  mean  field  theory  which  succeeds  in  pre- 
dictinj  the  exponents  characterising  the  divergence  of  the  compressibility  and 
tM  shape  cf  the  coexistence  curve.  The  mean  field  description  is  consistent 
with  the  existence  e*  long-range  interactions  between  the  polar  molecules  of 
peataaecylic  acid.  This  work  represents  the  first  measurement,  of  the  equation 
of  state  of  a  two-dimensional  monolayer  and  should  stimulate  the  investigation 
of  the  rich  range  of  phase  transitions  which  are  possible  in  such  two-dimensional 
systems  - 

This  work  has  just  keen  wrictei  up  for  publication  as  a  Physical  Review 
Letter.  The  manuscript  was  submitted  in  Dec.  1973.  Dr.  Hawkins  is  now  a 
SHli-r  Instil  e  Fellow  in  Physics  at  the  Bniveicity  o.  ilifornia  at  Berkeley. 

cl  Mr.  lies  Eoss  and  Professor  G.  Senedek  arc  now  reviewing  the  literature 
os  experimental  end  theoretical  studies  of  'he  transition  between  laminar  and 
Cwrh&Iest  field  flow.  Since  the  development  of  light  scattering  spectroscopy, 
it  haa  beet  dear  that  *  new  technical  instn*r*ent  now  exists  for  the  accurate 
•easuresKst  of  velocity  and  velocity  correlations  in  fluid  flow.  That  instrument 
:s  tSe  dec*  ws  of  Icppler  shifts  in  the  flcwing  fluid  Ly  the  methods  of  optical 
eJlsj  1 y  croscopy.  The  subject  of  the  Laminar- turbulent  transition  is  one  of 
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great  practical  and  theoretical  importance.  We  are  now  endeavoring  to  deter¬ 
mine  what  is  known  about  the  complex  phase  boundaries  of  the  transition,  and 
the  characteristics  of  the  flow  whose  measurement  is  of  principle  importance. 

At  this  stage  we  are  in  a  literature  search  and  discussion  phase. 

Report  Summary 

In  the  period  July  1  -  Dec.  31  we  have  completed  two  separate  investigation 
The  first  is  a  study  of  the  spectrum  of  light  scattered  quasi-elastically  from 
gels.  The  second  is  the  accurate  measurement  of  the  equation  of  state  of  a  two 
dimensional  gas  of  pentadecylic  acid  water  molecules  on  the  surface  of  water. 

We  are  now  in  the  process  of  analysis  of  the  problem  of  the  accurate  characteri 
zation  of  the  laminar-turbulent  rlow  transition  in  high  speed  flow.  The  exper¬ 
imental  method  employed  would  be  that  of  optical  mixing  spectroscopy,  or  if  the 
characteristic  size  of  fluctuation  requires  it,  we  will  use  microwave  or  ultra¬ 
sonic  waves  for  the  scattering  spectroscopy. 


